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Parti 

Preface 



In this working group we have investigated a number of aspects of searches for new physics beyond the Standard 
Model (SM) at the running or planned TeV-scale colliders. For the most part, we have considered hadron 
colliders, as they will define particle physics at the energy frontier for the next ten years at least. The variety 
of models for Beyond the Standard Model (BSM) physics has grown immensely. It is clear that only future 
experiments can provide the needed direction to clarify the correct theory. Thus, our focus has been on exploring 
the extent to which hadron colliders can discover and study BSM physics in various models. We have placed 
special emphasis on scenarios in which the new signal might be difficult to find or of a very unexpected nature. 
For example, in the context of supersymmetry (SUSY), we have considered: 

• how to make fully precise predictions for the Higgs bosons as well as the superparticles of the Minimal 



Supersymmetric Standard Model (MSSM) (parts ^ and |IV|); 

MSSM scenarios in which most or all SUSY particles have rather large masses (parts ^and |V^ ); 

the ability to sort out the many parameters of the MSSM using a variety of signals and study channels 



(part|VII[); 

whether the no-lose theorem for MSSM Higgs discovery can be extended to the next-to-minimal Super- 
symmetric Standard Model (NMSSM) in which an additional singlet superfield is added to the minimal 
collection of superfields, potentially providing a natural explanation of the electroweak value of the pa- 



rameter n (part VIII); 



sorting out the effects of CP violation using Higgs plus squark associate production (part|IX[); 
the impact of lepton flavor violation of various kinds (part 



• experimental possibilities for the gravitino and its sgoldstino partner (part pC| ); 

• what the implications for SUSY would be if the NuTeV signal for di-muon events were interpreted as a 
sign of R-parity violation (part |XII| ). 

Our other main focus was on the phenomenological implications of extra dimensions. There, we considered: 



constraints on Kaluza Klein (KK) excitations of the SM gauge bosons from existing data (part XIII) and 



the corresponding projected LHC reach (part |XIVD ; 

techniques for discovering and studying the radion field which is generic in most extra-dimensional 



scenarios (part XV); 

the impact of mixing between the radion and the Higgs sector, a fully generic possibility in extra- 



dimensional models (part XVI); 



production rates and signatures of universal extra dimensions at hadron colliders (part pCVII[ ); 



• black hole production at hadron colliders, which would lead to truly spectacular events (part XVIII ). 
The above contributions represent a tremendous amount of work on the part of the individuals involved and 
represent the state of the art for many of the currently most important phenomenological research avenues. Of 
course, much more remains to be done. For example, one should continue to work on assessing the extent to 
which the discovery reach will be extended if one goes beyond the LHC to the super-high-luminosity LHC 
(SLHC) or to a very large hadron collider (VLHC) with ^/s ~ 40 TeV Overall, we believe our work shows 
that the LHC and future hadronic colliders will play a pivotal role in the discovery and study of any kind 
of new physics beyond the Standard Model. They provide tremendous potential for incredibly exciting new 
discoveries. 
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Part II 

Theoretical Developments 

/. Gunion, J. Hewett, K. Matchev, T. Rizzo 

Abstract 

Various theoretical aspects of physics beyond the Standard Model at hadron colliders 
are discussed. Our focus will be on those issues that most immediately impact the 
projects pursued as part of the BSM group at this meeting. 

1. Introduction 

The Standard Model (SM) has had a tremendous success describing physical phenomena up to energies ~ 100 
GeV. Yet some of the deep questions of particle physics are still shrouded in mystery - the origin of electroweak 
symmetry breaking (and the related hierarchy problem), the physics of flavor and flavor mixing, CP-violation 
etc. Any attempt to make further theoretical progress on any one of these issues necessarily requires new 
physics beyond the SM. 

It is generally believed that the TeV scale will reveal at least some of this new physics. Throughout 
history, we have never gone a whole order of magnitude up in energy without seeing some new phenomenon. 
Further support is given by attempts to solve the gauge hierarchy problem. Either there is no Higgs boson 
in the SM and then some new physics must appear around the TeV scale to unitarize WW scattering, or the 
Higgs boson exists, and one has to struggle to explain the fact that its mass is minute in (fundamental) Planck 
mass units. Very roughly, there are three particularly compelling categories of new physics that are capable of 
solving the hierarchy problem. 

• Supersymmetry (SUSY): 

Low energy supersymmetry eliminates the quadratic ultraviolet sensitivity of the Higgs boson mass, 
which arises through radiative corrections. Supersymmetry guarantees that these contributions cancel 
between loops with particles and those with their superpartners, making the weak scale natural provided 
the superpartner masses are 0(1 TeV). 

In its minimal version, a supersymmetrized standard model has only one additional free parameter - the 
supersymmetric Higgs mass /u. However, supersymmetry has to be broken, which leads to a proliferation 
of the number of independent input parameters. There are many different models on the market, differing 
only in the way SUSY breaking is communicated to "our world". Furthermore, one can go beyond the 
minimal supersymmetric extension of the Standard Model (MSSM), e.g. to the Next-to-Minimal Super- 
symmetric Standard Model (NMSSM) where an extra singlet superfield is added to the MSSM matter 
content. Then the so-called R-parity breaking models introduce additional Yukawa-type couplings be- 
tween the SM fermions and their superpartners; there are models with multiple extra U(l) gauge groups, 
etc. (for a recent review, see [|I]]). Garden varieties of all of these models have been extensively studied. 
In this report, our focus will be on models which yield unusual signatures and/or make discovery/study 
of SUSY more difficult. 

• Technicolor (TC): 

Technicolor (for a recent review, see ||^) has made a resurgence through models where the heavy top 
quark plays an essential role, such as the top-color assisted technicolor model and models in which 
an extra heavy singlet quark joins with the top-quark to give rise to electroweak symmetry breaking 
(EWSB). Very little work was done on this class of models at this workshop and so we will not discuss 
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such models further. It should, however, be noted that in most of these models, an effective low-energy 
Higgs sector emerges that typically is equivalent to a general two-Higgs-doublet model (2HDM). Light 
pseudo-Nambu-Goldstone bosons can also be present. 
• Extra dimensions: 

Extra dimensions at or near the TeV^^ scale may bring the relevant fundamental particle physics scale 
down to a TeV and thus eliminate the hierarchy problem [§,0. If this scenario were true, it would have a 
profound influence on all types of physics at the LHC and other future colliders. Extra dimensions impact 
the Higgs sector and can even give rise to EWSB. They can also lead to Kaluza Klein (KK) excitations of 
normal matter. The production of small black holes at the LHC becomes a possibility. Such black holes 
would promptly decay to multiple SM particles with a thermal distribution, giving striking signatures. A 
number of the many possibilities and the related experimental consequences were explored during this 
workshop and are reported here. 



2. SUSY and expectations for hadron colliders 

Even within the context of the minimal supersymmetric model (MSSM) with R-parity conservation, there are 
103 parameters beyond the usual Standard Model (SM) parameters. Different theoretical ideas for soft-SUSY 
breaking can be used to motivate relations between these parameters, but as time progresses more and more 
models are being proposed. In addition, one cannot rule out the possibility that several sources of soft-SUSY 
breaking are present simultaneously. 

Typically, any theoretical model will provide predictions for the soft-SUSY breaking parameters at a 
high scale, such as the GUT scale. For example, in mSUGRA, the minimal supergravity model (sometimes also 
called the constrained MSSM - cMSSM), the universal GUT-scale scalar mass Mq, the universal GUT-scale 
gaugino mass A/1/2, the universal trilinear term Aq, the low-energy ratio tan /3 of Higgs vacuum expectation 
values, and the sign of the fi parameter, 

Mo,Mi/2,Ao,tan/3,sign(/u) (1) 

fully specify all the soft-SUSY breaking parameters once the renormalization group equations (RGE) are re- 
quired to yield correct EWSB. More generally, the RGEs provide a link between the experimentally observed 
parameters at the TeV scale and the fundamental physics at the high-energy scale. The amount of information 
we can extract from experiment is therefore related to the precision with which we can relate the values of the 
parameters at these two vastly different scales. Precise predictions require multi-loop results for the RGE and 
the related threshold corrections, and a careful assessment of all systematic uncertainties. This is the focus 



of a couple of the contributions to this report (parts III and IV). At the meeting, there was also considerable 



discussion of the extent to which a given set of low energy parameters could be ruled out or at least discrim- 
inated against by virtue of constraints such as: requiring that the LSP be the primary dark matter constitute; 
correct b —>■ S7; 'correct' gi^j — 2; etc. Currently there are many programs available for evaluating the impact 
of such constraints, and they tend to give diverse answers. In some cases, numerically important effects have 
been left out, e.g. certain co-annihilation channels, higher-order terms in the RGE equations, and so forth. In the 
remaining cases, the spread can be taken as an indication of the theoretical uncertainty involved in relating the 
TeV and unification scales. While progress in this area has been made, as summarized in [^, no summary of 
the status was prepared for this report. However, one important conclusion from this effort is clear. There are 
regions of parameter space, even for the conventional mSUGRA case of Eq. ([I]), for which very high sparticle 
masses could remain consistent with all constraints. This observation led to renewed focus on LHC sensitivity 
to SUSY models with very high mass scales (parts ^ and VI), as possibly also preferred by coupling constant 



unification with as{mz) < 0.12. For example, naturally heavy squark masses are allowed in the focus point 
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scenario [ 6|] and would ameliorate any possible problems with flavor-changing neutral currents (FCNC) related 
thereto [^]. 

More generally, it would be unwise for the experimental community to take too seriously the predictions 
of any one theoretical model for soft-SUSY breaking. It is important that convincing arguments be made that 
TeV-scale SUSY (as needed to solve the hierarchy problem) can be discovered for all possible models. Much 
work has been done in recent years in this respect, and such efforts were continued during the workshop and 
are reported on here. In general, the conclusions are positive; TeV-scale SUSY discovery at the LHC will be 
possible for a large class of models. Further, after the initial discovery, a multi-channel approach, like the one 
presented in part VII, can be used to determine the soft-SUSY-breaking parameters with considerable precision. 

An important aspect of verifying the nature of the SUSY model will be a full delineation of its Higgs 
sector. In the MSSM, the Higgs sector is a strongly constrained 2HDM. In particular, in the MSSM, there is 
a strong upper bound on the mass of the lightest CP-even Higgs boson {mh < 130 GeV) and strong relations 
between its couplings and the CP-odd Higgs mass parameter niA- As a result, there is a 'no-lose' theorem for 
MSSM Higgs discovery at the LHC (assuming that Higgs decays to pairs of SUSY particles are not spread 
out over too many distinct channels). However, if tua > 300 GeV and tan (3 has a moderate value somewhat 
above 3, then existing analyses indicate that it will be very hard to detect any Higgs boson other than the light 
CP-even h (which will be quite SM-lrke). The H, A and (all of which will have similar mass) might well 
not be observable at the LHC. Further work on extending the high-tan /3 r signals for the H, A, to the 
lowest possible tan (3 values and on finding new signals for them should be pursued. 

However, an even bigger concern is the additional difficulties associated with Higgs discovery if the 
MSSM is extended to include one or more additional singlet superfields (leading to additional Higgs singlet 
scalar fields). The motivation for such an extension is substantial. First, such singlets are very typical of string 
models. Second, it is well-known that there is no convincing source for a weak-scale value of the fj, parameter 
of the MSSM. The simplest and a very attractive model for generating a weak-scale value for /i is the NMSSM 
in which one singlet superfield is added to the MSSM. The superpotential term XSH^Hu (where S is the singlet 
superfield and Hd, Hu ai"e the Higgs superfields whose neutral scalar component vevs give rise to the down and 
up quark masses, respectively) gives rise to a weak-scale value for /i provided A is in the perturbative domain 
and (5) = O {mz)- Both of these conditions can be naturally implemented in the NMSSM. This simple and 
highly-motivated extension of the MSSM leads to many new features for SUSY phenomenology at the LHC 
and other future colliders. However, its most dramatic impact is the greatly increased difficulty of guaranteeing 
the discovery of at least one of the NMSSM Higgs bosons (there now being 3 CP-even Higgs bosons, 2 CP-odd 
Higgs bosons and a charged Higgs pair). Very substantial progress was made as part of this workshop in filling 
previously identified gaps in parameter space for which discovery could not be guaranteed. However, remaining 
additional dangerous parameter regions, and the new relevant experimental discovery channels, were identified. 
Substantial additional effort on the part of the LHC community will be required in order to demonstrate that 
Higgs discovery in these new channels will always be possible. Part VIII| of this report discusses these issues 
in some depth. 

In the simplest models of soft-SUSY-breaking, it is generally assumed that the soft-SUSY-breaking pa- 
rameters will not have phases (that cannot be removed by simple field redefinitions). Even in the MSSM, the 
presence of such phases would be an essential complication for LHC SUSY phenomenology, and most par- 
ticularly for Higgs sector discovery and study. In general, many things become more difficult. An exception 
would be if one can simultaneously produce a pair of squarks in association with a Higgs boson. Such signals 
would allow a first determination of the non-trivial phases of the theory, since the production of the CP-odd A 
in association with two light top squarks, A + ti + ti, is an unequivocal signal of non-trivial phases for the /i 
and A (soft tri-linear) parameters of the MSSM. Some aspects of this are explored in part LX. The experimental 
viability of such signals will require further study. 
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In many SUSY models, lepton flavor violating (LFV) decays of various particles can occur. Lepton- 
flavor- violating interactions can easily arise as a result of a difference between the flavor diagonalization in the 
normal fermionic leptonic sector as compared to that in the slepton sector. Typically this is avoided by one of 
two assumptions: a) a common leptonic flavor structure for the lepton and slepton sectors (alignment) or b) 
flavor-blind mechanism of SUSY breaking, which yields slepton mass matrices which are diagonal in flavor 
space. No convincing GUT-scale motivation for either of these possibilities has been expounded. In fact, many 
string models suggest quite the contrary (see, e.g. [||]). Further, neutrino masses and mixing phenomenology 
could be indicating the presence of lepton flavor violating interactions, especially in the context of the see-saw 
mechanism. In particular, as shown in part ^ of this report, expectations based on neutrino mixing phenomenol- 
ogy lead to rates for t ^ fij decays at high tan (3 (which enhances these decays in the MSSM) that are very 
similar to existing bounds on such decays, implying that they might be observed in the next round of exper- 
iments. If one wishes to suppress LFV decays in the most general case, very large slepton masses would be 
required. This would, of course, fit together with the large squark masses needed for guaranteed suppression of 
FCNC decays. 

One parameter that is not conventionally included in the 103 MSSM SUSY parameters is the goldstino 
mass (which determines the mass of the spin-3/2 gravitino). The gravitino mass is related to the scale of SUSY 
breaking F by 




Further, the interactions of the goldstino part of the gravitino (and of its spin zero sgoldstino partners) are 
proportional to l/F. (The masses of the goldstinos ms,mp are not determined.) In mSUGRA models and the 
Mke, F is sufficiently large that the goldstino and sgoldstino masses are so large, and their interaction strengths 
so small, that they are not phenomenologically relevant. However, in some models of SUSY breaking F is 
relatively small. A well-known example is gauge-mediated SUSY breaking for which F can be small enough 
for the goldstino to be the true LSP into which all more massive SUSY particles ultimately decay. In such 
a case, all of SUSY phenomenology changes dramatically. The sgoldstinos might also be light, with masses 
anywhere below 1 TeV being reasonable. In this case, for \/F < 1 TeV, they could yield some very significant 



experimental signals, discussed in part XI . For example, they might appear in rare decays of the J/ ip and T 



or lead to FCNC interactions. For small enough F, direct production of sgoldstinos becomes significant at the 
LHC for masses up to about a TeV (in particular via a gg ^ S vertex of the form -^^F^J^F^^^°^S) and would 
yield some unique signatures. 

The possibility of R-parity violation in SUSY models has been extensively considered There are 
three possible sets of RPV couplings as specified in the superpotential: 

KjkLiLjEk + X'iji^LiQjDk + X'ljjJJiDjDk , (3) 

where SU(2) and color-singlet structures are implied. Here, Xijk iX'ljf^) must be antisymmetric under i ^ j 
(j ^ k). For proton stability, we require that either the A'/^^ = or that Xijk = X!-^ = 0. One of the most 
under-explored possibilities for the LHC is that one or more of the A'"s is non-zero. This would imply that 
the neutialino ultimately decays to 3 jets inside the detector. There would be no missing energy. If the mass 
difference between the x? and xt is small (as possible, for example, for anomaly mediated SUSY-breaking and 
in some types of string-motivated boundary conditions) or if the leptonic branching fractions of the charginos 
and heavier neutralinos are small, then there might also be few hard leptons in the LHC events. The main 
SUSY signature would be extra events with large numbers of jets. Whether or not such events can be reliably 
extracted from the large QCD background, and especially the maximum SUSY particle mass for which such 
extraction is possible, is a topic awaiting future study. The leptonic type of RPV would lead to very clear LHC 
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signals for SUSY, in which events would contain extra leptons as well as some missing energy from the extra 
neutrinos that would emerge from decays. For example, A212 would lead to decays of the neutralino LSP such 
as Xi At/"i^- 



It is just possible that the NuTeV dilepton events [ |10| ] could be a first sign of R-parity violation. The 



explanation proposed in part pClI| requires A232 7^ (leading to the decays Xi /^l^b^t ^^'^ ^1 ~^ ^^R^l 



and conjugates thereof). The explanation proposed for the Tevatron events, in which the light neutralinos are 
produced in B^, decays) would also require the existence of a mixed leptonic-hadronic RPV coupling A'^]^3. 
In general, the weakness of the constraints on couplings involving the 3rd generation and the large size of the 
similar Yukawa couplings related to quark mass generation both favor signals related to 3rd generation leptons 
and quarks. 

3. Extra Dimensions 

An alternative to SUSY for explaining the hierarchy problem is that the geometry of space-time is modified at 
scales much less than the Planck scale, Mpi. In such models, which may still be regarded as rather speculative, 
but have attracted a lot of attention recently, the 3-spatial dimensions in which we live form a 3-dimensional 
'membrane', called 'the wall', embedded in a much larger extra dimensional space, known as 'the bulk', and 
that the hierarchy between the weak scale ~ 10^ GeV and the 4-dimensional Planck scale M pi ~ 10^^ GeV is 
generated by the geometry of the additional bulk dimensions. This is achievable either by compactifying all the 
extra dimensions on tori, or by using strong curvature effects in the extra dimensions. In the first case, Arkani- 
Hamed, Dimopoulos, and Dvali (ADD) [||,^, 12| used this picture to generate the hierarchy by postulating a 



large volume for the extra dimensional space. In the latter case, the hierarchy can be established by a large 
curvature of the extra dimensions as demonstrated by Randall and Sundrum (RS) [Q]. It is the relation of 
these models to the hierarchy which yields testable predictions at the TeV scale. Such ideas have led to extra 
dimensional theories which have verifiable consequences at present and future colliders. 

There are three principal scenarios with predictions at the TeV scale, each of which has a distinct phe- 
nomenology. In theories with Large Extra Dimensions, proposed by ADD [^, II, T^, gravity alone propagates 



in the bulk where it is assumed to become strong near the weak scale. Gauss' Law relates the (reduced) Planck 
scale Mpi of the effective 4d low-energy theory and the fundamental scale Mp), through the volume of the 
6 compactified dimensions, Vs, via Mpi = VsMp, . Mpi is thus no longer a fundamental scale as it is 
generated by the large volume of the higher dimensional space. If it is assumed that the extra dimensions are 
toroidal, then setting Mp) ^ TeV to eliminate the hierarchy between M pi and the weak scale determines the 
compactification radius R of the extra dimensions. Under the further simplifying assumption that all radii are 
of equal size, Vs = {2ttRY, R then ranges from a sub-millimeter to a few fermi for 5 = 2 — 6. Note that the 
case of (5 = 1 is excluded as the corresponding dimension would directly alter Newton's law on solar-system 
scales. The bulk gravitons expand into a Kaluza-Klein (KK) tower of states, with the mass of each excitation 
state being given by = v? /R^. With such large values of R the KK mass spectrum appears almost contin- 
uous at collider energies. The ADD model has two important collider signatures: (z) the emission of real KK 
gravitons in a collision process leading to a final state with missing energy and {ii) the exchange of virtual KK 
graviton towers between SM fields which leads to effective dim-8 contact interactions. Except for the issue of 
Black Hole (BH) production to be discussed below, we will say no more about the ADD scenario as work was 
not performed on this model at this workshop. 

A second possibility is that of Warped Extra Dimensions; in the simplest form of this scenario [||] gravity 
propagates in a 5d bulk of finite extent between two (3 + 1) -dimensional branes which have opposite tensions. 
The Standard Model fields are assumed to be constrained to one of these branes which is called the TeV brane. 
Gravity is localized on the opposite brane which is referred to as the Planck brane. This configuration arises 
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from the metric ds^ = e~'^^yr]^ydx^dx^ — dy^ where the exponential function, or warp factor, multiplying the 
usual 4d Minkowski term produces a non-factorizable geometry, and y G [0, vri?] is the coordinate of the extra 
dimension. The Planck (TeV) brane is placed at y = 0(7ri?). The space between the two branes is thus a slice 
of AdS^: 5d anti-deSitter space. The original extra dimension is compactified on a circle so that the wave 
functions in the extra dimension are periodic and then orbifolded by a single discrete symmetry Z2 forcing the 
KK graviton states to be even or odd under y —y. Here, the parameter k describes the curvature scale, which 
together with M£, {D = 5) is assumed [Q] to be of order Mpi, with the relation Mp; = M^^/k following 
from the integration over the 5d action. Note that that there are no hierarchies amongst these mass parameters. 
Consistency of the low-energy description requires that the 5d curvature, = — 20fc^, be small in magnitude 
in comparison to Md, which implies k/Mpi <, 0.1. We note that mass scales which are naturally of order 
M PI on the y = brane will appear to be of order the TeV scale on the y = ttR brane due to the exponential 
warping provided that ttR ~ 11 — 12. This leads to a solution of the hierarchy problem. 

The 4d phenomenology of the RS model is governed by two parameters, A^^ = M pie^^^'^ , which 
is of order a TeV, and k/Mpi. The masses of the bulk graviton KK tower states are m„ = Xnke~^^'^ = 
Xn^nk/Mpi with the Xn being the roots of the first-order Bessel function Ji. The KK states are thus not 
evenly spaced. For typical values of the parameters, the mass of the first graviton KK excitation is of order a 



TeV. The interactions of the bulk graviton KK tower with the SM fields are [ ]13[ ] 

1 1 °° 

Mpi ^ ^ M^^^ ^ 



AC = -^r^-(x)/.W(x) - —T>^^{x) 5^ h(;J{x) , (4) 



n=l 



where T^^'^ is the stress-energy tensor of the SM fields, /i^fj is the ordinary graviton and h^^J are the KK graviton 
tower fields. Experiment can determine or constrain the masses m„ and the coupling A^. In this model KK 
graviton resonances with spin-2 can be produced in a number of different reactions at colliders. Extensions of 



this basic model allow for the SM fields to propagate in the bulk []14|-|18|]. In this case, the masses of the bulk 
fermion, gauge, and graviton KK states are related. A third parameter, associated with the fermion bulk mass, 
is introduced and governs the 4d phenomenology. In this case, KK excitations of the SM fields may also be 
produced at colliders. 

One important aspect of the RS model is the need to stabilize the separation of the two branes with 



kR ~ 11 — 12 in order to solve the hierarchy problem. This can be done in a natural manner [19] but leads 
to the existence of a new, relatively light scalar field with a mass significantly less than A^ called the radion. 
This is most likely the lightest new state in the RS scenario. The radion has a flat wavefunction in the bulk 
and is a remnant of orbifolding and of the graviton KK decomposition. This field couples to the trace of the 
stress-energy tensor, ~ T^/A^, and is thus Higgs-like in its interactions with SM fields. In addition, it may 
mix with the SM Higgs altering the couplings of both fields. Searches for the radion and its influence on the 
SM Higgs couplings will be discussed below. 

The possibility of TeV~^-sized extra dimensions arises in braneworld models [20|. By themselves, they 
do not allow for a reformulation of the hierarchy problem but they may be incorporated into a larger structure 
in which this problem is solved. In these scenarios, the Standard Model fields may propagate in the bulk. This 
allows for a wide number of model building choices: 

• all, or only some, of the SM gauge fields are present in the bulk; 

• the Higgs field(s) may be in the bulk or on the brane; 

• the SM fermions may be confined to the brane or to specific locales in the extra dimensions. 

If the Higgs field(s) propagate in the bulk, the vacuum expectation value (vev) of the Higgs zero-mode, the 
lowest lying KK state, generates spontaneous symmetry breaking. In this case, the gauge boson KK mass 
matrix is diagonal with the excitation masses given by [Mq + n ■ n/i?^]^/^, where Mq is the vev-induced 
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mass of the gauge zero-mode and n labels the KK excitations in 5 extra dimensions. However, if the Higgs 
is confined to the brane, its vev induces off-diagonal elements in the mass matrix generating mixing amongst 
the gauge KK states of order (MqR)'^. For the case of 1 extra dimension, the coupling strength of the bulk 
KK gauge states to the SM fermions on the brane is V^g, where g is the corresponding SM coupling. The 
fermion fields may (a) be constrained to the (3 + 1) -brane, in which case they are not directly affected by the 
extra dimensions; (b) be localized at specific points in the TeV^^ dimension, but not on a rigid brane. Here the 
zero and excited mode KK fermions obtain narrow Gaussian-like wave functions in the extra dimensions with 
a width much smaller than R^^. This possibility may suppress the rates for a number of dangerous processes 



such as proton decay [gJJ. (c) The SM fields may also propagate in the bulk. This scenario is known as 
universal extra dimensions [^. (4 + 5) -dimensional momentum is then conserved at tree-level, and KK parity, 
(— 1)", is conserved to all orders. TeV extra dimensions lead to an array of collider signatures some of which 
will be discussed in detail below. 

Theories with extra dimensions and a low effective Planck scale (Mo) offer the exciting possibility 
that black holes (BH) somewhat more massive than Mo can be produced with large rates at future colliders. 
Cross sections of order 100 pb at the LHC have been advertised in the analyses presented by Giddings and 



Thomas [ |23| ] and by Dimopoulos and Landsberg [Q]. These early analyses and discussions of the production 
of BH at colliders have been elaborated upon by several groups of authors [ 25 ] and the production of BH by 



cosmic rays has also been considered [[32|-|39|]. A most important question to address is whether or not the BH 
cross sections are actually this large or, at the very least, large enough to lead to visible rates at future colliders. 

The basic idea behind the original collider BH papers is as follows: consider the collision of two high 
energy SM partons which are confined to a 3-brane, as they are in both the ADD and RS models. In addition, 
gravity is free to propagate in 6 extra dimensions with the 4+6 dimensional Planck scale assumed to he Md ~ 1 
TeV. The curvature of the space is assumed to be small compared to the energy scales involved in the collision 
process so that quantum gravity effects can be neglected. When these partons have a center of mass energy in 
excess of ~ Af and the impact parameter of the collision is less than the Schwarzschild radius, Rs, associated 
with this center of mass energy, a 4+(5-dimensional BH is formed with reasonably high efficiency. It is expected 
that a very large fraction of the collision energy goes into the BH formation process so that Mbh — 
The subprocess cross section for the production of a non-spinning BH is thus essentially geometric for each 
pair of initial partons: a ~ e7ri?|, where e is a factor that accounts for finite impact parameter and angular 
momentum corrections and is expected to be <^ 1. Note that the 4 + (5-dimensional Schwarzschild radius 



scales as Rs 



Mbh 



, apart from an overall 6- and convention-dependent numerical prefactor. This 



approximate geometric subprocess cross section expression is claimed to hold when the ratio M bh I^'Id is 
"large", i.e., when the system can be treated semi-classically and quantum gravitational effects are small. 



Voloshin [ |40| , |41p has provided several arguments which suggest that an additional exponential suppres- 
sion factor must be included which presumably damps the pure geometric cross section for this process even in 
the semi-classical case. This issue remains somewhat controversial. Fortunately it has been shown [ p^ that the 
numerical influence of this suppression, if present, is not so great as to preclude BH production at significant 
rates at the LHC. These objects will decay promptly and yield spectacular signatures. A discussion of BH 
production at future colliders is presented in one of the contributions. 
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Part III 

FeynSSG v.1.0: Numerical Calculation of the 
mSUGRA and Higgs spectrum 

A. Dedes, S. Heinemeyer, G. Weiglein 

Abstract 

FeynSSG v . 1 . is a program for the numerical evaluation of the Supersymmet- 
ric (SUSY) particle spectrum and Higgs boson masses in the Minimal Supergravity 
(mSUGRA) scenario. We briefly present the physics behind the program and as an 
example we calculate the SUSY and Higgs spectrum for a set of sample points. 

In the Minimal Supersymmetric Standard Model (MSSM) no specific assumptions are made about the 
underlying SUSY-breaking mechanism, and a parameterization of all possible soft SUSY-breaking terms is 
used. This gives rise to the huge number of more than 100 new parameters in addition to the SM, which 
in principle can be chosen independently of each other. A phenomenological analysis of this model in full 
generality would clearly be very involved, and one usually restricts to certain benchmark scenarios, see Ref. [Q] 
for a detailed discussion. On the other hand, models in which all the low-energy parameters are determined in 
terms of a few parameters at the Grand Unification (GUT) scale (or another high-energy scale), employing a 
specific soft SUSY-breaking scenario, are much more predictive. The most prominent scenario at present is the 
minimal Supergravity (mSUGRA) scenario [p3|-52]. 



In this note we present the Fortran code FeynSSG for the evaluation of the low-energy mSUGRA spec- 
trum, including a precise evaluation for the MSSM Higgs sector. The high-energy input parameters (see below) 
are related to the low-energy SUSY parameters via renormalization group (RG) running (taken from the pro- 



gram SUITY [53, 54]), taking into account contributions up to two-loop order. The low-energy parameters are 



then used as input for the program FeynHiggs [ [55| ] for the evaluation of the MSSM Higgs sector. 



The simplest possible choice for an underlying theory is to take at the GUT scale all scalar particle masses 
equal to a common mass parameter Mq, all gaugino masses are chosen to be equal to the parameter and 
all trilinear couplings flavor blind and equal to Aq. This situation can be arranged in Gravity Mediating SUSY 
breaking Models by imposing an appropriate symmetry in the Kahler potential [p^-|52|], called the minimal 
Supergravity (mSUGRA) scenario. In order to solve the minimization conditions of the Higgs potential, i.e. in 
order to impose the constraint of REWSB, one needs as input tan [3{Mz) and sign(/i). The running soft SUSY- 
breaking parameters in the Higgs potential, ruHi ^^id > are defined at the EW scale after their evolution 
from the GUT scale where we assume that they also have the common value Mq. Thus, apart from the SM 
parameters (determined by experiment) 4 parameters and a sign are required to define the mSUGRA scenario: 

{ Mo , Ml/2 , Ao , tan/3 , sign(/i) } . (1) 

In the numerical procedure we employ a two-loop renormalization group running for all parameters involved, 
i.e. all couplings, dimensionful parameters and VEV's. We start with the MS values for the gauge couplings 
at the scale Mz, where for the strong coupling constant a trial input value in the vicinity of 0.120 is used. 
The MS values are converted into the corresponding DR ones ||5^. The MS running h and r masses are run 
down to mb = 4.9 GeV, = 1.777 GeV with the SU(3)(^ x U(l)g,„ RGB's [|^] to derive the running 
bottom and tau masses (extracted from their pole masses). This procedure includes all SUSY corrections at the 
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one-loop level and all QCD corrections at the two-loop level as given in [58]. Afterwards by making use of the 
two-loop RGB's for the running masses m^, nir, we run upwards to derive their MS values at Mz, which are 
subsequently converted to the corresponding DR values. This procedure provides the bottom and tau Yukawa 
couplings at the scale Mz- The top Yukawa coupling is derived from the top-quark pole mass, = 175 GeV, 
which is subsequently converted to the DR value, rrTi{mt), where the top Yukawa coupling is defined. The 
evolution of all couplings from Mz running upwards to high energies now determines the unification scale 
Mqut and the value of the unification coupling acuT by 

ai(MGUT)lDR = «2(^^gut)Idr = «GUT • (2) 

At the GUT scale we set the boundary conditions for the soft SUSY breaking parameters, i.e. the values for Mq, 
Ml/2 and ^0 are chosen, and also a3(MGUT) is set equal to acuT- All parameters are run down again from 
Mqut to Mz- For the calculation of the soft SUSY-breaking masses at the EW scale we use the "step function 
approximation" [53, 54|]. Thus, if the equation employed is the RGE for a particular running mass m{Q), then 



Qq is the corresponding physical mass determined by the condition m((5o) = Qo- After running down to Mz, 
the trial input value for as has changed. At this point the value for tan /3 is chosen and fixed. The parameters 
ImI and B are calculated from the minimization conditions 



tan2/3-l 2' 



f^\Q) = H \ - t^mUQ) , (3) 



(mi(Q)^ + m2(Q)^)sin2/3(Q) 

Only the sign of the /i-parameter is not automatically fixed and thus chosen now. This procedure is iterated 
several times until convergence is reached. 

In Q is the renormalization scale. It is chosen such that radiative corrections to the effective 

potential are rather small compared to other scales. In tan/3 = V2/V1 is the ratio of the two vacuum 

expectation values of the Higgs fields H2 and Hi responsible for giving masses to the up-type and down-type 
quarks, respectively. In @,(^, tan /3 is evaluated at the scale Q, from the scale Mz, where it is considered as 
an input parameter]]. By mjj, = m'jj_ + S^- in (|3|),(^ we denote the radiatively corrected "running " Higgs 
soft-SUSY breaking masses and 

2 _ ^2 I ,,2 i ip — ^2 I ,,2 



m, = m 



.+^' + S,^ = mjj^ {i = 1, 2) , (5) 



where E^,; are the one-loop corrections based on the 1-loop Coleman- Weinberg effective potential AV, S^. 

1 dAV 
2vi dvi ' 



2 



M 
Q2 



(6) 



Here .J a is the spin of the particle a, Ca are the color degrees of freedom, and = 1(2) for real scalar 
(complex scalar), = 1(2) for Majorana (Dirac) fermions. Q is the energy scale and the Ma are the field 
dependent mass matrices. Explicit formulas of the S^,- are given in the Appendices of [58,^]. In our analyses 



contributions from all SUSY particles at the one-loop level are incorporated^. With M| here we denote the 



' See for example the discussion in the Appendix of [ 



^ The corresponding two-loop corrections are now available for a general renormalizable softly broken SUSY theory [^]. Assuming 
the size of these higher-order corrections to be of the same size as for the Higgs-boson mass matrix, the resulting values of ^ and B 
could change by ~ 5 — 10%. The possible changes would hardly affect the results in the Higgs-boson sector but could affect to some 
extent the analysis of SUSY particle spectra, especially when Mo and M1/2 are lying in different mass regions. 
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tree level "running" Z boson mass, M|(Q) = \{g\ ^ si)^^ (^^ = f f + v^, extracted at the scale Q from 



its physical pole mass Mz 
conditions (|),(|)[} 



91.187 GeV. The REWSB is fulfilled, if and only if there is a solution to the 



For the predictions in the MSSM Higgs sector we use the code FeynHiggs []55|], which is implemented 
as a subroutine into FeynSSG. The code is based on the evaluation of the low-energy Higgs sector parameters 
in the Feynman-diagrammatic (FD) approach [^-^] within the on-shell renormalization scheme. Details 
about the conversion of the low-energy results from the RG running, obtained in the DR scheme, to the on- 
shell scheme can be found in Ref. [^5|]. In the FD approach the masses of the two CP-even Higgs bosons, 
m/i and vtlh, are derived beyond tree level by determining the poles of the h — -propagator matrix, which is 
equivalent to solving the equation 



m 



H,tme 



0, 



(V) 



where T^s, s = h, H, hH denotes the renormalized Higgs boson self-energies. Their evaluation consists of the 
complete one-loop result combined with the dominant two-loop contributions of 0{aas) []62|-^ and further 



subdominant corrections [56, 57|, see Refs. [62-64, 68] for details 



An analysis employing FeynSSG for the constraints on the mSUGRA scenario from the Higgs boson 
search at LEP2 and the corresponding implications for SUSY searches at future colliders has been presented 



in Refs. []65ip9|]. As another example we present here the results of the low-energy SUSY spectrum for some 
sample points [[70|]. (Some of these sample points are now included in the "SPS" (Snowmass Points and 
Slopes) [^ that have recently been proposed as new benchmark scenarios for SUSY searches at current and 
future colliders.) 

The sample points are presented in Table 1. For these results we have set the 1-loop corrections S^,^ 
equal to zero and all the thresholds are switched on. Thus for the points considered here a one loop improved 
tree level analysis is done. If we switch on the full 1-loop corrections then the points E,F,H,J,K, and M, 
fail to satisfy electroweak symmetry breaking, /i^ from (||) is negative. In addition, the weak mixing angle, 
sin^ Owi^z), has been set to 0.2315. An updated version which employs the effective weak mixing angle as 
a boundary condition at the electroweak scale is under way (in fact such an analysis had been done in the past 
using the program SUITY, see [^, 54, 7l|]). It is intended to regularly update FeynSSG with the upcoming 
new versions of the SUITY and FeynHiggs programs. 



^ Sometimes in the literature, the requirement of the REWSB is described by the inequality m\{Q)m\{Q) — \^{Q)B{Q)\^ < 0. 
This relation is automatically satisfied here from and from the fact that the physical squared Higgs masses must be positive. 
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Model 


A 


B 


c 


D 


E 


F 


G 


11 


1 


J 


K 


E 


M 


'"'1/2 


624 


258 


415 


549 


315 


1090 


390 


1585.5 


364 


785 


1006 


471 


1600 




137 


100 


90 


120 


1500 


2970 


123 


459 


188 


320 


1000 


330 


1500 


tan l3 


5 


10 


10 


10 


10 


10 


20 


20 


35 


35 


40.3 


45 


48 


sign(^) 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 


^0 









































III I 


175 


175 


175 


175 


175 


175 


175 


175 


175 


175 


175 


175 


175 


Masses 




























\ut\ 


811 


362 


551 


705 




941 


515 


1719 


480 


936 




595 


1660 




114 


113 


116 


116 




118 


117 


121 


117 


121 




119 


123 


HO 


947 


414 


629 


769 




3171 


580 


2065 


502 


1003 




578 


1709 


AO 


947 


414 


629 


769 




3171 


580 


2065 


502 


1003 




578 


1709 


H± 


939 


420 


625 


789 




3151 


569 


1920 


472 


867 




461 


818 


Y? 
Al 


260 


101 


169 


229 




475 


158 


693 


148 


332 




196 


705 




484 


185 


314 


429 




853 


295 


1273 


274 


618 




363 


1293 


y§ 

A3 


813 


368 


555 


707 




942 


520 


1720 


485 


938 




599 


1661 


y5 

A4 


827 


387 


570 


713 




985 


534 


1728 


499 


948 




611 


1670 


Al 


483 


185 


314 


429 




852 


295 


1273 


274 


618 




362 


1293 




826 


387 


570 


715 




985 


534 


1728 


500 


948 




612 


1670 


.9 


1382 


619 


953 


1228 




2371 


901 


3266 


847 


1713 




1074 


3301 


BL, Ml 


437 


206 


295 


386 




3038 


292 


1127 


311 


610 




456 


1818 


e/j, Mi? 


273 


146 


184 


241 




2991 


195 


744 


236 


435 




376 


1609 




431 


190 


284 


378 




3037 


281 


1125 


300 


605 




449 


1816 


Tl 


271 


137 


176 


234 




2966 


168 


702 


165 


351 




261 


1228 


T2 


438 


209 


297 


387 




3026 


299 


1118 


322 


602 




449 


1673 


Vr 


430 


189 


283 


377 




3025 


277 


1112 


289 


584 




419 


1666 


Ul, Cl 


1261 


575 


874 


1122 




3546 


831 


2958 


794 


1581 




1028 


3293 


Ur, Cr 


1216 


559 


845 


1082 




3507 


805 


2835 


770 


1524 




997 


3183 


dh, SL 


1264 


581 


877 


1125 




3547 


835 


2959 


798 


1583 




1031 


3294 


dR, Sr 


1211 


559 


843 


1078 




3503 


803 


2820 


768 


1517 




994 


3169 


tl 


971 


419 


663 


874 




2465 


630 


2340 


596 


1237 




779 


2534 




1211 


604 


864 


1076 




3077 


820 


2735 


772 


1457 




953 


2826 


bi 


1167 


531 


807 


1037 




3071 


754 


2711 


686 


1393 




859 


2739 


b2 


1211 


560 


842 


1075 




3481 


799 


2772 


752 


1460 




941 


2833 



Table 1: Mass spectra in GeV for mSUGRA points calculated with program FeynSSG vl . (see text for details). Points (E) and 
(K) fail to pass the Radiative Electroweak Breaking requirement, i.e., jj,^ < 0. Points (F) and (M) exhibit instabiUty, i.e., the program 
reaches a poor convergence. The charged Higgs Boson (i?*) mass is given at tree level. 
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Part IV 

Theoretical Uncertainties in Sparticle Mass 
Predictions and SOFTSUSY 

B. C. Allanach 

Abstract 

We briefly introduce the SOFTSUSY calculation of sparticle masses and mixings 
and illustrate the output with post-LEP benchmarks. We contrast the sparticle 
spectra obtained from ISASUGRA7.58, SUSPECT2.004 with those obtained from 
S0FTSUSY1.3 along SNOWMASS model lines in minimal supersymmetric stan- 
dard model (MSSM) parameter space. From this we gain an idea of the uncertainties 
involved with sparticle spectra calculations. 

Supersymmetric phenomenology is notoriously complicated. Even if one assumes the particle spectrum 
of the minimal supersymmetric standard model (MSSM), fundamental patterns of supersymmetry (SUSY) 
breaking are numerous. It seems that there is currently nothing to strongly favor one particular scenario above 



all others. In ref. [72], it was shown that measuring two ratios of sparticle masses to 1% could be enough to 
discriminate different SUSY breaking scenarios (in that case, mirage, grand-unified or intermediate scale type 
I string-inspired unification). Thus, in order to discriminate high energy models of supersymmetry breaking, 
it will be necessary to have better than 1% accuracy in both the experimental and theoretical determination of 



some superparticle masses. An alternative bottom-up approach [ l73[ ] is to evolve soft supersymmetry breaking 
parameters from the weak scale to a high scale once they are 'measured'. The parameters of the high-scale 
theory are then inferred, and theoretical errors involved in the calculation will need to be minimized. 



We now briefly introduce SOFTSUSYl . 3 ||74p, a tool to calculate the masses and mixings of MSSM 
sparticles. It can be downloaded from the URL 

http://allanach. home .cern.ch/allanach/softsusy.html 

It is vahd for the R-parity conserving MSSM with real couplings and includes full 3-family particle or sparticle 



mixing. The manual [74] can be consulted for a more complete description of approximations and the algo- 
rithm used. Low energy data (together with tan (5{Mz)) set the Standard Model gauge couplings and Yukawa 
couplings: Gp, a, as{Mz) and the fermion masses and CKM matrix elements. The user provides a high- 
energy unification scale and supersymmetry breaking boundary conditions at that scale. The program derives 
the MSSM spectrum consistent with both of these constraints and radiative electro weak symmetry breaking at a 
scale MsusY = y/f^HJ^H^- Below Mz, three-loop QCD(8)one-loop QED is used to evaluate the MS Yukawa 
couplings and gauge couplings at Mz- These are then converted into the DR scheme, including finite and 
logarithmic corrections coming from sparticle loops. All one-loop corrections are added to the top mass and 
gauge couplings, while the other Standard Model couplings receive approximations to the full one-loop result. 
The radiative electroweak symmetry breaking constraint incorporates full one-loop tadpole corrections. The 
gluino, stop and sbottom masses receive full one-loop (logarithmic and finite) corrections, with approximations 
being employed in the one-loop corrections to the other sparticles. In the CP-even Higgs sector, the calculation 



is FEYNHIGGSFAST-like [|75|, 76], with additional two-loop top/stop corrections. The other Higgs' receive full 



one-loop radiative corrections, except for the charged Higgs, which is missing a self-energy correction. Cur- 
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rently, the MSSM renormalization group equations (used above Mz) are two-loop order except for the scalar 
masses and scalar trilinear couplings, which are all one-loop order equations. 



A series of points in MSSM universal supersymmetry breaking parameter space were identified [77] 
as being relevant for study, taking the results of the LEP2 collider searches (and dark matter considera- 
tions) into account. For this workshop, the parameters of each benchmark were changed until the output of 



ISASUGRA7 . 51 matched that of ref. [|77|]. The standard of these parameters is used to compare the output 
of several codes in these proceedings. We illustrate the S0FTSUSY1.3 calculation by presenting its output 
of these modified "post-LEP benchmark points" in table |l[ We use as {Mzf^^ = 0.119, mt = 175 GeV, 
mh{mi)Y'^^ = 4.2 GeV. We note that four of these points do not break the electroweak symmetry consistently. 
However, many of the points were picked specifically in order to be close to the electroweak symmetry-breaking 
boundary and so this feature is perhaps not so surprising. 

Studies of the ability of future colliders to search for and measure supersymmetric parameters have often 
focused on isolated 'bench-mark' model points [|77[-[79|] such as the post-LEP benchmarks. This approach, 
while being a start, is not ideal because one is not sure how many of the features used in the analyses will 
apply to other points of parameter space. Collider signatures typically rely upon identifying decay products 
of produced sparticles through cascade decay chains. The resulting signatures of different scenarios of SUSY 
breaking are not only highly dependent upon the scenario that is assumed, but also upon any model parame- 



ters []79|]. As a supersymmetry breaking parameter is changed, the ordering of sparticle masses can change, 
switching various sparticle decay branches on and off. In an attempt to cover more of the available parameter 
space, the Direct Investigations of SUSY Subgroup of SNOWMASS 2001 has proposed eight bench-mark model 
lines for study [||]. 

The lines were defined to have the spectrum output from the ISASUGRA program (part of the 
ISAJET7 . 51 package [^]) for mt = 175 GeV. Knowledge of the uncertainties in this calculation will be 
important when data is confronted with theory, i.e. when information upon a high-energy SUSY breaking 
sector is sought from low-energy data. Here, we intend to investigate the theoretical uncertainties in sparti- 
cle mass determination. To this end, we contrast the sparticle masses predicted by three modern up-to-date 



publicly available and supported codes: ISASUGRA7 .58*, SOFTSUSYl . 3 [JTj and SUSPECT2 .004 [|81|]. 
The asterisk indicates a changed version of ISASUGRA7 . 5 8, as detailed below. 

Each of the three packages calculates sparticle masses in a similar way, but with different approxima- 



tions [|82[]. In certain model line scenarios, we calculate the fractional difference for some sparticle s 

™ SOFTSUSYl. 3 _ ™CODE 

J-CODE _ IJj^ '"'S /IN 

•'■5 j^SOFTSUSYl.3 ' ^ ' 

where CODE refers to ISASUGRA7 . 58*, or SUSPECT2 .0 04. /™™ then gives the fractional difference of 
the mass of sparticle s between the predictions of CODE and SOFTSUSYl . 3. A positive value of /™™ then 
implies that s is heavier in SOFTSUSYl . 3 than in CODE. 

We focus upon model lines in scenarios which are currently supported by all three packages, i.e. super- 
gravity mediated supersymmetry breaking (mSUGRA). At a high unification scale Mqut = 1-9 x 10^^, the 
soft-breaking scalar masses are set to be all equal to mo, the universal scalar trilinear coupling to Aq and each 
gaugino mass Mi 2,3 is set. tan/3 is set at Mz- The three choices of model lines are displayed in Table |^ 
Model line A displays gaugino mass dominance, ameliorating the SUSY flavor problem. Model line B has 
non-universal gaugino masses and model line F corresponds to focus-point supersymmetry close to the 
electroweak symmetry breaking boundary. 



The differences in the output between three earlier versions of the codes has already been discussed [|83|]. 
Ref. [ ^ ] showed significant order 1 % numerical round-off error in the gluino and squark masses. Even worse, 
along model line F there were 10%, 3% numerical round-off errors in the lightest neutralino and chargino 
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Model 


A 


B 


C 


D 


E 


F 


G 


H 


I 


J 


K 


L 


M 
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624 


258 


415 


549 


315 


1090 
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785 
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mo 
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90 
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2970 


123 


459 


188 
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10 


10 


10 


10 


10 
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20 
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322 
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- 
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461 
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429 


847 
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531 


- 




118 


114 


119 


119 


- 


- 


119 


126 


118 


123 


- 


119 


- 




877 


379 


575 


708 


- 


- 


528 


1884 


452 


905 


- 


440 


- 
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558 
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392 
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289 




H 
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566 
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410 


013 




331 






Xi 
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99 
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654 
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319 
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n 

X2 


465 


176 


301 


411 






282 


1211 


262 


593 




347 




xl 


740 


328 


498 


636 






465 


1582 


433 


847 




530 




n 


756 


351 


516 


644 






482 


1591 
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859 




546 
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282 


1211 


262 


593 




347 
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483 
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617 
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894 
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840 


1684 
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202 


287 


376 


- 


- 


283 


1072 


300 


584 


- 


464 


- 


e-R, HR 


269 


144 


181 


238 


- 


- 


190 


703 


227 


414 


- 


391 


- 




420 


186 


277 


368 


- 


- 


272 


1069 


290 


579 


- 


458 


- 


n 


427 
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289 


376 


- 
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289 
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576 
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267 
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174 


232 


- 
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- 


240 
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458 






1252 


510 
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WR, CR 


1200 


551 


830 


1066 






791 


2767 


756 


1487 




985 
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1254 


576 


867 


1114 






825 


2905 


788 


1555 
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dR, dR 


1193 


550 


827 


1060 






787 


2748 


753 


1479 




981 




h 


1174 


583 


834 


1044 






791 


2632 


742 


1397 




903 






949 


415 


649 


856 






617 


2252 


583 


1192 




755 




hi 


1146 


523 


790 


1018 






740 


2632 


672 


1353 




884 




b2 


1190 


548 


822 


1053 






776 


2692 


722 


1400 




811 





Table 1: Post-LEP Benchmark points. Mass spectra in GeV for minimal SUGRA models calculated with program SOFTSUSYl . 3 and 
unification scale Mx — 1-9 x 10^^ GeV, Ao = 0. Columns with dashes for spectra indicate points which did not break electroweak 
symmetry correctly. All massive parameters are quoted in units of GeV. 



Table 2: Model lines in mSUGRA investigated here, rrit = 175 GeV, Mqut = 1-9 x 10^*^ GeV and a»(Mz)^^ = 0.119 are used. 



Model line 


tan (3 
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Ml 


M2 


M3 


mo 


sgn/n 


A 


10 
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Ml/2 


Ml/2 
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+ 


B 


10 





I.6M2 


M2 


M2 


M2/2 


+ 


F 


10 





Ml/2 


Ml/2 


Ml/2 


2M1/2 + 800 GeV 


+ 
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Fig. 1 : Fractional differences between tlie spectra predicted for mSUGRA model line A 
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Fig. 2: Fractional differences between the spectra predicted for model line B 



masses respectively. These numerical round-off errors were due to the ISASUGRA calculation, but this was 
not obvious because ISASUGRA was used for the normalization in the equivalent of eq. |I[ Stop masses were 
not examined. The lightest stop mass could be very important for SUSY searches, for example at the Tevatron 
collider. We now perform the comparison again, with the following differences: the output of SOFT SUSY is 
used for the normalization, up-to-date and bug-fixed versions of each code are used, we include the lightest 
stop mass in the comparison and the ISASUGRA7 .58* package is hacked to provide better accuracy in the 
renormalization group evolution^ 

We pick various sparticle masses that show a large difference in their prediction between the three cal- 
culations. For model line A, Fig. |I|a shows fl^^^^^^lp'^o ^+ g lightest stop, sbottom, squark, neutral Higgs, 
neutralino, chargino and gluino mass difference fractions respectively). Fig. |l|b shows the equivalent results 
for the output of SUSPECT. Model line B differences are shown in Fig. ^. Jagged curves in the figures are a 
result of numerical error in the SUSPECT calculation, and are at an acceptable per-mille level level for squarks, 

"We re-set two parameters in subroutine SUGRA to DELLIM=2 . Oe-3 and NSTEP=2000. 
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Fig. 3: Fractional differences between the spectra predicted for model line F 



gluinos and the lightest neutraUno. The lightest Higgs and lightest chargino do not display any appreciable 
numerical error. 

Figs, [l],^ share some common features. The largest discrepancies occur mostly for low where 
the super-particle spectrum is lightest. The gluino and squark masses are consistently around 5% lower in 
ISASUGRA than the other two codes, which agree with each other to better than 1% with the exception of 
the hghtest stop, which SUSPECT finds to be less than 4% heavier than SOFTSUSY. We note here that this 
uncertainty is not small, a 3% error on the lightest stop mass at = 700 GeV in model line A corresponds 
to an error of 35 GeV, for example. The lightest CP-even Higgs is predicted to be heaviest in SOFTSUSY, 
SUSPECT gives a value up to 6% lighter for large Mi/2> whereas ISASUGRA gives a value up to 2% lighter 
(again for large Mi/2)- This could be to some degree due to the fact that SOFTSUSY uses a FEYNHIGGSFAST 
calculation of the neutral Higgs masses with important two-loop effects added [[7^], which predicts masses that 
tend to be higher than the one-loop calculation (as used in ISASUGRA or SUSPECT). The gaugino masses 
display differences between the output of each of these two codes and SOFTSUSY, up to 4% at the lighter end 
of the model lines. 

The focus-point scenario (model line F) is displayed in Fig. ^ Fig. ^a is cut off for low /2 because 
ISASUGRA does not find a consistent solution that breaks electroweak symmetry there, contrary to the other 
two codes. The overall view of spectral differences is similar to that in model lines A and B except for the 
masses of the lightest chargino and neutralino. They display large 10-100% differences in Fig. ^. In focus 
point supersymmetry, the bilinear Higgs mass parameter // is close to zero and is very sensitive to threshold 



corrections to rrit [184|]. For small fi < Mz, the lightest chargino and neutralino masses become sensitive to 
its value. The predicted value of ^{Mz) differs by 10%- 100% between ISASUGRA and the other two codes' 
output. SUSPECT and SOFTSUSY have closer agreement, the largest differences being that the chargino is 
predicted to be 4% hghter at low M1/2 and the lightest CP even Higgs to be 4% heavier in SOFTSUSY. Only 
a few of the threshold corrections to mt are included in the ISASUGRA calculation, whereas SOFTSUSY, 
for example, includes all one-loop corrections with sparticles in the loop. SUSPECT also adds many of the 
sparticle loop corrections to mt- Because model line F has heavy scalars, another possibility for the large 
discrepancy with ISASUGRA could potentially be that ISASUGRA employs two-loop renormalization group 
equations for scalar masses, whereas the other two codes use one-loop order for them. This explanation seems 
unlikely because of the relative agreement observed in the scalar masses, which ought to be more sensitive to 
this effect. 
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To summarize, with the current technology, we do not yet have the desired accuracy for discrimination of 
supersymmetry breaking models or measurement of their parameters from the sparticle spectrum. We note that 
possible future linear colliders could determine some sparticle masses at the per-mille level [85|. An increase 
in accuracy of the theoretical predictions of sparticle masses by about a factor 10 will be necessary. 
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CoUiders 
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Abstract 

While it is natural for supersymmetric particles to be well within the mass range 
of the large hadron collider, it is possible that the sparticle masses could be very 
heavy. Signatures are examined at a very high energy hadron collider and a very high 
luminosity option for the Large Hadron Collider in such scenarios. 

1. Introduction 

If supersymmetry is connected to the hierarchy problem, it is expected [86,^] that sparticles will be sufficiently 
Ught that at least some of them will be observable at the Large Hadron Collider (LHC) or even at the Tevatron. 
However it is not possible to set a rigorous bound on the sparticle masses. As the sparticle masses rise, the fine 
tuning problem of the standard model reappears, but the sparticle masses become large enough so that they are 
difficult to observe at LHC. 



It is also possible that SUSY is the solution to the dark matter problem []88|-|90|], the stable, lightest su- 
persymmetric particle (LSP) being the particle that pervades the universe. This constraint can be applied to the 
minimal SUGRA [45,p[-94] model and used to constrain the masses of the other sparticles. Recently sets of 



parameters in the minimal SUGRA model have been proposed [ [7^ ] that satisfy existing constraints, including 
the dark matter constraint and the one from the precise measurement of the anomalous magnetic moment of 



the muon []95|], but do not impose any fine tuning requirements. This set of points is not a random sampling of 
the available parameter space but is rather intended to illustrate the possible experimental consequences. These 
points and their mass spectra are shown in Table [1} Most of the allowed parameter space corresponds to cases 
for which the sparticles have masses less than 1 TeV or so and is accessible to LHC. Indeed some of these 
points aie quite similar to ones studied in earlier LHC simulations [96, 0]. Points A, B, C, D, E, G, J and L 



fall into this category. As the masses of the sparticles are increased, the LSP contribution to dark matter rises 
and typically violates the experimental constraints. However there are certain regions of parameter space where 
the annihilation rates for the LSP can be increased and the relic density of LSP's lowered sufficiently. In these 
narrow regions, the sparticle masses can be much larger. Points F, K, H and M illustrate these regions. This 
paper considers Point K, H and M at the LHC with a luminosity upgrade to 1000 fb^^ per year (SLHC) and at 
a possible higher energy hadron collider (VLHC). We assume an energy of 40 TeV for the VLHC and use the 
identical analysis for both machines. Point F has similar phenomenology to Point K except that the squark and 
slepton masses are much larger and consequently more difficult to observe. For the purposes of this simulation, 
the detector performance at 10'^^ cm~^s~^ and at the VLHC is assumed to be the same as that of ATLAS for the 
LHC design luminosity. In particular, the additional pileup present at higher luminosity is taken into account 
only by raising some of the cuts. Isajet 7.54 [|80|,|98|] is used for the event generation. Backgrounds from ti, 
gauge boson pairs, large pT gauge boson production and QCD jets are included. 
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Table 1: Benchmark SUGRA points and masses from Ref. 
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2. Point K 

Point K has Ma ~ 2M^o and gluino and squark masses above 2 TeV. The strong production is dominated 
by valance squarks, which have the characteristic decays — > xf 9; ~^ XiQ- The signal can be 

observed in the inclusive effective mass distribution. Events are selected with hadronic jets and missing 
and the following scalar quantity is formed: 

^eff = Et + ^TJet + ^ ET^lepton 
jets leptons 

where the sum runs over all jets with Ex > 50 GeV and < 5.0 and isolated leptons with Et > 15 GeV 
andlr/l < 2.5 . The following further selection was then made: events were selected with at least two jets with 
Pt > O.lMeff, Et > O.SMeff, A(/)(jo, -Kt) < vr - 0.2, and A</)(jo, ji) < 2tt/3. These cuts help to optimize 
the signal to background ratio. The distributions in Mgj f for signal and background are shown in Figure It 
can be seen that the signal emerges from the background at large values of Me//. The LHC with 3000 fb~^ 
of integrated luminosity has a signal of 510 events on a background of 108 for M^g > 4000 GeV. These rates 
are sufficiently large so that a discovery could be made with the standard integrated luminosity of 300 fb~^. 
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Fig. 1: Mcft distribution for SLHC (left) and VLHC (right) for Point K. Solid: signal. Shaded: SM background. 



However the limited data samples available will restrict detailed studies. 
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Fig. 2: Pt distribution of hardest jet in 2jet + Et events for SLHC (left) and VLHC (right) for Point K. 



Production of qnqR followed by the decay of each squark to qx\ gives a dijet signal accompanied by 
missing Et- In order to extract this from the standard model background, hard cuts on the jets and Et are 
needed. Events were required to have two jets with pT > 700 GeV, Et > 600 GeV, and A0(ji, < O.Svr. 
The resulting distributions are shown in Figure 0. Only a few events survive at the LHC with 3000 fb^^. The 



transverse momentum of the hardest jet is sensitive to the mass [97]. The mass determination will be limited 
by the available statistics. 

The decay X2 ~^ Xi^ is dominant so we should expect to see Higgs particles in the decay of qi (qi 
X2^ ~^ Xihq)- The Higgs signal can be observed as a peak in the bb mass distributions. In order to do this, it is 
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Fig. 3: Mtb distribution for SLHC (left) and VLHC (right) for Point K. 



essential that 6— jets can be tagged with good efficiency and excellent rejection against light quark jets. There 
is a large background from ti that must be overcome using topological cuts. Events were selected to have at 
least three jets with > 600,300, 100 GeV, Et > 400 GeV, M^s > 2500 GeV, A(t>{ji,ET) < O.Ott, and 
A0(ji, ^2) < O.Gvr. The distributions are shown in Figure ^assuming the same 6-tagging performance as for 
standard luminosity, i.e., that shown in Figure 9-31 of Ref. [ ^7| ] which corresponds to an efficiency of 60% 
and a rejection factor against light quark jets of ~ 100. This 6— tagging performance may be optimistic in the 
very high luminosity environment. However our event selection is only ~ 10% efficient at SLHC and might 
be improved. There is much less standard model background at VLHC. However, there is significant SUSY 
background from g — > bib, iit which becomes more important at the higher energy. At the VLHC and possibly 
a the SLHC, it should be possible to extract information on the mass of by combining the Higgs with a jet 
and probing the decay chain ql 



X2Q 



qhxi (see e.g. [99]). 



3. Point M 

Point M has squark and gluino masses around 3.5 TeV and is beyond the reach of the standard LHC. Only 375 
SUSY events of all types are produced for 1000 fb~^ at LHC, mainly valence squarks {uL-,dL,uji,dii) and 
gauginos (xf 1 X2)- The VLHC cross section is a factor of 200 larger. About half of the SLHC SUSY events are 
from electro weak gaugino pair production mostly X2 xf ■ The dominant decays of these are X2 ~* Xi^ 
and xf Xi^^- Rates are so small that no signal close to the Standard Model backgrounds could be found 
for the SLHC. 

The effective mass distributions for Point M at SLHC and VLHC are shown in Figure ^ using the same 
cuts as for Point K. As expected, the SLHC signal is very marginal: there are only 20 signal events with 10 
background events for Mcs > 5000 GeV and SOOOfb^^. Several attempts to optimize the cuts did not give any 
improvement. Requiring a lepton, a hadronic r, or a tagged 6 jet did not help. We are forced to conclude that it 
is unlikely that a signal of any type could be observed. The VLHC signal is clearly visible and could be further 
optimized. 

The dilepton rates are shown in Fig Events are selected that have Mg// > SOOOGeV Et > 0.2Me// 
and two isolated leptons with pT > 15GeV and the mass distribution of the dilepton pair is shown. As 
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Fig. 4: Mcff distribution for SLHC (left) and VLHC (rigiit) for Point M. 
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Fig. 5: Dilepton mass distribution for SLHC (left) and VLHC (right) for Point M. Solid: l+T. Dashed: e±^T. 



expected, nothing is visible at SLHC. The distribution at VLHC is dominated by two independent decays (e.g. 
xfxf ~^ Xi^^Xi^^)' so that e+e^ + ji^ and e^/x^ rates are almost identical except for the Z peak in 
the former which arises mainly from q — > gxf ~^ Qxf^- 

On the basis of this preliminary, study it seems unlikely that Point M can be detected at 14 TeV even 
with 3000 fb"^ Higher energy would be required. 

4. Point H 

Point H is able to accommodate very heavy sparticles without overdosing the universe as the destruction rate 
for the Xi is enhanced by coannihilation with a stau. This implies a very small splitting between the fi and 



29 




Fig. 6: pt distribution of fi at SLHC (left) and VLHC (riglit) for Point H. Daslied: all fi . Solid: fi with At > 7 ns 



the Xi- In this particular case, fi -/-^ Xv^ so it must decay by second order weak processes, fi Xi^^ei^r, 
giving it a long lifetime. The dominant SUSY rates arise from the strong production of valance squarks, with 
— > x^q, X2q and qji XiQ- The staus which are produced from cascade decays of these squarks, then exit 
the detector with a signal similar to a "heavy muon". 
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Fig. 7: Calorimetric Et distributions in fi events for SLHC (left) and VLHC (right) for Point H. 

The Pt spectrum of these quasi-stable fi for 1000 fb^^ is shown in Figure ^. The ATLAS muon system 
1 97] has a time resolution of about 0.7 ns for time of flight over a cylinder of radius 10 m and half-length 20 m. 
The spectrum with a time delay At > 10(T(7ns) is also shown. Notice that this signal could be observed at 
the LHC with ~ 300 fb~^. Triggering on a slow fi may be a problem since the time- window for the trigger 
chambers is limited. However, the Et in SUSY events as measured by the calorimeter is quite large as shown 
in Figure 0. It probably is possible to trigger just on jets plus Et, the distribution for which is shown in 
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Fig. 8: Mcff distribution for SLHC (left) and VLHC (right) for Point H. Solid: signal. Shaded: SM background. 



Figure ^. The mass of the stable stau can be measured by exploiting the time of flight measurements in the 
muon measurement system. Studies of such quasi stable particles at somewhat smaller masses carried out at the 
ATLAS detector showed a mass resolution of approximately 3% given sufficient statistics (see Section 20.3.4.2 
of Ref [97]). A precision of this order should be achievable with 3000 fb^^ at either the LHC or VLHC. One 



can then build on the stable stau to reconstruct the decay chain using techniques similar to those used for the 



GMSB studies [g7|,|100p. This is not pursued here. 

The stable ri signature is somewhat exceptional so we explore other signatures that do not require it and 
would be present if the stau decayed inside the detector. For such high masses the strong production is mainly 
of u and d. Events are selected with hadronic jets and missing Et and the effective mass formed as in the case 
of Point K. To optimize this signature, events were further selected with at least two jets with pT > O.lMes, 
Et > O.SMeff, A(l){jo,ET) < n - 0.2, and A(/)(jo, ji) < 2-k/3. The M^s distributions after these cuts for 
the SLHC and the VLHC are shown in Fig Note that at the SLHC the number of events in the region where 
S/B > 1 is very small. Given the uncertainties in the modeling of the standard model backgrounds via the 
shower Monte Carlo, it is not possible to claim that the SLHC could see a signal using this global variable. The 
VLHC should have no difficulty as there are several thousand events for Me/ / > 5 TeV 

Dileptons arise from the cascade qi qx2 Q^^^'Xv The dilepton mass distributions should have a 
kinematic endpoint corresponding to this decay. Figure ^ shows the distribution for same flavor and different 
flavor lepton pairs. Events were required to have Mcs > 3000 GeV and Et > 0.2McfT and to have two isolated 
opposite sign leptons with Et > 15 GeV and \rj\ < 2.5. The structure at the VLHC is clear; the edge comes 
mainly from X2 ^ ^l^^' which has a branching ratio of 15% per flavor. This gives an endpoint at 



(M|o-M|)(M| 



Xi 



Mi 



447.3 GeV 



consistent with the observed endpoint in Figure ^ Of course this plot does not distinguish ii and in. In the 
case of the upgraded LHC, the signal may be observable, but it should be noted that the background is uncertain 
as only three generated events passed the cuts. 
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Fig. 10: Mte distribution for SLHC (left) and VLHC (right) for events containing a n for Point H. Solid: . Dashed: e='=/i'F. 



If the stable stau is used then the situation improves considerably. The dilepton mass for events contain- 
ing a fi with a time delay 7 < At < 21.5 ns is shown in Figure [l^. Since At > 10a, the standard model 
background is expected to be negligible. The SLHC signal is improved and a measurement should be possi- 
ble. The acceptance for VLHC is somewhat worse than the inclusive sample, but having the correlation of the 
dileptons with the fi should be useful. 

The VLHC gives a gain of ~ 100 in statistics over the LHC for the same luminosity at this point, which 
is at the limit of observability at the LHC. If the VLHC luminosity were substantially lower, the improvement 
provided by it would be rather marginal. The cross section increases by another factor of ~ 100 at 200 TeV. 
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5. Conclusions 



We have surveyed the signals at hadron colliders for the SUGRA models proposed by [77] concentrating on 
the cases where the sparticle masses are very large. While the masses of the sparticles at Point K are such that 
SUSY would be discovered at the baseline LHC, the event rates are small and detailed SUSY studies will not 
be possible. The reach of the LHC would be improved by higher luminosity where the extraction of specific 
final states will become possible. The cross section at a 40 TeV VLHC is approximately 100 times larger than 
that at LHC. This leads to a substantial gain, but it is important to emphasize that this gain requires luminosity 
at least as large as that ultimately reached by the LHC and detectors capable of exploiting it. Point H has a 
special feature in that the stau is quasi-stable. This feature would enable a signal to be extracted at SLHC. If the 
tau mass were raised slightly so that its lifetime were short, then only the VLHC could observe it. The masses 
in the case of Point M are so large that the VLHC would be required for discovery. Point F has a gluino mass of 
order 2 TeV and should be observable at the LHC exploiting the production of gluinos followed by the decays 
to Xi and hence to leptons. 

The Points A, B, C, D, G, I, and L which are much less fine tuned have similar phenomenology to the 
"Point 5" or "Point 6" analysis of [97] in that lepton structure from the decay X2 ~^ ^R^ ~^ ^^^ X? and/or 



X2 — > ^ T^r Xi is present. In most cases decay X2 ~^ is also allowed, so that a more complicated 
dilepton mass spectrum is observable. This should enable the extraction of m^^ in addition (for an example 



see Fig 20-53 of [|97|]). Points A, D and L have higher squark/gluino masses and will require more integrated 
luminosity. Nevertheless one can have confidence that the baseline LHC will make many measurements in all 
of these cases. 
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Part VI 

SUSY with Heavy Scalars at LHC 

/. Hinchlijfe and RE. Paige 

Abstract 

Signatures at the LHC are examined for a SUSY model in which all the squarks and 
sleptons are heavy. 

1. Introduction 

SUSY models may give new contributions to flavor changing neutral currents, CP violation, etc., through 
loops involving squarks and sleptons. These effects are reduced if the scalars are heavy. The "inverted hi- 



erarchy" [jlOlj] and "focus point" [ ]102| ] scenarios provide examples of ways in which heavy scalars could be 
accommodated naturally. 

This note examines the LHC signatures for a minimal SUGRA model with 

mo = 1500 GeV, mi/2 = 300 GeV, = 0, tan/3 = 10, sgn^ = +. 
The gaugino masses are similar to those considered previously, e.g., 

M(x?) ~ 109 GeV, M{xt) ~ 161 GeV, M(x^) « 289 GeV, M{g) ^ 782 GeV. 
Most of the scalars have masses around 1500 GeV; the light Higgs mass is 116 GeV. 

The SUSY production cross section at the LHC is dominated by gluino pairs. The two largest branching 
ratios are 

B{g Xitb + h.c.) B{g x^t^ + h-c) ^ 23%. 

However, decays into both charginos and all four neutralinos with all allowed combinations of quarks are 
significant. This leads to many complex signatures. 

IS AJET 7.51 was used to generate events for the signal and for all the Standard Model (SM) backgrounds. 
The detector response to these events was simulated using a parameterized simulation with parameters appropri- 
ate to the ATLAS detector. Jets were found using a simple cone algorithm with R = 0.4. Lepton identification 
efficiency and b and r jet tagging and misidentification were included with parameterized efficiencies and back- 
grounds based on full simulation of ATLAS. A micro-DST was saved and subsequently analyzed using Root as 
a framework. The statistics for the signal correspond to approximately 100 fb^^. The statistics for the largest 
SM background samples correspond to a much smaller luminosity but are sufficient to show that the Standard 
Model backgrounds are small after cuts. 



2. Effective Mass Distribution 

An inclusive signature based on multiple jets plus missing energy Et was useful at many of the SUSY points 
considered previously and remains so here. Since the jet multiplicity is higher here, the effective mass was 
defined to include all jets and leptons, not just the four hardest jets: 

Afjot Mop 

1=1 i=l 

Events were selected to have 
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Fig. 1: Left: Mgff distribution for signal (curve) and SM background (shaded). Right: Same with > 1 6 tags. 

• At least six jets with pT > 100, 50, 30, 30, 30, 30 GeV; 

• Et> max(100GeV,0.2Meff); 

• Transverse sphericity St > 0.2; 

• Meff > 1000 GeV. 

The resulting M^s distribution for signal and background is shown in Figure [I[ In contrast to many previous 
cases, the signal emerges from the SM background well past its peak, but nevertheless the S/ B ratio is large 
for large enough A/eff- Thus, discovery of a deviation from the SM is easy, although not quite so easy as in 
earlier cases. 

This signal can be improved by requiring at least one b jet. A b tagging probability of 70% was chosen, 
and the corresponding light jet rejection was taken from full simulation results for ATLAS. This distribution is 
also shown in Figure [l[ As expected, the S/B ratio is improved with only a small loss of signal. 

3. Top Reconstruction 

Given the large branching ratios for ^ —>■ x^t^ + he, it is natural to tiy to reconstruct hadronic top decays. If 
everything decays hadronically, the jet multiplicity from each gluino is 6 for xf and 8 for , giving a total 
of 12 to 16 jets without any gluon radiation. This produces a severe combinatorial background; lepton-based 
signatures are considerably easier. Work on top reconstruction is continuing. A more sophisticated jet algorithm 
might work better for these complex events. 



4. e+e- Signature 

The three-body decay X2 ~^ Xi^~^^~ has a kinematic endpoint at M(x2) — -^(x?) = 61.5 GeV. Events 
satisfying the cuts given in Section ^ were required to have two OS,SF leptons with pT > 15 GeV and \r]\ < 
2.5. The reconstruction efficiency was assumed to be 90% for both e and /x. Figure § shows the resulting 
and e^^^ mass distributions. Any contribution from two independent decays should cancel in the difference 
of these. This difference shows both a continuum with an endpoint at the expected place and a Z peak coming 
from decays of heavy gauginos. 
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Fig. 2: Left: Dilepton mass distributions for I'^E from signal (solid), e^fi^ from signal (dashed), and i'^i from SM background 
(shaded) after cuts. Right: £^£~ — £^£^ distribution (solid) and with > 1 6 tag (dashed). 

The largest single source of X2^^ 9 ^ X2*^; while the largest sources of heavy gauginos are g ^ xt^^ 
and g —>■ x'^ti. Thus one expects a large fraction of dileptons to be accompanied by a b. Figure |^ also shows 
the subtracted distribution without and with at least one b tag. 

5. e^fi^ — e^fj,^ Signature 

Two independent leptonic decays of the same gluino, e.g., g — > Xj^tb with xi —>^e~X and t — > fi~^X, gives 
an OS dilepton signature. Since the ^ is a Majorana fermion, any contribution from leptonic decays involving 
both gluinos will cancel in the combination e^fi^ — e^/i^. (Equal acceptance for e and fi is assumed here. In 
reality one would have to correct for acceptance; this correction can be checked using Z ^ e+e^ , fi^n^ data.) 
The resulting distribution using the cuts described above is shown in Figure |3[ 

While the e^/x^ — e^fi^ dilepton distribution should have a true kinematic endpoint corresponding to the 
maximum possible mass from gluino decay, this is not useful because many particles are unobservable (u, x?) 
or not included (jets). The largest contributions to this channel should come from g — > xf ^^^^ 2 ~^ 
and t ii/b. Three samples of gxi events with 200 < pr < 600 GeV (the typical range for the gluino pt) 
and with one of the three decay chains forced were generated. The same analysis was applied except that the 
required number of jets was reduced from 6 to 3. The mass distributions for the three possibilities are shown 
in Figure ^ All three are qualitatively similar to that shown in Figure ^ The shapes are somewhat different 
and presumably could be distinguished with sufficient statistics after a detailed analysis. This has not yet been 
attempted. 

The sign-subtracted e/i pair was next combined with each of the three hardest jets (with px > 100 GeV) 
in the event. The distribution of the minimum of the three masses is shown in Figure |5[ The distribution in 
the case that the jet giving the minimum is tagged as a 6 is also shown in the Figure. If one of the jets is 
from the same gluino as the dilepton pair, then this distribution should have a kinematic endpoint related to the 
gluino mass. The choice of three jets is a compromise between including the right jet and including too many. 
The expected shape from a single gluino was again determined using the g + Xi sample; this is also shown in 
Figure A similar analysis combining the e/i with two jets found too much combinatorial background. 
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Fig. 3: e^fj,^ — e^fi^ distribution for signal (solid) and SM background (siiaded). 
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Fig. 4: Left: e^fi^ — e^/i* mass distribution for three possible signal contributions. Right: Sum weighted by branching ratios. 



If one of the leptons is from t Wb, then the smallest £j mass should be less than the kinematic limit 
for this decay, ^ (ml — m^)/2 = 110 GeV. This minimum mass is plotted in Figure |6| and has the expected 
shape. However, a rather small fraction of the jets so selected are tagged as 6's, while the b tagging efficiency 
is about 60%. 



6. i^i^ + jets Signature 

If both gluinos decay via g — > xf with xf — ^ Xi^^^' the signature is four hard jets plus two leptons. 
Requiring the leptons to be the same sign causes the loss of half of the signal but greatly reduces the SM 
background. Events were selected to have 
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Fig. 5: Left: Minimum efij mass formed from the sign-subtracted e/i pair, Figure ^ and one of the three hardest jets with pr 
100 GeV. Solid: all jets. Dashed: with b tag. Shaded: SM background. Right: Same for g + Xi events with forced decays. 
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Fig. 6: Minimum £j mass from sample in Figure ^. Solid: all jets. Dashed: with b tag. 



• Four jets with pr > 40 GeV, the first with pr > 100 GeV; 

• Mcff > 500 GeV; 

• > max(100GeV,0.1Meff); 

• St> 0.2; 

• At least 2 leptons; 

• Less than 2 tagged b jets. 

The two hardest leptons were required to be the same sign. For each of the three possible ways of pairing 
the jets, the larger of the dijet masses was taken, so the minimum of the three masses should be less than 
the dijet endpoint for gluino decay, M{g) — M(xf ) = 620 GeV. The distributions for all three and for the 
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Fig. 7: Dijet mass distributions for like-sign dilepton events for signal (solid) and SM background (shaded). Left: all three combinations. 
Right: minimum combination. 




Fig. 8: Same as Figure |7| for forced decays. 

minimum are shown in Figure 0. The expected distribution obtained by forcing the decays g —>■ xt'^d and 
xf — > Xie~^i^ is shown in Figure |8[ The endpoint has the expected value. However, the sample is not very 
pure: Figure ^ shows that there are other contributions. Harder cuts on extra jets did not help significantly to 
improve the purity. 



7. j Signature from g 

The decays g —>■ X29^ X39 have branching ratios of 



has a branching ratio of ~ 100%. The £'^1 
as a b. Figure 10 shows the resulting (e+e 



L% and 2% respectively at this point; X3 —>■ Xi^ 
pair was combined with any jet with pT > 200 GeV not tagged 
+ — e^fi^) + j mass distributions for the Z peak and 
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Fig. 9: SUSY particle content for events selected in Figure ^ using the ISAJET particle numbering scheme. 
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Fig. 10: Minimum mass of e^e~ +/^'*'/^~ — e^/x^ combined with jet with pT > 200 GeV. Left: Z peak. Right: Below X2 ^ Xi^^^ 
endpoint. 



for the X2 ~^ Xi^^^ continuum. The Z distribution should have an endpoint at 602.3 GeV that can be 
calculated in terms of the masses involved. The continuum is more comphcated since the dilepton mass also 
has a distribution. 

8. xt + X2 Signature 

Direct production of gauginos is only a factor of ~ 10 smaller than gg production at this point. Events were 
selected to have three leptons with pT > 50, 20, 20 GeV. A jet veto of 30 GeV was imposed. The OS,SF 
dilepton mass distribution is shown in Figure 11. The dilepton endpoint is known from Figure^ Requiring an 



OS,SF pair below this endpoint gives the £^ distribution in the same figure. Clearly the SM background 



40 




M|| (GeV) M||| (GeV) 

Fig. 1 1 : Mass distributions for £'^£~' (left) and i'^l'i^ (rigiit) witii jet veto for signal (curves) and SM background (shaded). 

is still comparable to the signal; it would be worse if the Z width in SM WZ events were properly taken into 
account. Also, the effect of pileup on the jet veto has not been included. Thus this channel does not seem very 
promising. 
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Part VII 

Inclusive study of MSSM in CMS 

S. Abdullin, A. Albert, F. Charles 

Abstract 

The Minimal Supersymmetric Standard Model is an extension of the Standard Model, 
the most economical one in terms of new particles and new couplings. Many studies 
have been performed on the observation of supersymmetry, but mostly limited to the 
mSUGRA model. Here we consider the possibility of a broader test of SUSY, using 
a less constrained model than mSUGRA, the pMSSM (phenomenological MSSM). 
This study is made in an inclusive way in the framework of the CMS experiment. 
We first show the ability of CMS to discover SUSY in a large domain of pMSSM 
parameter values. We then attempt to estimate the uncertainties in the determination 
of MSSM parameter values using essentially kinematical measurements. 



1. Aim 

The MSSM is a good candidate for the new Physics expected at the TeV scale. Experiments at both LEP 



| |103| ] and Tevatron | |104[ ] have been looking for evidence of SUSY, but for the moment no signal has been 



observed. The results of the searches give limits on sparticle masses. Some indirect measurements like the 



6 — > S7 branching ratio, the anomalous muon magnetic moment g-2 [ 177| , |105| ], dark matter searches [ |106| ] or 
the Z width also provide constraints on SUSY parameter values. But all these results still leave a large MSSM 
parameter space unexplored. 



The goal of this study is to evaluate the ability of the CMS detector [107] to observe signals of super- 
symmetry in a large domain of MSSM parameter values. The mSUGRA [181, |108| | (minimal SUper GRAvity) 
model, with its only five free parameters (mo, m'i/2. ^o. tan/3 and sign(/i)) is very popular and has been 



the subject of many studies up to now. A study of mSUGRA, performed in a similar way [ |109| ] to the work 
presented here, concludes that for both low and high values of tan (3, and for both positive and negative /^i, the 
mass reach for gluinos and squarks is up to « 2.5 to 2.7 TeV for lOOfb"^. 

Two reasons motivate us to extend the mSUGRA study to a less constrained model, the "phenomeno- 
logical" MSSM (pMSSM). On the one hand, mSUGRA is a rather constrained model, very specific and not 
illustrative of the variety of all possible supersymmetric models. On the other hand, contrary to mSUGRA, 
the pMSSM has no fixed hierarchy of masses. In this case, some extreme mass hierarchies could show a sig- 
nificantly different kinematical behavior than in the case of mSUGRA, which could prevent the discovery of 
supersymmetry even for relatively low values of the spaiticle masses. Moreover, in the case of pMSSM we have 
various types of cascades which produce many types of final states, with similar signatures, but not the same 
types of particles. For example, are we going to be able to observe supersymmetry with a final state containing 
multiple jets; taus instead of muons and electrons; and jets produced by c quarks instead of b quarks? The type 
of particles produced is really important in order to discover supersymmetry in CMS and to identify the SUSY 
scenario at work. 

In this pMSSM framework, we are going to show that supersymmetry could be discovered over a large 
scale of masses in the versus nig plane. Next we show that there are ways to estimate the values of the 
MSSM parameters using kinematical quantities measured by the CMS detector and event rates. An advantage 
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of this approach is its model independence - the only dependence comes from the hierarchy of masses [IIC] 



Finally, we estimate the statistical uncertainties due to this method of extraction of the MSSM parameter values. 

2. Theoretical framework 
2.1 MSSM 

Supersymmetry (SUSY) is a symmetry between fermions and bosons. Some of the motivation for SUSY has 



been reviewed in the Introduction. The MSSM D81| , |108| ] is the Supersymmetric extension of the Standard 



Model which introduces the minimal number of new particles (only one per SM particle and 4 additional Higgs 
bosons) and no new couplings. The MSSM contains 124 independent parameters, including the 19 ones of the 
Standard Model. 

2.11 pMSSM (phenomenological MSSM) 

Some phenomenological constraints allow to reduce the number of MSSM free parameters: 

• no new sources of CP violation, 

• no Flavor Changing Neutral Current effects, 

• universality of the first two generations. 

These three constraints leave only 19 free parameters : 

• tan (3 : ratio of the vacuum expectation values of the two Higgs doublets fields, 

• Ma '■ mass of the pseudoscalar Higgs boson, 

• ^ : SUSY preserving Higgs mass parameter, 

• Ml, M2, M3 : bino, wino and gluino mass parameters, 

• Mq Mujf,M^-^, M^, Mejf : unified first and second generation sfermion masses, 

• Mf^, M^-^, Mq, Mj^, M-f^ : third generation sfermion masses, 

• At, Afj, At-: third generation trilinear couplings. 

2.12 A restricted pMSSM 

The model used in our study is a pMSSM, but with a further reduction in the number of free parameters. 
It is an intermediate model between mSUGRA and the pMSSM, relaxing the constraints of mSUGRA but 
still more constrained than the pMSSM. Reference gives some examples of such models, taking into 
account more constraints than pMSSM. We take into account, respectively, the mass unification of squarks and 
sleptons (universality of the three generations of sparticles) assuming that the mixing is not too large for the 
third generation. We also consider the unification of trilinear coupling At = A^ = A-^-. This leads to 9 free 
parameters: tan/3. Ma, ^, Mi, M2, Mg, Mq, M^, A3. This constrained model allows us to perform simpler 
simulation, while keeping the diversity of signatures of MSSM events. 

2.2 Examples of signal events 

2.21 An example of MSSM cascade 

Figure [l] shows an MSSM event of the type gq —>■ q^g, with 5 jets including 2 b quark jets and 3 leptons 
including 2 r's in the final state. The 2 neutralinos x? produce missing transverse energy E^"^^^. 
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Fig. 1: An example of an MSSM cascade resulting from gqL production. 




(a) Another example of an mSUGRA cascade. 



(b) GEANT output for the mSUGRA event shown in 
Fig.ra 



Fig. 2: Another detailed example of an mSUGRA event. 



2.22 A more detailed mSUGRA example 

An mSUGRA event of the type gq — > gq is shown in figure 2(a) , while figure 2(b) shows the corresponding 
event display in CMS obtained from GEANT 1 1 11 1 for this event. We used the following parameter values: 

mo = 1000 GeV, mi/2 = 500 GeV, Aq = 0, tan /3 = 35, /i > 0. 

The final state is made of 6 jets including 4 b-quark jets, and 2 neutralinos which produce E"^^^^ . 
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3. Simulation procedure 
3.1 Signal production 

We use a model with 9 parameters, which make up a hyperspace in 9 dimensions. To simplify the analysis we 
use a discretization of the parameter values. The choice of the number of values for each parameter depends 
on the parameter sensitivity. We used a grid for squark and gluino masses with 9 values evenly spaced between 
600 and 3000 GeV, because the event characteristics at LHC depend primarily on these two masses. On the 
other hand, many observables are not very strongly dependent on the parameter tan f3. We thus use only two 
values, to distinguish the behavior at large and small values of this parameter. The values selected for each 
parameter in this analysis aie the following: 

• Mq : 600, 900, 1200, 1500, 1800, 2100, 2400, 2700, 3000 GeV 

• Mg : 600, 900, 1200, 1500, 1800, 2100, 2400, 2700, 3000 GeV 

• : 200, 1000, 3000 GeV 

• Ml : 100, 500, 1000, 2000 GeV 

• Ms : 100, 500, 1000, 2000 GeV 

• Ma : 200, 1000, 3000 GeV 

• As : 0, 2000 GeV 

• /i : 200, 500, 2000 GeV 

• tan (3: 2, 50 

We end up with a total of 140000 different sets of parameter values. For each set, we generate 1000 
events, a compromise between the limits imposed by the handling of the data flow and sufficiently small statis- 
tical errors. The theoretical and experimental constraints make it possible to reduce the number of combinations 
toatotalof 17.10^. 

These imposed constraints are the following: 

• constraint on the Higgs mass > 100 GeV, 

• lightest chargino mass > 100 GeV, 

• lightest neutralino mass > 50 GeV, 

• the lightest neutralino is the LSP. 

Signal events were generated using the ISAJET program [|98|]. The CMS detector response was obtained 



from the fast Monte-Carlo code (non-Geant) CMS JET 4.51 0TT2| ]. Chai-acteristics of CMS JET software are 
given in figure ^ 

3.2 Background production 

The background production to this g, q SUSY search was estimated using Standard Model events leading to 



similar signatures as the MSSM events. The background was produced using PYTHIA [ |1 13| ]. We consider 
here SUSY signals with the following event characteristics: 

• production from to n isolated leptons (electrons or muons), 

• a large value for the average missing transverse energy (> 200 GeV), 

• more than 2 jets with large transverse energy (> 40 GeV). 

Therefore we must consider as potential Standard Model backgrounds all processes with large rest 
masses which yield large transverse missing energy, producing energetic jets, and possible isolated leptons. 
Thus the backgrounds we consider are the following: pp — > tt,W + jet, Z + jet, WW, ZZ, ZW. We also 
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Fig. 3: Details of CMSJET characteristics. 



consider QCD events with several high energy jets, including heavy flavors (b and c). The missing transverse 
energy in this type of event originates either from semi-leptonic b,c decay or from imperfections and fluctu- 
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Fig. 4: Feynman diagrams of a few baclcground examples. 



ations in the response of the detector which may fake missing transverse energy. Some of the background 
mechanisms are shown in figure g The primary contributions to the background after cuts is at large values of 
the transverse momentum px of the produced particles in the 2-body final state. To obtain sufficient statistics, 
we generate events independently for several intervals of pT (see table |l]). A total of two hundred million events 
have been generated. 

The QCD event sample generated for low px is tiny compared to the required one but fortunately, there 
is correlation between pT and the maximal produced E'^*^* value, so one does not expect high values of E'™*''** 
for low pt- To be confident in the simulation, we apply some preliminary cuts during the generation: 

• EJp^'^^ > 200 GeV limit, below which the QCD jet background become dominant; 

• at least 2 jets with E^"^* > 40 GeV in l??-?^*] < 3. 

The isolation of the leptons is given by the following requirement 

• muon with pj, > 10 GeV within the muon acceptance or electron with p^ > 20 GeV within \r]'^\ < 2.4; 

• no charged particles with pT > 2 GeV in a cone of R = 0.3 around the direction of the lepton; 

• ^ E!^'-'- in a cone ring 0.05 < R < 0.3 around the lepton impact point has to be less than 10% of the 
lepton transverse energy. 



3.21 Pile-up 

We also take into account event pile-up, i.e. 25 inelastic pp events on average per bunch crossing with a Poisson 
distribution. The two upper graphs in figure || illustrate the ratio between the lepton isolation efficiency with 
and without pileup, as a function of pseudorapidity, and transverse momentum (the definition of lepton isolation 
is given in [ |1 12 |). The efficiency is reduced to 85% due to the multiplicity of particles produced in each bunch 
crossing. In the two lower plots giving the event missing transverse energy, and the scalar sum of the event 
transverse energy, the solid curves are without pile-up and the dashed ones with pile-up. Pile-up does not make 
a very significant difference for total missing transverse energy, but increases the total transverse energy. 



3.3 Different selection criteria 

To distinguish signal from background, we are led to apply kinematical cuts on the observables we extract from 
the CMS detector. 
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processes 


interval (ueV) 


a (pb) 
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% oi neeaea lor lUO ib 
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100 - 200 


1.37-10'* 


6.000 • 10'^ 
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7.15 • 10^ 


3.229 • lO'^ 


0.45 


QCD 


400 - 800 


2740 


3.259 • 10"^ 


11.9 




800 - 1200 


60.0 


6.033 • 10** 


100.0 




> 1200 


4.8 


4.947 • 10^ 


103.1 


total 






2.342 • 10*^ 





Table 1: Background repartition by pr interval. For each pr interval we give the process cross-section in pb and the number of 
generated events. The last column shows the percentage of events we have generated compared to the expected number of events for 
100 fb"'^ of integrated luminosity. 



Table ^ gives all the different values we use for the selection. This yields approximately 10000 combi- 
nations of cuts to optimize signal to background ratio. 

3.4 Signal significance estimator 

We make a systematic search for all sets of parameters and thus define limits of discovery through calculations 
of the significance of the signals. The definition of the significance we use is the following one: 

significance = —=, (1) 

V B 

where S is the number of signal events and B the number of background events. A significance exceeding 
5 indicates that the corresponding set of MSSM parameter values is experimentally accessible. In order to 
optimize the significance, we used some cuts which are listed in the Table ^ 

To show the importance of cuts to separate the signal from the background, tables |3| and ^ give examples 
of the number of events for signal and background, before and after cuts. We notice the very important effect 
of the cuts, the number of background events decreasing from 2.4 10® to ^ 100, i.e. by a factor roughly equal 
to 2 X 10®, and the significance increasing by a factor « 100. 
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Fig. 5: Importance of pile-up for lepton isolation and energy measurements. 



parameter 


different cut values 


Number of jets 
(with a minimal transverse energy of 40 GeV) 


2, 3,4,5,6, 7, 8, 9, 10 jets 


Transverse momentum of the highest energy jet 


40, 150, 300, 400, 500, 600, 700, 800, 900, 1000 GeV 


Transverse momentum of the second highest energy jet 


40, 80, 200, 300, 400, 500 GeV 


Missing transverse momentum 


200, 300, 400, 500, 600, 700, 800, 

900, 1000, 1100, 1200, 1300, 1400 GeV 


Total transverse momentum 


700, 900, 1100, 1300, 1500, 1700, 1900, 2100, 
2300, 2500, 2700, 2900, 3100, 3300, 3500 GeV 


Angle $ between the missing transverse 
momentum and the momentum of the isolated lepton 


0, 20 degrees 



Table 2: All sets of cut values for significance optimization. 



4. Analysis 

4.1 Calculation of significance 
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before cuts 


after cuts 


number of signal events 


6152 


431 


number of background events 


240 10^ 


124 


significance 


0.397 


38.614 



Table 3: An example of the effect of cuts on the number of signal and background events. We give also the significance obtained. The 
MSSM parameter values are the following: A/,- = 1000 GeV, Ah = 500 GeV, M2 = 1000 GeV, Mg = 2700 GeV, Mg = 900 GeV, 
Ma = 200 GeV, tan /3 = 50, p = 500 GeV, ^3 = GeV. 





before cuts 


after cuts 


number of signal events 


6121 


355 


number of background events 


240 10^ 


112 


significance 


0.395 


33,5 



Table 4: An example of the effect of cuts on the number of signal and background events. We give also the significance obtained. The 
MSSM parameter values are the following: Mj = 1000 GeV, Mi = 100 GeV, M2 = 500 GeV, Mg = 2100 GeV, Mg = 2100 GeV, 
Ma = 1000 GeV, tan/3 = 50, p = 2000 GeV, A3 = 2000 GeV, which correspond to Fig.|6(a)[ 



4.11 Illustration of the analysis for some specific parameter values 



Figure 5(a) gives distributions of signal (in black) and backgrounds (in gray) for some kinematical quantities 
before any cuts are applied; the signal is not easily distinguishable from the background at this stage, as the 
cross section is too much smaller. But quantities such as E'^^'^^ and E'|"™ have a very different shape, thus 
cutting on these variables would greatly enhance the signal to background ratio. 

The specific example shown in figure 5(a)| corresponds to the following values of MSSM parameters: 



Ml = 1000 GeV, Ml = 100 GeV, M2 = 500 GeV, Mq = 2100 GeV, Mg = 2100 GeV, 

Ma = 200 GeV, /x = 2000 GeV, A3 = 2000 GeV, tan /? = 50. (2) 

The cross section for this set of parameters is o" = 58 fb, and the significance after applying all cuts is equal to 
33. 

4.12 Example with either very broad or very narrow hierarchy of masses 

We have investigated in some detail one of the major points of difference between mSUGRA and pMSSM, 



namely the non fixed hierarchy of masses in case of pMSSM. In the first example (figure p(b)| ), with a very 
broad mass spectrum, the masses of neutralinos are chosen to be much lower than the masses of squarks, 
gluinos and sleptons. Sparticles production is therefore dominated by neutralinos and charginos. The specific 
parameter values are the following: 

Ml = 2000 GeV, Mi = 500 GeV, M2 = 500 GeV, Mg = 2000 GeV, Mq = 2000 GeV, 

Ma = 1000 GeV, = 200 GeV, A3 = GeV, tan (3 = 50. (3) 

The cross section for this set of parameters is it = 1.22 pb and despite the abundance of neutralinos and the 
low production rate of gluinos and squarks, one is still able to obtain after appropriate cuts (discussed in the 
following section) a significance equal to 10.2. 
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(a) Sample without specific hierarchy. AIj — 1000 GeV, 
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(b) Sample with broad hierarchy of masses. Mj — 2000 
GeV, Ml = 500 GeV M2 = 500 GeV, Mg = 2000 
GeV, Mij = 2000 GeV Ma = 1000 GeV, tan /? = 50, 
= 200 GeV, ^3 = GeV. 
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(c) Sample with narrow hierarchy of masses, low value. 
1500 GeV, Ml = 940 GeV, M2 = 2000 GeV, 
1020 GeV, Ma = 1000 GeV, 
tan 13 = 50, fi = 1050 GeV, A3 = GeV. 



Mi 



(d) Sample with narrow hierarchy of masses, medium value. 
1520 GeV, Ml = 1450 GeV, M2 = 2000 GeV, 



Mi 
Mi 



1500 GeV, Ms 



1520 GeV, Ma = 1000 GeV, 



tan P = 50,fi = 1500 GeV, A3 = GeV. 



Fig. 6: Distribution of signal and background for different observables before cuts. 
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As a second example (figure |6(c)D we chose a case where the masses of the neutralinos, gluinos, squarks 
and sleptons are comparable. The specific parameter values are the following: 

Ml = 1500 GeV, Mi = 940 GeV, Ms = 2000 GeV, Mg = 1000 GeV, Mg = 1020 GeV, 
Ma = 1000 GeV, /i = 1050 GeV, A3 = GeV, tan (3 = 50. (4) 

SUSY production now mostly proceeds via gluinos and squarks with a cross section a = 2.014 pb and a 
significance equal to 36.3 after selection cuts. 



For the third example (figure 5(d)) the parameter values are the following: 

Ml = 1520 GeV, Mi = 1450 GeV, M2 = 2000 GeV, Mg = 1500 GeV, Mg = 1520 GeV, 
Ma = 1000 GeV, = 1500 GeV, A3 = GeV, tan (3 = 50. (5) 

The masses of the neutralinos, gluinos, squarks and sleptons are similar again but heavier than in the previous 
example. The main sparticle production proceeds still via gluinos and squarks, with, in this case, a cross section 
a = 0.126 pb and a significance equal to 3.2 after cuts are applied. 

Even for the sets of parameter values which would seem difficult (either very similar masses or on the 
contrary very broad span of masses), applying cuts allows to obtain good results. However, in the case of a 
small spread of sparticle masses, the discovery limit is about 1.5 TeV instead of 2.5 TeV as obtained with the 
usual mSUGRA-type mass hierarchies. 



4.13 Discovery limits 

We now generalize this study to determine a Umit of discovery for the MSSM. For each MSSM point, we are 
looking for the set of cuts which gives the maximum value for the significance. With these collected values of 
the maximal significance, we can draw the discovery limits (estimated as the isocurve of a significance equal 
to 5). The significances are calculated for an integrated luminosity J Ldt = 100 fb^^ corresponding to one 
year of LHC at high luminosity. Our first result is the isocurve of significance equal to 5, given in the plane 
(Mg, Mg), which are the two most important parameters. The 7 other parameters have the same fixed value for 
all the MSSM points of the plane. One example of such a result is given by Figure ^ with the following values 
of the other parameters: 

Ml = 3000 GeV, Mi = 500 GeV, M2 = 1000 GeV, Ma = 200 GeV, 

^l = 500 GeV, A3 = GeV, tan [3 = 50. (6) 

Four isocurves are given in this figure, corresponding to a specific event topology selection according to 
the number of isolated leptons produced. The curves labelled "0 lepton", "1 lepton", "2 leptons" correspond to 
a calculation of significance using only events with respectively 0, 1, 2 lepton(s). The curves labelled "all" use 
all events to calculate the significance. For each point and each type of lepton selection, we manage to find the 
set of cuts which gives the largest significance. 

We can now try to compile all these results in one characteristic limit. We could first combine our results 
in a conservative fashion, establishing the region in the (Mg,Mg) plane in which any set of the orthogonal 
parameters will be accessible by the CMS detector. In other words, for each point in the (Mg, Mg) plane above 
the reach curve there exist at least one set of the 7 other parameter which has a significance under 5. Figure |8] 
shows this pessimistic mass reach in the (Mg, Mg) plane. In terms of statistics, there are a total of 2,962 
combinations (out of the original 35,000 combinations of parameters) which don't pass the cut on significance 
(> 5). We have also found a single point under the curve which does not pass the significance cut: 

Ml = 3000 GeV, Mi = 2000 GeV, M2 = 2000 GeV, M-g = 1500 GeV, Mg- = 1800 GeV, 
Ma = 3000 GeV, = 200 GeV, A3 = 2000 GeV, tan (3 = 50. (7) 
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Fig. 7: Graph of discovery limits in the Mq versus Mg plane. 

This particular point exhibits a narrow hierarchy of masses just like in the previous example, and provides a 
lower limit of discovery. In a typical configuration we have a discovery limit of about 2.5 TeV. This exact reach 
in parameter space depends on the magnitude of the background cross section within the kinematical cuts. Here, 
we assume that the PYTHIA cross section are correct. This is clearly invalid, as higher order QCD corrections 
to t i, W+jets, Z+jets are not incorporated. This is an aspect of systematic uncertainties to be addressed in a 
later study. 

4.2 From kinematical observables to MSSM parameters 

4.21 Choice of observables 

In a second exercise, we reverse the problem and try to see whether on the basis of event kinematical variables 
and event rates it would be possible to determine the MSSM parameter values. A total of 11 observables are 
used to separate the different sets of MSSM parameter values. 

• average number of leptons per event {Ni), 

• average number of jets per event {Nj), 

• mean value of jet momenta (Ptj), 

• mean value of lepton momenta {Pa), 

• mean value of missing transverse energy 

• number of events Ntot, 

• number of events with 0,1,2 or 3 leptons No, Ni, N2, N^, 

• mean value of total transverse energy 



- lepton (solid) 



- 1 lepton (dashed) 



- 2 leptons (dot-dashed) 




Mi=3000., M,=500., M2=1000., Mnj=200. 
tanp = 50., n = 500., A3 = 0. 
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Fig. 8: Discovery limits of the MSSM in the (Mg, Mg) plane. 



These observables are characteristic of the measurement to be done with the CMS detector, and corre- 
spond to signatures of MSSM events. In particular, there are correlations between these observables and the 
MSSM parameters. 

Figure ^ illustrates some of these correlations between observables and pMSSM parameters, more specif- 
ically, between {Ef^^), Ni, {Ni), {Nj) and the parameter fi. The values of the other parameters are: 

M|- = 3000 GeV, Mi = 500 GeV, Ms = 500 GeV, Mg = 2100 GeV, Mg = 2400 GeV, 

Ma = 200 GeV, A3 = 2000 GeV, tan p = 50. (8) 

It can be seen that the value of all these kinematical quantities decreases with increasing fi. This shows that 
these observables are sensitive to the value of fi. We have obtained the same behavior with the other pMSSM 
parameters, which is an argument to use pMSSM instead of mSUGRA. 

4.22 Separation of parameters 

After optimizing the cuts to achieve maximal significance, we extract the values of all observables for signal and 
background, and we calculate the statistical uncertainties cje associated with each observable; the uncertainties 
C7e are equal to the ratio of the standard deviation of the distribution and the number of events, thus : 

• o"e oc ^fN for N number of events (for example the number of events with leptons); 

• (je oc in general and for mean value ( like (E'j"**'^)) with ai the root mean square. 

In this way we take into account the uncertainties on the averages derived from the small number of 1000 events 
generated for each MSSM point. Can these values of the observables be linked to the values of the pMSSM 
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Fig. 9: Relationship between some observables and the fj, parameter for Mj = 3000 GeV, Mi = 500 GeV, M2 = 500 GeV, Mq = 2100 
GeV, Mg = 2400 GeV, A/a = 200 GeV, tan/3 = 50, A3 = 2000 GeV. 



parameters? In other words, are we able to distinguish two different sets of MSSM parameter values using only 
the values of the observables for each set? If this discrimination is possible, we associate each set of parameter 
values with the corresponding values of kinematical quantities, and also we use interpolation when we turn 
to the continuous case. Discrimination is carried out in the following way: one MSSM point is considered 
distinguishable from another one when the difference between values of at least one of the observables for 
the two points is greater than 5 standard deviations of the considered observable cje , calculated for the point 
we take as reference, i.e. if, for example, Nj^f — Nj^^ > 5ae'^-^ where Nj^f is the number of jets for the 

reference point, Nj^^ is the number of jets for another points and al'^-^ is the statistical uncertainty on the 
average number of jets of the reference point. By calculating, for each MSSM point, the difference between 
the MSSM reference point and any other point for each observable and by expressing these variations in terms 
of the value of respective uncertainties calculated for the reference point, we can discriminate between them. 
Figure [I^ shows an example of discrimination for the following values of parameters: 

Mj = 3000 GeV, Mi = 100 GeV, M2 = 1000 GeV, Mg = 2700 GeV, M-g = 2400 GeV, 
Ma = 200 GeV, /i = 500 GeV, ^3 = GeV, tan /? = 50. (9) 
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b) cut 1 + 2 
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Fig. 10: Separation of MSSM points (each entry corresponding to a set of MSSM parameter values): a) sets remaining after one cut ( 
{niep) variation < Scxe); b) sets remaining after two cuts ( (rijet) variation < 5cTe); c) sets remaining after three cuts ( EY''^"" variation 
< ScTe ) and d) sets remaining after four cuts ( £1"'" variation < 5cre )• 



By applying successively a cut at 5 standard deviations on the following observables: average number of 
leptons, of jets, missing transverse energy and total transverse energy, we manage to separate the reference 
point from all the other ones. This method works for all MSSM points ( the number of cuts needed for the 
separation varies from point to point). 

4.23 Evaluation of the statistical uncertainties 

We showed in the previous section that there is a possibiUty to distinguish one set of parameter values from 
the others. We are now going to estimate the statistical uncertainties for each MSSM parameter value. We 
estimate these uncertainties using the uncertainties calculated for each observable. For that, we consider one 
set of MSSM parameter values defined as the reference point. Then for any chosen parameter, we measure the 
number of standard deviations for each observable between the reference point and the MSSM point having the 
same parameter values except for the one in question. The value of this parameter has to be different by one 
unit on the grid of values. For example, if we want to obtain the statistical uncertainties on the /x parameter for 
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1 If* 

V aiLiu^ 


/T 1 


U — 


M(L) (GeV) 


3000. 


9.899 


9.899 


Ml (GeV) 


500. 


3.035 


4.259 


Mz (GeV) 


500. 


5.451 


5.451 


M(Q) (GeV) 


1800. 


12.384 


21.092 


M(G) (GeV) 


1800. 


2.481 


1.901 


A (GeV) 


200. 


4.749 


4.749 


tan [3 


50. 


0.575 


0.575 


fi (GeV) 


2000. 


17.331 


17.331 


A3 


2000. 


8.534 


8.534 



Table 5: Resolutions of the MSSM parameter values for intermediate masses M(g) and M(g) (large statistics) and significance = 30. 



fixed values of the other parameters, we consider the point which has the same values for the other parameters 
and the next higher value for /i in the grid of parameters and then we calculate 



N 



I rpMiss rpMiss I 



(10) 



is the missing transverse energy at the point we want to calculate the resolution, E^2^^ the 

the statistical uncertainty estimate for £^|;^****. We 



where E^{'-'' 

missing transverse energy at the other point and a^nns 
assumed that 



N 



(11) 



l-n 



where a^-^ is the statistical uncertainty we want to estimate. We could also take another observable or a linear 
combination of observables to calculate the uncertainties. 

In the following two tables, we have used EJ^^^^ as the observable for the calculation of the resolution. 



parameter 


value 


0-+ 




M(L) (GeV) 


1000. 


146.491 


146.491 


Ml (GeV) 


2000. 


55.387 


55.387 


M2 (GeV) 


500. 


23.740 


23.740 


M(Q) (GeV) 


2700. 


18.868 


16.289 


M(G) (GeV) 


2700. 


42.142 


40.814 


A (GeV) 


3000. 


118.390 


118.390 


tan/3 


50. 


4.411 


4.411 


IX (GeV) 


500. 


136.504 


136.504 


A3 


2000. 


95.730 


95.730 



Table 6: Resolution of the MSSM parameter values for high masses M(q) and M(3) (low statistics) and significance = 6. 

Tables |5| and ^ show some examples of resolution obtained by this method. These values are small 
since this is only the statistical error. The resolution should degrade after including systematic errors. Table ^ 
shows the calculation of statistical resolution for medium squark and gluino masses, where statistics are large. 
Table ^ shows a similar calculation for large masses for squarks and gluinos, where statistics are low and the 
uncertainties are thus much more important. 
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Fig. 11: An example of a resolution in the plane Mq versus Mg (the curve represents 1 a). 



Figure 1 1 shows an example of resolution in the plane mg versus rrig for the following reference point: 



Mj = 200 GeV, Mi = 100 GeV, Ms = 500 GeV, Mg = 1500 GeV, Mg = 2100 GeV, 
Ma = 200 GeV, /x = 200 GeV, ^3 = 2000 GeV, tan /3 = 50. 



(12) 



We have, for this point, a resolution of about 5 to 10 GeV (using the number of events with leptons, Nqi, as 
the discriminator). 



5. Conclusion 

We demonstrated the possibility to discover a phenomenological MSSM using an inclusive study in the MSSM 
parameters space. Once we have discovered SUSY, using the kinematical observables for parameter determi- 
nation proved to be an efficient method. Statistical uncertainties obtained are relatively small (< 50 GeV for 
squark and gluino masses). We could note, at the end of our study, that there was little difference between 
mSUGRA and the pMSSM. The discovery limit is to a large extent determined by the total cross section (and 
is around 2.7 TeV at CMS). The only difference appears for some points having a specific mass hierarchy. As 
an example, in the case of a compact hierarchy of masses, the limit we expect is about 1.5 TeV. 



58 



Part VIII 

Establishing a No-Lose Theorem for NMSSM Higgs 
Boson Discovery at the LHC 

U. Ellw anger, J.F. Gunion, C. Hugonie 

Abstract 

We scan the parameter space of the NMSSM for the observability of at least one 
Higgs boson at the LHC with 300 fb^^ integrated luminosity, taking the present 
LEP2 constraints into account. We restrict the scan to those regions of parameter 
space for which Higgs boson decays to other Higgs bosons and/or supersymmetric 
particles are kinematically forbidden. We find that if M/^VF-fusion detection modes 
for a light Higgs boson are not taken into account, then there are still significant 
regions in the scanned portion of the NMSSM parameter space where no Higgs bo- 
son can be observed at the ba level, despite the recent improvements in ATLAS and 
CMS procedures and techniques and even if we combine all non-fusion discovery 
channels. However, if the 1^1^-fusion detection modes are included using the cur- 
rent theoretical study estimates, then we find that for all scanned points at least one of 
the NMSSM Higgs bosons will be detected. If the estimated 300 fb~^ significances 
for ATLAS and CMS are combined, one can also achieve Scr signals after combin- 
ing just the non-H^iy -fusion channels signals. We present the parameters of several 
particularly difficult points, and discuss the complementary roles played by different 
modes. We conclude that the LHC will discover at least one NMSSM Higgs boson 
unless there are large branching ratios for decays to SUSY particles and/or to other 
Higgs bosons. 



1. Introduction 

Supersymmetric extensions of the standard model generally predict relatively light Higgs bosons. One of the 
most important tasks of the LHC is the search for Higgs bosons [79,^ 114]. An important milestone in 



understanding the potential of the LHC was the demonstration that at least one Higgs boson of the minimal 
supersymmetric standard model (MSSM) would be detectable at the > So" level throughout all of the MSSM 
parameter space so long as top squark masses do not exceed 1.5 to 2 TeV and so long as large branching 
fractions to decay channels containing supersymmetric particles are not substantial. 



In [115], we studied, subject to these same and a few other simplifying restrictions, the detectability 



of Higgs bosons in the next-to-minimal supersymmetric standard model (NMSSM). This short note presents 



the most relevant procedures and conclusions of [115]. In the NMSSM, one Higgs singlet superfield, S, is 



added to the MSSM in order to render unnecessary the bilinear superpotential term 11H1H2 by replacing it 
with XSH1H2, where the vacuum expectation value of the scalar component of S, (S), results in an effective 
bilinear Higgs mixing with = X{S). The detectability of the NMSSM Higgs bosons was first considered in 
a contribution to Snowmass 96 [ |1 16[ ]. The result, using the experimentally established modes and sensitivities 
available at the time, was that substantial regions in the parameter space of the NMSSM were found where 
none of the Higgs bosons would have been observable either at LEP2 or at the LHC even with an integrated 
luminosity of 600 fb^^ (two detectors with L = 300 fb^^ each). 



59 



Since then, progress has been made both on the theoretical and the experimental sides. On the theoretical 



side, the dominant two-loop corrections to the effective potential of the model have been computed [117 



118]. These lead to a modest decrease in the mass of the lightest Higgs scalar, holding fixed the stop sector 



pai^ameters. Inclusion of the two-loop corrections thus increases somewhat the part of the NMSSM parameter 
space excluded by LEP2 (and accessible at the Tevatron) [US], but is of less relevance for the LHC. On the 
experimental side the expected statistical significances have been improved since 1996 [^^,[n|]. Most 
notably, associated tih production with h ^ bb (originally discussed in [ |119| ]), which in the SM context is 
particularly sensitive to m/j < 120 GeV, has been added by ATLAS and CMS to the list of Higgs boson detection 
modes [[79|,^, 114 , 12C-123|. Analysis of this mode was recently extended [124] to rrih = 140 GeV, which, 
though not relevant in the SM case due to the decline in the bb branching ratio as the WW* mode increases, 
is highly relevant for points in our searches for which the WW* mode is suppressed in comparison to the SM 
prediction. In addition, techniques have been proposed [125- |130| ] for isolating signals for WW fusion to a 
light Higgs boson which decays to rr or WW^*\ 

It turns out that adding in just the tih process renders the no-Higgs-discovery parameter choices de- 
scribed and plotted in [ 116 |, including the "black point" described in detail there, visible [1311- I^i 1 115], we 
searched for any remaining parameter choices for which no Higgs boson would produce a > 5cr signal. In this 
search, we performed a scan over nearly all of the parameter space of the model, the only parameter choices 
not included being those for which there is sensitivity to the highly model-dependent decays of Higgs bosons 
to other Higgs bosons and/or superparticles. The outcome is that, for an integrated luminosity of 300 fb~^ at 
the LHC, there are still regions in the parameter space with < 5a expected statistical significance (computed as 
NsD = S/ \fB for a given mode) for all Higgs detection modes so far studied in detail by ATLAS and CMS, 
i.e. including the tih — > iihb mode but not the VFl^-fusion modes. On the other hand, the expected statistical 
significance for at least one of these detection modes is always above 3.6(T at 300 fb~^, and the statistical signif- 
icance obtained by combining (using the naive Gaussian procedure) all the non-T^H^-fusion modes is at least 
4.8(7. However, we find that all such cases are quite observable (at > 10. la) in one of the T^H^-fusion modes 
(using theoretically estimated statistical significances for these modes). For all points in the scan of parame- 
ter space, statistical significances obtained by combining all modes, including l^VF-fusion modes, are always 
>10.7(T. Thus, NMSSM Higgs discovery by just one detector with L = 300 fb~^ is essentially guaranteed for 
those portions of parameter space for which Higgs decays to other Higgs bosons or supersymmetric particles 
are kinematically forbidden. This represents substantial progress towards guaranteeing LHC discovery of at 
least one of the NMSSM Higgs bosons. 

In order to clarify the nature of the most difficult points in those portions of parameter space considered, 
we present, in sect. 4, examples of particularly difficult bench mark points for the Higgs sector of the NMSSM. 
Apart from the "bare" parameters of the model, we give the masses and couplings of all Higgs scalars, their 
production rates and branching ratios to various channels (relative to the SM Higgs) and details of the statistical 
significances predicted for each Higgs boson in each channel. The latter will allow an assessment of exactly 
what level of improvement in statistical significance will be required in the various different detection modes 
in order to render marginal modes visible. Of course, our estimates of the expected statistical significances 
are often somewhat crude (e.g. their dependence on the accumulated integrated luminosity). We believe that 
our procedures always err in the conservative direction, leading to statistical significances that might be a bit 
small. Thus, the LHC procedures for isolating Higgs boson signals could provide even more robust signals for 
NMSSM Higgs boson detection than we estimate here. 

The detection modes, which serve for the searches for standard model or MSSM Higgs bosons, include 
(using the notation h, a for CP-even, CP-odd Higgs bosons, respectively): 

^) 99 ^ 77; 

2) associated Wh or tih production with 77^^ in the final state; 
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3) associated tth production with h bb; 

A) gg ^ h/a or associated bbh/a production witli h/a tt; 

5) gg^h^ ZZ^*'> 4 leptons; 

6) gg^h^ WW^*'^ l+l~uV; 

7) LEP2 e+e^ Zh and e^e~ —>■ ha; 

8) WW -^h^TT; 

9) WW -^h^ WW^*\ 

where 8) and 9) are those analyzed at the theoretical level in [125- |130| ] and included in the NMSSM analysis 



for the first time in this paper. The above detection modes do not employ the possibly important decay channels 
i) /i — > hh, ii) h — > aa, iii) h —>■ h^h^, iv) h —>■ aZ, w) a ha, vi) a — > hZ, vii) h,a h^W^, 
viii) h,a ^ ti and ix) t h^b. The decay modes i)-vii) give high multiplicity final states and deserve a 



dedicated study [132], while the existing analyses of the tt final state signatures are not very detailed. Further, 
when kinematically allowed, the t —>■ h'^b signal would be easily observed according to existing analyzes. 
Thus, in this paper we restrict our scan over NMSSM parameter space to those parameters for which none of 
these decays are present. In addition, we take the constraints of LEP2 [via the mode 7)] into account, and only 



accept points for which 5a discovery at LEP2 would not have been possible [133, 134] . 

The Higgs sector of the NMSSM consists of 3 scalars, denoted /ii, /12, h^ with ruh^ < nih^ < nih^, 2 
pseudo-scalars, denoted ai, 02 with nia^ < maj, and a charged Higgs pair, denoted h^. Mixing of the neutral 
doublet fields with the gauge singlet fields in the scalar and in the pseudo-scalar sector can be strong. The scalar 
mixing can lead to a simultaneous suppression of the couplings of all the hi to gauge bosons, and hence to a 
suppression of many of the detection modes above. (Of course, the have no tree-level couplings to gauge 
boson pairs and the one-loop couplings are too small to yield useful event rates.) The couplings of the Higgs 
bosons to t- or b-quarks can be amplified, reduced or even change sign with respect to the standard model 
couplings. Hence negative interferences can occur among the (loop-) diagrams contributing to gg — > hi and 
hi 77, leading again to suppressions of the above detection modes. A complete simultaneous annulation of 
all detection modes is not possible, but simultaneous reduction of all detection modes is possible and it is for 
such parameter choices that NMSSM Higgs boson discovery is most difficult. 

In the next section, we define the class of models we are going to consider, and the way we perform 
the scan over the corresponding parameter space. In section 3 we describe our computations of the expected 
statistical significances of the detection modes 1) - 9) above. In section 4, we present six particularly difficult 
bench mark points (in table [l]) and details regarding their statistical significances in channels l)-9) in table |^ 
with a summary of overall statistical significances in table ^. Using these tables, we give a discussion of the 
properties of these points. 

2. NMSSM Parameters and Scanning Procedure 



In this paper, we consider the simplest version of the NMSSM [|135|-|48|], where the term ^HiH2 in the 
superpotential of the MSSM is replaced by (we use the notation A for the superfield and A for its scalar 
component field) 

\H1H2S + p3 ^ 

so that the superpotential is scale invariant. We make no assumption on "universal" soft terms. Hence, the five 
soft supersymmetry breaking terms 

m]j^Hl + ml^Hl + mlS^ + \AxH1H2S + '^A^S'^ (2) 
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Point Number 


1 


2 


3 


4 


5 


6 


Bare Parameters 




A 


0.0340 


0.0450 


0.0230 


0.0230 


0.1330 


0.0230 


K 


0.0198 


0.0248 


0.0129 


0.0069 


0.1459 


0.0114 


tan/? 


6.00 


5.25 


-5.5 


5.75 


-8 


-6 


/Ueff(GeV) 


140 


-110 


115 


-235 


100 


150 


Ax{GeV) 


-35 


25 


-95 


40 


-135 


-100 


A^iGeV) 


-150 


70 


-90 


80 


-75 


-110 


Scalar VTasses and Coimlinffs 








1 MM 


1 M'i 


1 1 '2 


1 1 A 


1 IZ 


Cv 


-u.oo 


M '20 


M '2/1 


M ^1 
U.O / 


M OT 
-U.O / 


M "71 
-U. / 1 


Ct 


-U.OJ 


U.jU 




M ^< 
U.OJ 


M Q1 

-U.oi 


M 

-U.OO 


Cb 


-i .U / 


M 

u.oo 


-L.Z / 


1 1 A 
1.10 


A ^M 

-4. jU 


O /IM 

-Z.4U 


gg Production Rate 


f\ 'JO 


M MO 


M MO 


M '20 

yj.jy 


M <A 
U.JO 


U.JO 


BR-n 




u.zo 


n no 
u.uy 


U.JO 


U.Uj 


nil 

U.l 1 


BRbh = BRtt 


1 1 o 


1 MQ 

i.Uo 


1 1 M 

l.iU 


1 1 o 


110 

l.lo 


11^ 

1.1 J 


BRWW^*^ 






OR 

\J.\JO 


^7 

V/. J / 


OA 

\J.\J'-r 




m J, rfieVI 

//t/l2 V,VJC V ^ 


125 


114 


114 


126 


144 


122 


Cv 


-0.74 


-0.83 


0.79 


-0.73 


0.46 


0.59 


Ct 


-0.72 


-0.74 


0.70 


-0.71 


0.57 


0.54 


Ch 


-1.49 


-3.28 


3.46 


-1.47 


-6.66 


2.24 


gg Production Rate 


0.46 


0.44 


0.40 


0.45 


1.18 


0.23 


BR^n 


0.33 


0.08 


0.07 


0.34 


0.01 


0.10 


BRbb = BRtt 


1.30 


1.18 


1.18 


1.32 


3.06 


1.31 


BRWW^*^ 


0.32 


0.08 


0.06 


0.33 


0.01 


0.09 




205 


153 


148 


201 


202 


155 


Cv 


-0.14 


-0.46 


-0.51 


-0.15 


0.18 


-0.39 


Ct 


-0.30 


-0.63 


-0.67 


-0.32 


0.17 


-0.55 


Cb 


5.80 


4.17 


4.20 


5.53 


0.68 


5.12 


gg Production Rate 


0.31 


0.84 


0.95 


0.33 


0.02 


0.80 




0.13 


0.05 


0.05 


0.15 


0.98 


0.03 


BRhh = BRtt 


308.66 


5.83 


3.92 


274.41 


13.97 


8.12 


BRWW^*1 


0.18 


0.07 


0.06 


0.21 


0.96 


0.05 


Pseudo-Scalar Masses and Couplings 




m„, (GeV) 


191 


112 


130 


130 


113 


145 


Ct 


0.03 


-0.03 


-0.10 


-0.01 


-0.10 


-0.16 


Cb 


1.16 


-0.83 


-2.95 


-0.19 


-6.55 


-5.77 


gg Production Rate 


0.00 


0.00 


0.03 


0.00 


0.31 


0.08 


ma, (GeV) 


206 


141 


137 


198 


174 


158 


Ct 


0.16 


0.19 


-0.15 


0.17 


-0.07 


-0.05 


Cb 


5.89 


5.18 


-4.64 


5.75 


-4.59 


-1.65 


gg Production Rate 


0.02 


0.07 


0.06 


0.02 


0.03 


0.00 


Charged Higgs Mass 




mc(GeV) 


221 


162 


157 


213 


157 


167 



Table 1 : We tabulate the input bare model parameters, the corresponding Higgs masses, and the corresponding Higgs couplings, relative 
to SM Higgs boson coupling strength, for 6 bench mark points. Also given for the CP-even hi are ratios of the gg production rate and 
various branching fractions relative to the values found for a SM Higgs of the same mass. For the CP-odd Ui, "gg Production Rate" 
refers to the value relative to what would be found if both the bb and the it 75 couplings had SM-like strength. 
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Point 


1 


2 


3 


4 


5 


6 


Channel 






hi Higg 


s boson 






Nsd{1) 


3.74 


0.35 


0.13 


3.18 


0.62 


0.83 


NsDi2) 


4.37 


0.59 


0.22 


3.92 


0.85 


1.22 


Nsoi^) 


2.79 


0.85 


0.85 


3.03 


4.83 


3.30 


Nsd{4) 


0.08 


0.07 


0.76 


0.09 


4.52 


0.40 


Nsd{5) 


0.83 


0.00 


0.00 


0.64 


0.12 


0.16 


Nsd{Q) 


1.10 


0.09 


0.03 


0.90 


0.16 


0.22 


Nsd{7) 


0.00 


3.37 


3.40 


3.29 


0.00 


4.79 




9.29 


Y.ll 


1.59 


8.93 


16.78 


1 A AO 

10.08 




z..jy 


no 

W. WW 


w.wu 


1 lA 


zLI 


zLQ 




6.54 


1.09 


1.17 


5.99 


6.69 


3.65 




fi 54 


3 55 


3 59 


6 84 


6 69 

\j.\jy 


6 0? 




11.61 


1.64 


1.97 


10.89 


18.07 


10.73 




1 1 . 1 


J . 1 J 




1 1 . Jl o 


1 o.W / 


1 1 . / J 


Channel 






/i2 Higg 


s boson 






Nsd{1) 


3.69 


0.83 


0.61 


3.62 


0.22 


0.55 




4.01 


1.25 


0.92 


3.93 


0.05 


0.74 


NsoiS) 


2.49 


3.95 


3.58 


2.30 


0.99 


1.77 


Nsd{4) 


0.16 


2.76 


2.93 


0.16 


3.62 


2.99 


NsDi5) 


1.84 


0.16 


0.11 


1.94 


0.56 


0.20 


Nsd{G) 


1.44 


0.22 


0.16 


1.46 


0.38 


0.18 


Nsd{7) 


0.00 


0.00 


3.31 


0.00 


0.00 


0.00 


Nsd{8) 


15.39 


15.17 


13.46 


15.05 


7.41 


9.89 




^ 7Q 


W.U J 


AA 


fi 0^ 


1 Q 

W. i " 


89 




6.44 


5.05 


4.76 


6.31 


3.82 


3.61 




6 44 


5 05 


5 80 

^ .0\J 


6 31 


3 82 


3 61 




17.65 


16.00 


14.28 


17.40 


8.34 


10.56 




1 / .Oj 


1 f\ no 

10. WW 


1 '4.00 


1 7 40 
1 / .^U 


O. 


IU.jO 


Channel 






hz Higg 


s boson 






Nsd{1) 


0.00 


0.59 


0.66 


0.01 


0.00 


0.32 


Nsd{2) 


0.00 


0.21 


0.25 


0.00 


0.00 


0.08 


Nsni^) 


0.00 


0.00 


1.13 


0.00 


0.00 


0.00 




3.79 


3.43 


3.62 


3.56 


1.55 


4.86 


Nsd{5) 


3.65 


2.51 


2.07 


4.46 


1.54 


1.66 


Nsd{6) 


0.80 


2.13 


1.52 


1.17 


0.38 


1.55 


Nsd{7) 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


NsoiS) 


0.00 


0.00 


9.06 


0.00 


0.00 


0.00 


NsDi9) 


0.00 


0.77 


0.79 


0.00 


0.00 


0.43 




5.32 


4.80 


4.64 


5.83 


4.76 


5.37 




5.32 


4.80 


4.64 


5.83 


4.76 


5.37 


V'E.=1-6,8,9[^Sd(*)]2 


5.32 


4.86 


10.21 


5.83 


4.76 


5.39 




5.32 


4.86 


10.21 


5.83 


4.76 


5.39 



Table 2: Scalar Higgs statistical significances, Nsu = S/\/B, in various channels for the 6 bench mark points. For each individual 
Higgs, we give (in order): Nsd for the channels 1) - 9) described in the text; Gaussian combined Nsd for non-lWF-fusion LHC 
channels; combined Nsd for non-M^M^-fusion LHC chaimels plus LEP2; combined Nsd for all LHC chaimels, including the fusion 
channels WW h ^ tt and WW — » /i — » IVW*^*^ channels; and combined Nsd for all LHC channels plus LEP2. 
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Point Number 


1 


2 


3 


4 


5 


6 


Best non-M^W^ fusion Nsd 


4.37 (hi) 


3.95 {h2) 


3.62 (hs) 


4.46 (/i3) 


4.83 (/ii) 


4.86 (hs) 


Best WW fusion Nsd 


15.39 (/i2) 


15.17 (/i2) 


13.46 (/i2) 


15.05 (/i2) 


16.78 (/ii) 


10.08 (/ii) 


Best combined Nsd w.o. 
M^H^-fusion modes 


6.54 (hi) 


5.05 (/i2) 


4.76 ih2) 


6.31 (/i2) 


6.69 (hi) 


5.37 (/13) 


Best combined Ns d with 
M^M^-fusion modes 


17.65 ih2) 


16.00 (/i2) 


14.28 (/i2) 


17.40 (/i2) 


18.07 (/ii) 


10.73 (hi) 



Table 3: Summary for all Higgs bosons. The entries are: maximum non-WH^ fusion LHC Nsd', maximum LHC WW fusion Nsd', 
best combined Nsd after summing over all non-M^W-fusion LHC channels; and best combined Nsd after summing over all LHC 
charmels. The Higgs boson for which these best values are achieved is indicated in the parenthesis. One should refer to the preceding 
table in order to find which channel(s) give the best values. 



are considered as independent. The masses and/or couplings of sparticles are assumed to be such that their 
contributions to the loop diagrams inducing Higgs production by gluon fusion and Higgs decay into 77 are 
negligible. In the stop sector, which appears in the radiative corrections to the Higgs potential, we chose the 
soft masses mq = rriT = Mgusy = 1 TeV, and varied the stop mixing parameter 



,2 



As in the MSSM, the value Xt = \/6 - so called maximal mixing - maximizes the radiative corrections to the 
Higgs masses, and we found that it leads to the most challenging points in the parameter space of the NMSSM. 

Assuming that the Higgs sector is CP conserving, the independent parameters of the model are thus: 
X,K,m'jj ,m'jj^,m'^, Ax and A^. For purposes of scanning and analysis, it is more convenient to eliminate 
m'jj^, ^11^ "^1 favor of Mz, tan/3 and fieS = ^{S) through the three minimization equations of the 



Higgs potential (including the dominant 1- and 2-loop corrections [118|) and to scan over the six independent 
parameters 

\,K,tanP,fie{i,Ax,Af,. (4) 

We adopt the convention X,k > 0, in which tan /3 can have either sign. The absence of Landau singularities 
for A and k below the GUT scale (~ 2 x 10^^ GeV) imposes upper bounds on these couplings at the weak 



scale, which depend on the value of ht and hence of tan/3 [ |135| - |139| ]. Using mf"'*^ = 175 GeV, one finds 



Amax ~ 0.69 and Kmax ~ 0.62 for intermediate values of tan /?. 

For each point in the parameter space, we diagonalize the scalar and pseudo-scalar mass matrices and 
compute the scalar, pseudo-scalar and charged Higgs masses and couplings taking into account the dominant 1- 



and 2-loop radiative corrections [1 18|. We then demand that the Higgs scalars satisfy the LEP2 constraints on 



the e'^e Zhi production mode (taken from [ |133| ], fig. 10), which gives a lower bound on nih^ as a function 
of the ZZhi reduced coupling. We also impose LEP2 constraints on e+e~ — > hiaj associated production (from 



[ |134| ], fig. 6), yielding a lower bound on m/j. + as a function of the Zhiaj reduced coupling. 

In order to render the above-mentioned processes i) - ix) kinematically impossible, we require the fol- 
lowing inequalities among the masses: 

m/j3 < 2mhi, 2mai, 2mft±, m^^ + Mz, mh± + Mw] 
"i-aa < iTT-hi + m-ai, ""^hi + Mz, mj^± + Mw', rnj^± > 155GeV. 
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In addition we require |/iofr| > 100 GeV; otherwise a light chargino would have been detected at LEP2. 
(The precise lower bound on | jj^cs \ depends somewhat on tan (3 and the precise experimental lower bound on 
the chargino mass; however, our subsequent results do not depend on the precise choice of the lower bound 
on l/icffl-) We further note that for the most challenging parameter space points that we shall shortly discuss, 
l^offl > 100 GeV is already sufficient to guarantee that the NMSSM Higgs bosons cannot decay to chargino 
pairs so long as the SU(2) soft-SUSY-breaking parameter M2 is also large. In fact, in order to avoid significant 
corrections to 77/14 and 770^ couplings coming from chargino loops it is easiest to take M2 ^ ^cfr (or vice 
versa). This is because the hixfXi coupling is suppressed if the xt is either pure higgsino or pure gaugino. 
Since the parts of parameter space that are challenging with regard to Higgs detection typically have ~ 
100 — 200 GeV, the validity of our assumptions requires that M2 be large and that the chargino be essentially 
pure higgsino. 



Using a very rough sampling, we determined, as expected from previous work [116], that it is only for 
moderate values of tan [3 that < 5(T signals might possibly occur. From this sampling, we determined the most 
difficult parameter space regions and further refined our scan to the following: 

• 4.5 < I tan/3| < 8 (both signs) in steps of 0.25; 

• 0.001 < A < min[0.21, Amax]. using 20 points; 

• 0.001 < K < min[0.24, Kmax]. using 20 points; 

• 100 GeV < |/ieff| < 300 GeV (both signs), in steps of 5 GeV; 

• < |Aa| < 160 GeV, with A\ opposite in sign to /icg, using steps of 5 GeV; 

• 25 GeV < \A^^^ < 170 GeV, with ^4^ opposite in sign to /iefr> using steps of 5 GeV. 

For those points sampled in this final scan which satisfy all the constraints detailed earlier, we compute the 
expected statistical significances for the processes 1) to 9) listed in section 1, as described in the next section. As 
a rough guide, from the ~ 10^ points detailed in the above list, we find about 250, 000 that pass all constraints 
and have Nsd < 5 (for L = 300 fb^^) in each of the individual discovery modes 1) - 7). We shall tabulate a 
number of representative points taken from this final set in section 4. 

3. Expected Statistical Significances 

^From the known couplings of the NMSSM Higgs scalars to gauge bosons and fermions it is straightforward 
to compute their production rates in gluon-gluon fusion and various associated production processes, as well 
as their partial widths into 77, gauge bosons and fermions, either relative to a standard model Higgs scalar or 
relative to the MSSM H and/or A. This allows us to apply "NMSSM corrections" to the processes 1) - 9) 
above. 

These NMSSM corrections are computed in terms of the following ratios. For the scalar Higgs bosons, 
cy is the ratio of the coupling of the hi to vector bosons as compared to that of a SM Higgs boson (the coupling 
ratios for hiZZ and hiWW are the same), and q, Cb are the corresponding ratios of the couplings to top and 
bottom quarks (one has Cr = c^). Note that we always have |cy| < 1, but q and C}, can be larger, smaller or 
even differ in sign with respect to the standard model. For the CP-odd Higgs bosons, cy is not relevant since 
there is no tree-level coupling of the Oj to the VV states; ct and are defined as the ratio of the 275 couphngs 
for ti and bb, respectively, relative to SM-like strength. 

We emphasize that our procedure implicitly includes QCD corrections to the Higgs production processes 
at precisely the same level as the experimental collaborations. First, the ATLAS and CMS collaborations 



employed Monte Carlo programs such as IS AJET [|98|] and PYTHIA [ |113| ] in obtaining results for the (MS)SM. 
These programs include many QCD corrections to Higgs production in a leading-log sense. This is the best 
that can cuiTcntly be done to implement QCD corrections in the context of experimental cuts and neural-net 
analyses. Clearly the more exact NNLO results for many of the relevant processes will slowly be implemented 
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in the Monte Carlo programs and increased precision for Higgs discovery expectations will result. Since our 
goal is to obtain NMSSM results that are completely analogous to the currently available (MS)SM results, we 
have proceeded by simply rescaling the available (MS)SM experimental analyses. In doing the rescaling of the 
Higgs branching ratios we have included all relevant higher-order QCD corrections [ ]149| ] using an adapted 
version of the FORTRAN code HDECAY [ |150| ]. Details regarding our rescaling procedures can be found in 
1 115 1 . Using the rescaling procedures, for each point in the parameter space of the NMSSM we obtain the 
statistical significances predicted for an integrated luminosity of 100 fb^^ for each of the detection modes 1) - 
9). In order to obtain the statistical significances for the various detection modes at 300 fb^^, we multiply the 
100 fb^^ statistical significances by \/3 in the cases 1), 2), 3), 5) and 6), but only by a factor of 1.3 in the cases 
4), 8) and 9). That such a factor is appropriate for mode 4), see, for example, fig. 19-62 in [[79[]. Use of this 
same factor for modes 8) and 9) is simply a conservative guess. 



4. Difficult Points 

As stated in the introduction we still find "black spots" in the parameter space of the NMSSM, where the 
expected statistical significances for all Higgs detection modes 1) - 7) are below 5cr at 300 fb'^ The reasons 



for this phenomenon have been described above; see also the corresponding discussion in [ |116| ]. However, 
after including the modes 8) and 9), the points that provide the worst 1) - 6) statistical significances typically 
yield robust signals in one or the other of the VFVF-fusion modes 8) and 9). 

In order to render the corresponding suppression mechanisms of the detection modes reproducible, we 
present the detailed properties of several difficult points in the parameter space in table [1} The notation is as 
follows: The bare parameters are as in eq. (2.5), with rrijj^, and m| fixed implicitly by the minimization 
conditions. (As noted earlier, with the convention A, k > in the NMSSM, the sign of tan /3 can no longer 
be defined to be positive.) For the reasons discussed below eq. (|3|) we chose in the stop sector uiq = uit = 
Msusy = 1 TeV and Xt = VQ for all of the points (1-6). We have also fixed m^op*^ = 175 GeV. For both scalar 
and pseudoscalar Higgs bosons, "gg Production Rate" denotes the ratio of the gluon-gluon production rate with 
respect to that obtained if q = Cb = 1, keeping the Higgs mass fixed. For scalar hi, this is the same as the ratio 
of the gg production rate relative to that predicted for a SM Higgs boson of the same mass. For the scalar hi, 
BR'y^ denotes the ratio of the 77 branching ratio with respect to that of a SM Higgs boson with the same mass. 
(A verification of the reduced gluon-gluon production rates or 77 branching ratios would sometimes require 
the knowledge of the couplings to higher precision than given, for convenience, in table |l|.) Also given for the 
scalar hi are the ratios BRbb and BRWW* of the bb and WW* branching ratios relative to the SM prediction 
(as noted above, one has BRtt = BRbb). 

In table ^, we tabulate the statistical significances for the hi in all the channels 1) - 9); production of the 
CP-odd Oj turns out to be relevant only when they add to the hi signals in process 4). Also note that, all these 
problematical points are such that m/j^ + nia^ > 206 GeV, so that e+e" — > /ii + oi followed by /ii, ai bb 
would have been kinematically forbidden at the highest LEP2 energy. Hence, for LEP2 mode 7) we only give 
the statistical significance for e^e~ —>■ Zhi. Also tabulated in tableware four statistical significances obtained 
by combining various channels. This combination is done in the Gaussian approximation: 



AT-combined 



1/2 



where ^ • runs over the channels i being combined. 

a) NsD obtained by combining LHC channels 1) - 

b) NsD obtained by combining LHC channels 1) - 



We give results for the following combinations: 
6); 

6) and LEP2; 
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c) NsD obtained by combining LHC channels 1) - 6) with the W^l^-fusion channels 8) and 9); 

d) NsD obtained by combining all LHC channels and LEP2, i.e. by combining all channels 1) - 9). 

In those cases where there is no LEP2 signal, a)=b) and c)=d). In addition, in our point selection we have 
required a mass difference of at least 10 GeV between scalar Higgses, so that they yield well separated signals 
and no statistical significance combination of two different scalar Higgses is needed. All parameter choices 
for which Higgs boson masses differ by less than 10 GeV yield stronger signals than the cases retained. (The 
increased net signal strength of overlapping Higgs signals in those channels with limited mass resolution arises 
as a result of Nfj^{l + 2) ~ {Si + S2)/7b > y'Sf + S'i/VB.) 

As summarized in table |3[ all of the tabulated "bench mark points" have statistical significances below 
5(7 for all of the detection modes 1) - 6) at 300 fb^^ and 7) at LEP2. In more detail, as tabulated in table ^ and 
summarized in table ^, the best signals in the modes 1) - 6) for the points #1 - #6 at the LHC are: 

• point #1, NsD =4.37 for mode 2) and hi; 

• point #2, NsD =3.95 for mode 3) and /12; 

• point #3, NsD =3.62 for mode 4) and h^; 

• point #4, NsD =4.46 for mode 5) and /13; 

• point #5, NsD =4.83 for mode 3) and hi, 

• point #6, NsD =4.86 for mode 4) and h^; 

Further, for point #3, the combined statistical significance of modes 1) - 6) (also tabulated in table ^ would 
still be below 5 for any one hi, although \/2Ng^ > 5 (as is hkely to be relevant by combining ATLAS and 
CMS data once each detector has accumulated L = 300 fb~^) for at least one of the /ij. However, for all these 
"difficult" points the l^l^/'-fusion modes 8) and/or 9) provide (according to theoretical estimates) a decent 
(sometimes very strong) signal. 

The points #1 - #4 differ as to which of the modes 1) - 6) and which hi yields the largest statistical 
significance should the l^VK-fusion mode 8) not provide as strong a signal as suggested by the theoretical 
estimates. To render these points observable without the l^/'VF-fusion mode 8) would require improvements of 
all detection modes 2) - 5). 

As in [116], we find that difficult points in the parameter space generally have | tan/3| ~ 5. This is 
the region of tan (3 for which the bbh, bba signals are still not very much enhanced but yet the gg —>■ h,a and 
tih, tia signals have been suppressed somewhat. In a few cases, however, difficulties also arise for | tan f3\ as 
large as 8, as shown in the case of point #5. Also as in [ ]116| ], the most difficult points are those in which the 
masses of the hi and are relatively close in magnitude, typically clustered in a ~ 60 GeV interval above 
~ 105 GeV. Such clustering maximizes the mixing among the different Higgs bosons and thereby minimizes 
the significance of the discovery channels for any one Higgs boson. In particular, it is for strong mixing among 
the hi that the statistical significance for discovery modes based on a large VV coupling for any one hi are 
most easily suppressed. 

Finally, for point #6, we have minimized the statistical significances for the TV VF-fusion modes over 
the parameter space, while keeping the statistical significances of modes 1) - 6) below 5. One can see that it 
still gives a strong lO.lcr signal in mode 8). [Smaller Nsd for mode 8) would have been possible if we had 
allowed stronger signals in modes 1) - 6), in particular had we allowed smaller mass separation, < 10 GeV, 
between the two lightest Higgs bosons.] In addition, for point #6 = 112 GeV and the ZZ coupling of hi 
is sufficiently large that it would have yielded a A.8a signal at LEP2. Had we taken a top quark mass slightly 
larger, mf^p*^ = 178 GeV, we would have found a very similar point with a hi mass of ~ 115 GeV, which could 
have been responsible for the excess observed at LEP2 [151-154]. 
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5. Discussion and Conclusions 

In this paper, we have addressed the question of whether or not it would be possible to fail to discover any 
of the Higgs bosons of the NMSSM using combined LEP2 and LHC data, possibly resulting in the erroneous 
conclusion that Higgs bosons with masses below 200 GeV have been excluded. We have demonstrated that, 
assuming that the decay channels i) - ix) are either kinematically disallowed or render a Higgs boson observ- 
able, this is unlikely (at the > 5a level) to happen. Certainly, there are points in NMSSM parameter space for 
which the statistical significances for the individual detection modes 1) - 6) (i.e. those analyzed in detail by 
ATLAS and CMS) are all well below So" for integrated luminosity of 300 fb^^. However, by combining several 
of the modes 1) - 6) and 300 fb^^ data from both ATLAS and CMS, a > 5o" signal can be achieved based 
just on modes 1) - 6). Further, we have found that throughout all of the NMSSM parameter space (scanned 
subject to the earlier listed restrictions) for which such weak signals in modes 1) - 6) are predicted, the theo- 
retical estimates for the VFl^/'-fusion modes indicate that an easily detected WW —>■ h ^ tt signal should be 
present. Thus, our conclusion is that for all of the parameter space of the NMSSM compatible with reasonable 
boundary conditions for the parameters at the GUT scale (with, of course, non-universal soft terms in general) 
and such that Higgs pair and SUSY pair decays of the Higgs bosons are kinematically forbidden, at least one 
of the NMSSM Higgs bosons will be detected at the LHC. This is a big improvement over the results from the 
earlier Snowmass 1996 study which was somewhat negative without the inclusion of the tih tibb mode 3), 
and the VFVF-fusion modes 8) and 9). 

It is amusing to note that all of our bench mark points for which Higgs discovery is most difficult at the 
LHC include a Higgs scalar with mass close to 115 GeV (with, however, reduced couplings to the Z boson). 



which could be responsible for the excess observed at LEP2 [|151j-|154p. 

Another important point that appears from our analysis is the fact that the full L = 300 fb^^ of inte- 
grated luminosity (per detector) is needed in order to have robust NMSSM Higgs discovery in the portion of 
parameter space considered here. Of course, as in the MSSM, it is very possible that only one of the CP-even 
NMSSM Higgs bosons might be detected at the LHC but that, as studied by Kamoshita et al. in []140|-|148[], the 



observation of all the CP-even Higgs bosons of the NMSSM would be possible at the LC by virtue of all having 
some non-negligible level of ZZ coupling and not having very high masses. Even at the LC, the CP-odd Higgs 
bosons might escape discovery, although this would not be the case for the parameter choices that we have 
found which make LHC discovery of even one NMSSM Higgs bosons most challenging. This is because, for 
such parameters, the Oj are relatively light and could be readily seen at the LC in the processes e~^e~ — > /ijOj, 
e+e^ —>■ vVaiai and e^e~ ^ Z* Zaiai, assuming an integrated LC luminosity of 1000 fb~^ and energy 



^/s> 500 GeV |155| 



This study makes clear the importance of continuing to expand the sensitivity of existing modes and 
continuing to develop new modes for Higgs detection at the LHC in order not to have to wait for construction 
of a linear e+e^ collider for detection of at least one of the SUSY Higgs bosons. In particular, study of modes 
i) - ix) and SUSY pair channels should all be pushed. The problematical points that we have emphasized 
here are unlikely to be substantially influenced by it or SUSY decays since all the Higgs masses are below 
~ 200 GeV so that it decays will be kinematically highly suppressed (one of the top quarks would have 
to be virtual) and SUSY pair decays are quite unlikely to be significant given LEP2 limits on the masses of 
SUSY particles. However, by allowing Higgs (in particular, pseudoscalar) masses such that one or more of the 
channels i)-vii) are kinematically allowed we have found points for which discovery in modes l)-9) will not be 



possible [132]. Thus, a full "no-lose" theorem for NMSSM Higgs boson discovery at the LHC will require 
exploring additional discovery modes sensitive to those portions of parameter space for which Higgs decays 
to other Higgs bosons are important, and might necessitate combining results from both the ATLAS and CMS 
detectors and/or accumulating more integrated luminosity. 



68 



6. Acknowledgments 

We wish to acknowledge helpful discussions with M. Sapinski and D. Zeppenfeld. C.H. would like to thank the 
late theory group of the Rutherford Appleton Laboratory and the Theoretical Physics Department of Oxford, 
where part of this work was achieved, for their kind hospitahty. This work was supported in part by the U.S. 
Department of Energy. 



69 



Part IX 

Effects of Supersymmetric Phases on Higgs 
Production in Association with Squark Pairs in the 
Minimal Supersymmetric Standard Model 

A. Dedes, S. Moretti 

Abstract 

We show how the Supersymmetric (SUSY) CP-violating phases can induce new final 
states in associated production of Higgs bosons with squark pairs of identical flavor 
(for the A^) as well as modify substantially those already present when the soft SUSY 
parameters are real (in the case of and h^). Hence these processes, particularly 
for light stop squarks, ii, are good candidates for phenomenological investigation, 
in order to confirm or disprove the existence of complex soft SUSY parameters. We 
illustrate this in the context of a general Minimal Supersymmetric Standard Model 
(MSSM), for a choice of SUSY parameters accessible at the Large Hadron Collider 
(LHC). 



It has recently been shown [ 156| -162] that, if CP-violating effects are manifestly inserted into the MSSM 
Lagrangian, by allowing the Higgsino mass term, /i, and the trilinear couplings, to be complex, thereby 
introducing two independent CP- violating phases [ 163 , |164 ], cj)^ and 0^, such that e**^^ = /u/|/i| and e**^-* = 
then the strength of the Higgs couplings to (s)particles can drastically be modified, inducing sizable 
effects, e.g., in the dominant production mode of neutral Higgs bosons at the LHC, i.e., gg (where 
= , and A^), through the squark-squark-Higgs vertices involving stops and sbottoms. These effects 
are a consequence of large values attained by (l)n and/or (j)A consistent with cancellations taking place in the 



SUSY contributions to the Electric Dipole Moments (EDMs) of neutron and electron [165-168] 



These same interactions also affect the associated production of Higgs bosons with third generation squark 



pairs. Since this process is expected to be accessible at the LHC, see Refs. [ ]169| - |175p , we investigate here some 
aspects of its phenomenology in the presence of complex parameters in the MSSM Lagrangian. Schematically, 
the production mechanism is the following: 



(1) 



where q = t,b, x^'^ = 1,2 and = /i*^, H'^, A'^, in all possible combinations, as appropriate in the MSSM. 
Notice that in such processes the existence of CP-violating effects in the SUSY Lagrangian would immediately 
be manifest from the detection of three particle final states involving a pseudoscalar Higgs boson and two 
identical squarks. In fact, if (p^ = (pA = 0, even in presence of mixing between the third-generation squarks, 
the q-f^Qy^A'^ couplings, with q = t,b and x = 1> 2, are identically zero [171]. Depending on the relative value 



of the final state masses in ([l]), rrig^, m^^, and M$o, the production of Higgs particles can be regarded as taking 
place either via a (anti)squark decay or via a Higgs-strahlung. 



^For simplicity, we assume A = Au ~ Ad at the electroweak (EW) scale, i.e., 0{A4z), where it and d refer to all flavors of up- 
and down-type (s)quarks. 
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We work in the theoretical framework provided by the 'complex' MSSM, the latter including explicitly the two 
CP- violating phases, (p^ and (pA, and assuming universality of the soft gaugino masses only at the Grand Unifi- 
cation (GUT) scale. We define its fundamental parameters without making any assumptions about the structure 
of the SUSY breaking dynamics at the Planck scale, whether driven by Supergravity (SUGRA), gauge/anomaly 
mediated (GMSB/AMSB) or proceeding via other mechanisms as we treat the MSSM as a low-energy effective 
theory. Among the possible setups of the MSSM parameter space compatible at one-loop with the EDM data, 
we choose here the one presented in Table |l|, in terms of fi, A, M^o (the mass of the pseudoscalar Higgs state), 
tan (3 (the ratio of the vacuum expectation values of the two Higgs doublet fields), Mq^ ^ 3 (the soft squark 
masses of the three squark generations) and Mg (the soft gluino mass). For this choice of MSSM inputs, the 
derived masses of the scalar are barely consistent with the latest bounds on Mffi from LEP2, of about 91 



GeV [ ]176| ] {mi is instead set above the Tevatron limit from Run 1, around 140 GeV [ ]177[ ] for our tan j3). This 
means that our specific parameter point will readily be probed at the LHC, as in correspondence of the inputs 
given in Table [l] and by varying both cf)^ and 4>a between and vr, one has that the lightest Higgs mass spans 
from 91 to 100 GeV and the hghtest stop one is between 140 and 240 GeV or so. Notice that the MSSM setup 
given in Table |l| should be regarded as one possible example of the rich phenomenology that can be induced by 
the CP- violating phases in the MSSM with a low tan /3 and squark/Higgs masses small enough to be produced 
at detectable rates via process d). In fact, we have found many others but refrained from showing them here 
for reasons of space. 

In the remainder, we will denote the regions of the ((/>^, 4)a) plane excluded from Higgs and squark direct 
searches by a shaded area. In addition, the inputs in Table [I] comply with the constraints deduced from the two- 
loop Barr-Zee type contributions to the fermionic EDMs [ ]178| , |179| ] (green "<&" symbols in the following) for 
most choices of (p^j^ and (pA- Finally, some {cpf^, (Pa) points will further be neglected following the requirement of 
positive definiteness of the squared squark masses: see eqs. (5)-(6) of Ref. [ 158| - 162] (magenta "x" symbols 



in the forthcoming plots). The top-left corner of Fig. [I] shows the phase dependence of M^o and m^^, outside 
such experimentally excluded areas. 

In the course of our discussion, we shall make only one simplification, which will not alter the conclu- 
sions of our work. That is, we will neglect one-loop mixing effects among the three neutral Higgs bosons of 



the MSSM [ ]18(]| - |185| ], on the ground that for our choice of parameters they turn out to be of the order of a few 
percent at most (see discussion in Refs. [156 -[162 1). Indeed, much larger effects will remain unaccounted for, 
such as the higher-order QCD corrections to the production process (|l]), which are likely to induce i^T-factors 
of the order 1.5-2 and whose calculation is presently not available. 











Mao 


tan P 


600 


2500 


300 


1000 


200 


2.7 



Table 1: One possible parameter setup of the MSSM satisfying the one-loop EDM constraints (all quantities in GeV, apart 
from the dimensionless tan (3) and yielding cross sections for process ([I]) manifestly dependent on the CP-violating phases 
(t)u and (j)A. 



We start our numerical analysis by referring to Fig. 1 of Ref. [ 156| , 157 |, where one can find the contour 
plots for the minimum values of the modulus of the common trilinear coupling, \ A\, above which the mentioned 
cancellations work at a level which is compatible with the experimental accuracy achieved in both the neutron 
and electron EDM measurements, i.e., \dri 



< 6.3 X IQ-^*^ ecm [jlS^] and \de\ < 4.3 x 10"^^ ecm 01871 ], 
for the choice of parameters in Table I and as a function of cp^ and cpA- These A values are those entered in 
the production vertices of the processes we considered, alongside the two discrete quantities and tan /3. As 
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mj-j^ and mjj^o. tanjS=3.7 
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Fig. 1: Contour plots illustrating the 0^ and (j)A dependence of: the lightest Higgs and stop masses (top-left), in red and 
blue, respectively; the cross sections for gg, qq fif*(/i°)[i/°]{74°} (top-right) [bottom-left] {bottom-right}, in red. The 
meaning of the black shaded area and of the colored green/magenta symbols is given in the text. For the cross sections, 
we have used the CTEQ(4L) [190| Parton Distribution Functions with Q = ^/l as the factorisation scale (also in a^), the 
latter evolved at two loops with all relevant (s)particle thresholds onset within the MSSM (as described in [p3|,p4|]). 



for the Higgs masses, we keep Af^o fixed at 200 GeV and derive tlie values of M/^o and Mj^o at two-loop level 
(see the discussion in Refs. [ 156 -162|] concerning the residual theoretical error on the latter). 

Among the processes of the type ([T]), much emphasis has been put on the case in which both the squark 
and Higgs scalar states are the lightest, that is, on the mechanism gg,qq iii\h^ [|l69|,|l7|]. In fact, if the 
typical A scale is in the TeV regime, then two concurrent effects take place, that render light Higgs production in 
association with the lightest scalar top quarks a more favorable Higgs discovery channel than the corresponding 
SM-like one, gg, qq — > tihP [188, 189]. On the one hand, since the mixing angle 9^ is proportional to rritA, 
the ii squark becomes much lighter than the t-quark and all other squarks. On the other hand, by looking at the 
expression of the iiiih^ vertex (e.g., see eq. (3) of Ref. [169|), for large A values, it is clear that its strength 
can overcome that of the top-Higgs Yukawa coupling entering the SlM-like reaction. 

We thus continue our investigation of squark-squark-Higgs production by considering this particular final 
state. Our choice of tan P in Table I reflects the remark made in Refs. 1 169, that the iiiih^ production 
rates are larger at smaller tan/3 (see, e.g.. Fig. 3 of Ref [169]), this being the consequence of the increase 
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of M^o for large values of the latter [158-162]. (As representative of the low tan/3 regime, the two values 
tan /? = 2, 3 were used in Refs. [169, 170].) As for the A spectrum, we are bound to those values guaranteeing 



the EDM cancellations (again, see Fig. 1 in [ ]156| , |157[ ]): that is, between 100 and 700 GeV, depending on 
(j)f^ and (j)A- The top-right corner of Fig. [I] presents the cross sections for gg, qq —>■ tit^/i*^ as a contour plot 
over the (cp^, (/>^) plane. We see a strong dependence on the SUSY phases, as the production rates vary over 
several orders of magnitude. The maximum of the cross section occurs at large (pA values, when 0^ is slightly 
above 7r/2. This can be understood in the following terms. Here, it is where the lightest stop mass reaches 
its allowed minimum. However, the fact that does not grow similarly at small 0^ values, when (p^ 

is slightly below tt/2 - the other region of the (c/)^, cpA) plane where is minimal - implies that the iiiih^ 
vertex too plays an important role in determining the actual size of the cross section. The fact that this coupling 
is maximal where the stop mass is minimal can easily be understood by noticing the cos cp^ dependence of 
eqs. (5) and (A.53) of Ref. [ |158 - 162] when (pA — > vr. Besides, a large value of the iiiih^ coupling combined 
with a small value of implies that the two-loop contributions of the Barr-Zee type graphs to the EDMs 
can be sizable, so that it is not surprising to see that larger rates for gg, qq — > iiilh^ accumulate towards the 
correspondingly excluded area. However, in areas not yet removed through the EDM measurements, one can 
find production rates as large as 800 fb. Finally, notice that - for the choice of parameters in Table I and 100 
GeV <A< 700 GeV as in Fig. 1 of [ Il56i[r57| ] - in the 'phaseless' hmit, i.e., the standard MSSM case, 
(pi_i,cpA — > 0, the yield is several orders of magnitude smaller, indeed well below detection level. 

A strong hierarchy exists, M/^o <C Mfjo Mao, among the neutral Higgs masses, for our setup of 
the MSSM. This should naturally allow one to disentangle in the experimental samples iiilh^ from iii\H^ + 
iii^A^ events. Furthermore, under the assumption that the lightest scalar quark will promptly be discovered at 
the LHC from some other source of SUSY events than those in (|l|), and its mass measured, we believe that final 
states with two heavy objects of identical mass m^^ recoiling against a rather central one with mass M^o (that, 
again, we assume to be known and reconstructed through bb and/or 77 decays), should be distinguishable from 
other Qx^x'*^" channels (q = t, b), in which X or x' 7^ iri the transverse mass distributions of the visible 

ii and i\ decay products. 

Even in such circumstances though, iii\H^ and tii\A^ events would still obey a sort of degeneracy 
in their appearance (especially after accounting for large detector resolutions in reconstructing masses), that 
could render non-trivial the operation of separating data containing one type from those induced by the other. 
(Our choice of MSSM parameters, producing M^a ~ M^o, ought to be representative also of such extreme 
experimental conditions.) Nonetheless, we see that the maxima and minima of a{iii*iH^) and a{iii\A^) occur 
in very different regions of the i(p^,cpA) plane. Moreover, in the allowed areas, the difference between the 
production rates of the two processes can even be a factor of 10 or more. In other terms, although it could 
always be possible to attempt the above 'separation' on the basis of the different decay patterns of the two 
Higgs bosons (and/or their topology), this might not be needed after all. In fact, as long as fi, A and tan (3 
have been constrained to some extent through some other measurements, then for some specific values of , 
a clear excess (i.e., well above the size of the uncertainties induced by unknown higher order QCD effects) 
of tit^^>° events, with 200 GeV < M^o < 210 GeV, above the = rates, could only be explained if 
iiilA^ events have indeed been produced. 

By comparing the bottom-left to the bottom-right plots in Fig. |l|, one realises that this separation can 
happen over a large portion of the (cp^, (pA) plane. One should also note the very different shapes of the two 
contours, such that tit^A*^ rates are largest where the iii*H^ ones are smallest. For example, just outside the 
areas excluded by the EDMs (when cp a or tt and (p^ w 7r/2), the pseudoscalar Higgs channel can reach 
the 100 fb level, whereas the scalar Higgs rates are always around 10 fb. For — > and (p^ 7r/6 (or, 
quite symmetrically, for cpA ^ it and cp^ ^ tt — 7r/6), the two process rates are of the same order, about 10 fb. 
Finally, in the limit cp^, (pA — > 0, a{tii\H^) is about 10 fb and a{iii\A^) is, of course, zero. 
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In our numerical simulations, we also have considered the large tan (3 scenario. However, in this case, 
once the EDM constraints were taken into account, we have found that both the theoretical plausibility of the 
MSSM and the phenomenological impact of the CP-violating phases were much reduced. On the one hand, 
in order to obtain the mentioned cancellations also for tan /3 > 10, one would need to have the soft squark 
masses Mq^ ^ as large as 6 TeV or more and the gluino one Mg > 3 TeV, that is, a quite 'unnatural' hierarchy 
in the soft SUSY breaking sector, if one aims to maintain Mg.^ around 300 GeV (so that final states 

remain within the reach of the LHC energy). On the other hand, for large tan (3 values, the not yet excluded 
(by direct searches, two-loop effects in the EDMs and positive definiteness of the squark masses) area of the 
((/>^, (j)A) plane is much smaller (see Fig. 3 of Ref. [158-162]). (Besides, also mixing effects among neutral 
Higgs states start becoming relevant for tan f3 > 10). Thus, although some sporadic points over the allowed 
{4>pi,4>A) regions can still be found, these yielding cross sections significantly different from those obtained 
in the phaseless case, we would conclude that only the 2 < tan /3 < 10 region is relevant in the experimental 
analysis of squark-squark-Higgs production. 

To summarise, we have shown that the LHC production rates of the lightest Higgs boson of the MSSM 
in association with a pair of lightest stop scalars are strongly affected by the presence of complex parameters 
in the soft sector of the SUSY Lagrangian, even when EDM constraints are taken into account. As a matter of 
fact, the gg, qq — > iiilh^ mechanism has recently been advocated as a new possible discovery mode of the 
boson, at least for certain combinations of the MSSM parameters, that we have emulated here to some extent. 
Thus, our results in this case have a twofold meaning. On the one hand, they emphasise that more inputs than 
those pertaining to a phaseless MSSM could be needed to describe the phenomenology of iii\h^ events (further 
recall that we have limited ourselves to the case of only two independent phases, (j)^ and (pA, those associated to 
the Higgsino mass and the universal trilinear couplings, respectively). On the other hand, they make the point 
that such a mechanism can be useful in assessing whether or not soft CP-violating phases are present. 

In this last respect, however, it would be even more intriguing to detect final states involving the pseu- 
doscalar Higgs boson in place of the lightest scalar one at the LHC. In fact, no matter the actual setup of the 
MSSM parameters, the detection of a tit\A^ state would unequivocally mean that (j)^ and/or (I)a are finite. In 
fact, the corresponding interaction is prohibited at tree-level in a MSSM with real masses and couplings in the 
soft sector and EW effects at the one-loop level are unlikely to yield tii\A^ production rates as large as those 
shown here: up to 10000 events per year for some phase combinations at high collider luminosity. Further 
notice that, in our analysis, we deliberately have chosen a small value of tan (), so that mixing effects among 
the three neutral Higgs states are very small even when the values of Mfjo and M^o are rather close. Finally, 
despite this mass degeneracy between the heavy scalar and pseudoscalar Higgs bosons, the relative production 
rates of iii\H^ and iii\A^ events are very different, both in size and in shape, over most of the (0^, 0^) plane, 
so that the two samples could even be separated experimentally. 
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PartX 

Study of the Lepton Flavor Violating Decays of 
Charged Fermions in SUSY GUTs 

T. Blazek 

Neutrino oscillations clearly show that individual lepton flavor is violated in nature. Here we present results 
of a study of lepton flavor violating decays £^ ^ 7 in a class of SUSY GUT models assuming that the 
third right-handed neutrino couples equally to the second and third lepton doublets with a large coupling. This 
corresponds to the neutrino Yukawa matrix (in the left-right basis) of the form 



1^ ~ 1 . (1) 




The large 23 entry, responsible for the large atmospheric neutrino mixing angle, induces potentially large lepton 
flavor violating effects in the low-energy effective theory. The see-saw mechanism yields a physical neutrino 
with a mass about 5 x 10^^ eV consistent with the Super Kamiokande observation providing the third right- 
handed neutrino mass is A//j3 w 3 x 10^^ GeV, much greater than the masses of the other two right-handed 
neutrinos. 

We focus on the large tan/3 50 regime of these models with (1^6)33 = (^i/)33- The rate of a decay 
7 is enhanced by (tan /?)^ and can be, approximately, rescaled by this factor for lower values of 
tan/3. Thus our study provides for the upper estimate of the lepton flavor violating decay rates of charged 
fermions. In Figure la-b we present the results for the branching ratio r /i7 in a typical model of this 
class. Besides (|I|) all other Yukawa matrices are hierarchical with small off-diagonal entries (more details 
can be found in |jT9l[]). The model assumes the Pati-Salam S'C/(4) x SU{2)l x SU{2)r symmetry at the 
GUT scale Mq ~ 3 x 10^^ GeV but the results in the leptonic sector would be quite similar for models with 
different unifying gauge group leading to Eq. (|i|). The two plots were obtained for two different values of the 
yu parameter fixed to 120 GeV in plot (a) and 300 GeV in plot (b). is the soft universal scalar mass and 
Ml/2 is its gaugino analogue, both at Mq- The t —>■ contour lines in Figure la-b should be compared to 
the experimental upper bound BR{t + 7) < 1-1 • 10"^. The allowed region overlaps with the region 
preferred by the data on the muon — 2, as shown in plots Ic-d. Clearly, this decay rate should be very close 
to the present limits and thus presents an exciting opportunity to observe lepton flavor violation or to constrain 
substantially this class of models. 

The decays /i — > e 7 and r — > e 7 are much more model dependent since small Yukawa entries are 



necessarily involved. In the model studied in [191] they are found to be well below the experimental limit. 
Thus our main result is the correlation between the maximal atmospheric mixing angle and large r — > /i 7 
branching ratio related through the large off-diagonal entry in the Yukawa matrix in (|I|). 
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Fig. 1: Contour lines of BR{t ^7) x lO** and Sa^,{SUSY) x lO" for two different values of fi. (a) BR{t -> fij) x 10**, for 
H = 120GeV. (b) BR{t ^i-f) x 10^ for ^J, = 300GeV. (c) 5a^,{SUSY) x 10^°, for = 120GeV. (d) 5a^,{SUSY) x 10^^ 
for /i = SOOGeV. In (c) and (d) the dashed curves mark the central value for not accounted for by the Standard Model and the 
borderlines of 1-sigma region for this quantity. The experimental upper limit on BR{r /17) is 1.1 x 10^. In all plots tan/3 — 50 
and soft trilinear parameter A — 0. 
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Part XI 

Interactions of the Goldstino Supermultiplet with 
Standard Model Fields 

D.S. Gorbunov 

Abstract 

In a set of supersymmetric extensions of the Standard Model the masses of the sgold- 
stinos are of the order of the electroweak scale. Thus sgoldstinos are expected to 
be produced at future colliders. The sgoldstino interactions with the fermions and 
gauge bosons of the Standard Model are determined by the MSSM soft mass terms 
and the scale of supersymmetry breaking. These interactions have been included 
into the CompHEP package. On the other hand, the sgoldstino couplings to Higgs 
bosons depend on the parameters of the hidden sector responsible for mediation of 
supersymmetry breaking. The measurement of these coupling constants would offer 
a unique probe of the hidden sector. 



1. Goldstino supermultiplet 

In any supersymmetric extension of the Standard Model (SM) of particle physics spontaneous supersymmetry 
breaking occurs due to a non-zero vacuum expectation value of an auxiliary component of some chiral or vector 
superfield. As a simple case, let us consider a model where 

S = s + \/26'V' + O'^Fs (1) 

is the only chiral superfield which obtains a non-zero vacuum expectation value F for its auxiliary component, 

{Fs) = F. (2) 

Then ip is a two-component Goldstone fermion, goldstino, and its superpartners 

are respectively a scalar and a pseudoscalar sgoldstino. 

In the framework of supergravity becomes the longitudinal component of the gravitino, due to the 
super-Higgs effect. As a result, the gravitino acquires a mass which in realistic models with a vanishing 
cosmological constant is completely determined by the supersymmetry breaking parameter F: 



VStt F 

The sgoldstinos remain massless at tree level and become massive due to corrections from higher order terms 
in the Kahler potential. If these terms are sufficiently suppressed, the sgoldstinos are light and may appear 
in particle collisions at high energy colliders. Such pattern emerges in a number of non-minimal supergravity 



models 1 192, 193] as well as in gauge mediation models if supersymmetry is broken via a non-trivial superpo- 



tential (see, e.g. [|194|,[195|] and references therein). Here we shall consider the sgoldstino masses ms and nip 
as free parameters. 
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The gravitinos and sgoldstinos interact with the MSSM fields and the corresponding coupling constants 
are inversely proportional to the supersymmetry breaking parameter F. Eq. (0) then implies that the gravitino 
has to be very light, otherwise the gravitino as well as the sgoldstinos are effectively decoupled from the MSSM 
fields at the energy scale of the colliders of the neai-term future. 



2. Effective Lagrangian 



The effective Lagrangian for the gravitino is obtained from N=l supergravity [ ]196| ], and may be used to 
calculate scattering processes involving any of the helicity components of the massive gravitino. Meanwhile the 
energy scale attainable at the present and the nearest future generation of accelerators favors the longitudinal 
component of gravitino as the most promising to be studied in collision experiments. In this case the longitu- 
dinal component of the gravitino, (5^ ~ z9^V'/"^3/2' effectively behaves as a two-component fermion if), the 
goldstino, and the interaction between the goldstino and the other fields is plainly given by the Goldberger- 
Treiman relation 

with J SUSY being a supercurrent. 

In order to obtain the low-energy effective Lagrangian for sgoldstinos one can use the spurion 



method [197]. It exploits the fact that, by definition, sgoldstinos are scalar components of the very super- 
multiplet S whose auxiliary component acquires a non-zero vacuum expectation value due to supersymmetry 
breaking. Then one can consider a simple supersymmetric model with non-renormalizable interactions be- 
tween S and the MSSM superfields, which yield the MSSM soft terms when a non-zero {Fs) is generated. 
Consequently, the corresponding coupling constants are fixed by the ratios of the soft terms and F. 

In general, the supersymmetry breaking part of the Lagrangian has the form 

- CbreaHng = J] "^fc I '^^ I ' + ( ^ ^« TrA°A° + h.C. J 
k \ a J 

- {BWoh^^ + A^JielWj, + A^^qidlWj, + A^^g^ng/i^^ + h.c) , (5) 

where k {a) runs over all scalar 0^ (gaugino Xq) fields. Since supersymmetry is broken spontaneously, Eq. 
implies the following effective interaction between the MSSM superfields and the goldstino supermultiplet S: 

^S—MSSM ^S~Kdhler ~l~ ^S— gauge ~\~ ^S— superpotential > 



where 



Cs-KaHler = - j cP 9 dH S ■ J] || $t g^^i , 

all matter 
and Higgs fields 

Cs-gauge = ^ [d^O S ^TrWW" + h.C. , 



all gauge 
fields 



f f B A^ ■ A^ ■ ■\ 

^S- superpotential = Id 9 S ■ €ij{— H\)HIj + V^ElH\) + QlD^H}j + Q^U^HIj \ + h.C. 

These terms emerge if the fields from the hidden sector, where supersymmetry breaking occurs, are integrated 
out. The only remnant of the hidden sector is the goldstino supermultiplet, which may remain light. 
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Integrating over 6, 9 and taking into account Eqs. (|l]), @ we obtain the soft supersymmetry breaking 
terms as well as the interactions between the components of the goldstino supermultiplet and the components 
of the MSSM superfields: 



Cs-Kdhler = - i^S -^iFk+h.C.] , (6) 

all matter ^ ^ 

and Higgs fields 

all gauge ^ 



fields 



^"-s ■ F:,,F: - i^s ■ F^^^e^'^'PF:,^ + h.c)j , (7) 



- - a ^lu- a 8F 



B , , B 



C-S- super potential — ' XhXjJ ~ • yh^j^F^j^ + F}j^h-[ 



+ -f [liet ■ shh + lixh ■ set + sli ■ etxh " sFlelh^j, - sliFE^h^o " slielFlj^ 



^ A 



D 



F 



+ Q>t • shjj + sqi • ulxu + sul ■ qixu " sF^Qulh^^ - sqiFu^hfj - s^^ulF^^ 



+ h.c.j . (8) 

Here we presented only the leading order in terms; the convention for the Levi-Civita tensor is e^iss _ _-|^ 

Eliminating the auxiliary fields, one obtains the low-energy effective Lagrangian for the interactions 
between the components of the goldstino supermultiplets and the MSSM fields. 



3. Phenomenology of the model 

Until now, there is no experimental evidence for a gravitino or sgoldstinos. The study of their phenomenology 
places bounds on their coupling constants. Note that all sgoldstino coupling constants introduced in the previous 
sections are completely determined by the MSSM soft terms and the supersymmetry breaking parameter F, 
while the sgoldstino masses ms and mp remain arbitrary. Depending on the values of ms and mp, the 
sgoldstinos may show up in different experiments. The phenomenologically interesting models can be separated 
into four classes: 

• The sgoldstino masses are of order the electroweak scale, while \/F ~ 1 TeV — sgoldstinos may then be 



produced at high-energy colliders [ |198| , |199| ] (see section 3 of Ref . [ ]200| ] for a sketch of the sgoldstino 
collider phenomenology). 

The sgoldstino masses m^, mp ~ 1 MeV-^1 GeV, while \/F ~ 1 TeV — sgoldstinos may then emerge 
as products of rare meson decays [201, 202], such as T ^ S{P)^, J/tp —>■ S{P)^. 

Models with flavor violation in the soft trilinear couplings, Aab / A6ab — sgoldstino interactions then 
lead to flavor violating processes. In particular, sgoldstinos may contribute to FCNC (mass differences 
and/or CP- violation in the neutral meson systems) | 203[ , |204j ], and, if kinematically allowed, sgoldstinos 
appear in rare decays such as t ^ cS{P) [[205[], fi —>■ eS{P), K —>■ ttS [202], etc. 



The sgoldstinos are lighter than 1 MeV — these models may be tested in low energy experiments [ ]202| ], 
such as reactor experiments, conversion in a magnetic field, etc. Sgoldstinos may also play a very impor- 



79 



tant role in astrophysics and cosmology [202, 206-208]: they may change the predictions of Big Bang 



Nucleosynthesis, distort the CMB spectrum, affect SN explosions and the cooling rate of stars, etc. 
4. Incorporation into CompHEP 

To leading order in 1/F and to zero order in the MSSM gauge and Yukawa coupling constants, the interactions 
between the components of the golds tino supermultiplet and the MSSM fields are derived in Ref. [ ]200| ]. They 
correspond to the processes most attractive for collider studies - where only one of these new particles appears 
in a final state. In this case the light gravitino behaves exactly as a goldstino. For sgoldstinos, as they are R-even, 
only the sgoldstino-goldstino and sgoldstino-SM fields couplings have been included as the most interesting 
phenomenologically. All new coupling constants between the components of the goldstino supermultiplet and 
the MSSM fields are completely determined by the ratios of the soft supersymmetry breaking parameters and 
F. This Lagrangian has been incorporated (see Ref. [ 20C ] for details) into the CompHEP software package [J 



This package may be used in the calculation of any tree-level process with one on-shell gravitino or sgoldstino. 
The universality of the Lagrangian makes it possible to apply the package in studying the phenomenology 
of any supersymmetric extension of the Standard Model. Being included into CompHEP, this model should 
be regarded as an additional option allowed within the framework of the CompHEP/SUSY model. Currently 
accessible versions of this package operate only with real parameters and coupling constants. Likewise the 
trilinear soft couplings are assumed to be proportional to the corresponding Yukawa couplings, Aab = ■AabUab- 
The sgoldstino Lagrangian transformed in accordance with these rules reads 

M 

= - E ^S.F-^,Fr''-^y^fS{e^jlieth}, + h.c.) 

all gauge ^ ^ 
fields 

ad aU 

- ^yf, • S{e,,qidthh + h.c.) - -^y^, ■ S{e,,qythlj + h.c.) , (9) 

M 

= H ^P-F^,,e>^'''^F:,^-i-^y!^,.P{e.,lieth'n-h.^^^ 

all gauge ^ ^ 
fields 

- i^Vab ■ P{e^jqidth^D " h.c.) - i^vl ' P{^^l<ta<K " ^.c) . (10) 



5. Remarks 

It is worth noting that the independent direct measurement of the MSSM soft supersymmetry breaking terms 
and the gravitino or/and sgoldstino couplings offers a unique possibility to estimate the supersymmetry breaking 
scale \fF. 

The sgoldstino couplings to superpartners become relevant for models with heavy sgoldstinos, where the 
sgoldstinos also decay to SM superpartners. In the spurion approach the corresponding coupling constants are 
not completely determined by the MSSM soft breaking terms, but may depend on new parameters originating 
from the hidden sector. A similar situation happens with the sgoldstino interaction terms proportional to the 
MSSM gauge or Yukawa couplings. For instance, accounting for the F and D auxiliary fields (see Eqs. (^, 
(0), (|8|)) may cause sgoldstino-Higgs mixing. The corresponding coefficients also depend on new parameters 
from the hidden sector. Indeed, the additional interaction between the sgoldstino and the Higgs bosons arises 



The CompHEP package [209] automatically calculates tree-level particle decay rates and cross sections and is aimed to improve 



the accuracy, to cut down the efforts and to shorten the time usually required for studying high-energy collision processes. 
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from the effective superpotential 

C = Pl J(feS^ ■ eijH}jHij + h.c. 

with jl being a new dimensionful constant. Thus the sgoldstino couplings to Higgs bosons allow us to probe 
the hidden sector. The sgoldstino decay mode into the lightest MSSM Higgs boson becomes important in 
models with fairly Ught sgoldstinos (ms(^p) > 250 GeV) and the measurement of this partial width constrains 
some combination of hidden sector parameters, if the supersymmetry breaking scale is known (i.e., from the 
sgoldstino partial width into two photons and a measurement of the gaugino masses). 
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Part XII 

Attempts at Explaining the NuTeV Observation of 
Di-Muon Events 

A. Dedes, H. Dreiner, and R Richardson 

Abstract 

The NuTeV Collaboration has observed an excess in their di-muon channel, possibly 
corresponding to a long-lived neutral particle with only weak interactions and which 
decays to muon pairs. We show that this can not be explained by pair production of 
neutralinos as suggested in the literature. In the parameter region allowed by LEP the 
event rate is far too small. We propose instead a new neutralino production method 
via S-mesons, which can fully explain the observation. 

1. Introduction 

The NuTeV Collaboration has searched for long-lived neutral particles (N^) with mass M^o > 2.2 GeV and 



small interaction rates with ordinary matter [ ]10| , |210| , [211[ ]. They look for the decay of the neutral particles in 



a detector which is 1.4 km downstream from the production point and observe 3 /i/i events where they only 
expect to see a background of 0.069 ± 0.010 events. The probability that this is a fluctuation of this specific 
channel is about 4.6 a. The simple supersymmetric scenarios discussed previously can not lead to an excess at 
NuTeV, since the decisive supersymmetric parameter range has been excluded by LEP. We propose instead the 
production of light neutralinos via i?-mesons which could give a measurable excess. We briefly present the two 
possible models and then discuss them quantitatively. We have also considered the production rate for neutral 



heavy leptons but this is too low and does not lead to a viable explanation [212] 



2. The Rp Violating Model 

Only couplings, A212 and A232 give a di-muon signature. For A212 the neutralino will decay with equal probabil- 
ity to efiu and fifiu. No e/i-events are observed, we therefore propose one dominant R-parity violating coupling 
A232- This coupling corresponds to the two neutralino decay modes Xi f^Lf^n'^r and x? '^l ^^^^ 
as their complex conjugates. For a light neutralino the r/i decays are sufficiently phase space suppressed to give 
an expectation below one event. For the light neutralino production we shall only consider single neutralino 
production in the decay of bottom hadrons. The bottom hadrons are formed following the production of a bb 
pair. These hadrons can then decay via the R-parity couplings A'-^g, (i = 1, 2, 3). We will only consider the 
decays of the and via R-parity violation i.e. — > z/jX? and B^ — > ^fxi- This mechanism allows 
one to produce light neutralinos via a strong interaction process. 

Using results for the RPV branching ratios of the B mesons and the neutralino lifetime we can find 
regions in (A2325 ^113) parameter space, for a given sfermion mass, in which there are 3 ± 1 events inside the 
NuTeV detector, this is shown in Fig. [T[ 

This model can be tested at other experiments. At the NOMAD experiment for the same B^-meson 
branching ratio we obtain about an order of magnitude more events than at NuTeV. Thus our model can be 
completely tested by the NOMAD data! For neutralino production we are relying on a rare B -meson decay 
which can possibly be observed at a present or future B-factory although this may be difficult as the leptons 
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Fig. 1: Regions in A232, K13 parameter space in which we would expect 3 ± 1 events to be observed in the NuTeV detector. The limits 
[213] on the couplings A'113 (crosses) and A213 (diamonds) allow solutions between the two points. The region above the stars is ruled 
out for the coupling A213 by the limit on the product of the couplings A232A213. The hatched region shows the bound on the coupling 
Aii3 from perturbativity. 



produced will be very soft. 
3. Conclusions 

We have reconsidered the NuTeV di-muon observation in the light of supersymmetry with broken R-parity 
and neutral heavy leptons. We have shown that it is not possible to obtain the observed event rate with pair 
production of light neutralinos or via the production of neutral heavy leptons. However, we have introduced 
a new production method of neutralinos via B-mesons. Due to the copious production of 5-mesons in the 
fixed target collisions the observed di-muon event rate can be easily obtained for allowed values of the R-parity 
violating couplings. 

The model we have proposed can be completely tested using current NOMAD data. If the NOMAD 
search is negative our model is ruled out and the NuTeV observation is most likely not due to physics beyond 
the Standard Model. 
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Part XIII 

Kaluza-Klein States of the Standard Model Gauge 
Bosons: Constraints From High Energy Experiments 

K. Cheung and G. Landsberg 

Abstract 

In theories with the standard model gauge bosons propagating in TeV~^-size extra 
dimensions, their Kaluza-Klein states interact with the rest of the SM particles con- 
fined to the 3-brane. We look for possible signals for this interaction in the present 
high-energy collider data, and estimate the sensitivity offered by the next generation 
of collider experiments. Based on the present data from the LEP 2, Tevatron, and 
HERA experiments, we set a lower limit on the extra dimension compactification 
scale Mc > 6.8 TeV at the 95% confidence level (dominated by the LEP 2 results) 
and quote expected sensitivities in the Tevatron Run 2 and at the LHC. 



This contribution is a shortened version of the recent paper [214], with the focus on future high-energy 
facilities. The details of the formalism used to obtain the results presented here can be found in [214|. 

Recently, it has been suggested that the Planck, string, and grand unification scales can all be significantly 
lower than it was previously thought, perhaps as low as a few TeV [20, 215| -217]. An interesting model was 
proposed [218-221], in which matter resides on a p-brane {p > 3), with chiral fermions confined to the 
ordinary three-dimensional world internal to the p-brane and the SM gauge bosons also propagating in the 
extra 6 > dimensions internal to the p-brane. (Gravity in the bulk is not of direct concern in this model.) 
It was shown [218] that in this scenario it is possible to achieve the gauge coupling unification at a scale 
much lower than the usual GUT scale, due to a much faster power-law running of the couplings at the scales 
above the compactification scale of the extra dimensions. The SM gauge bosons that propagate in the extra 
dimensions compactified on S^/Z2, in the four-dimensional point of view, are equivalent to towers of Kaluza- 
Klein (KK) states with masses M„ = ^/Mf+n^]W (n = 1,2,...), where R = M^^ is the size of the 
compact dimension, Mc is the corresponding compactification scale, and Mq is the mass of the corresponding 
SM gauge boson. 

There are two important consequences of the existence of the KK states of the gauge bosons in collider 
phenomenology, (i) Since the entire tower of KK states have the same quantum numbers as their zeroth- 
state gauge boson, this gives rise to mixings among the zeroth (the SM gauge boson) and the nth-modes 
(n = 1, 2, 3, ...) of the W and Z bosons. (The zero mass of the photon is protected by the U(1)em symmetry 
of the SM.) (ii) In addition to direct production and virtual exchanges of the zeroth-state gauge bosons, both 
direct production and virtual effects of the KK states of the W, Z, 7, and g bosons would become possible at 
high energies. 

In this proceedings, we study the effects of virtual exchanges of the KK states of the W, Z, 7, and g 
bosons in high energy collider processes. While the effects on the low-energy precision measurements have 
been studied in detail 1 222| - 229| ], their high-energy counterparts have not been systematically studied yet. We 
attempt to bridge this gap by analyzing all the available high-energy collider data including the dilepton, dijet, 
and top-pair production at the Tevatron; neutral and charged-current deep-inelastic scattering at HERA; and the 
precision observables in leptonic and hadronic production at LEP 2. 
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We fit the observables in the above processes to the sum of the SM prediction and the contribution from 
the KK states of the SM gauge bosons. In all cases, the data do not require the presence of the KK excitations, 
which is then translated to the limits on the compactification scale Mc- The fit to the combined data set 
yields a 95% C.L. lower limit on Mc of 6.8 TeV, which is substantially higher than that obtained using only 
electroweak precision measurements. In addition, we also estimate the expected reach on Mc in Run 2 of the 
Fermilab Tevatron and at the LHC, using dilepton production. 



Interactions of the Kaluza-Klein States 



We use the formalism of Ref. [ ]219| - |221[ ], based on an extension of the SM to five dimensions, with the fifth 
dimension, x^, compactified on the segment /Z2 (a circle of radius R with the identification —x^). 
This segment has the length of -kR. Two 3-branes reside at the fixed points = and x^ = ttR. The SM 
gauge boson fields propagate in the 5D-bulk, while the SM fermions are confined to the 3-brane located at 
x^ = 0. The Higgs sector consists of two Higgs doublets, (pi and 02 (with the ratio of vacuum expectation 
values V2/V1 = tan /?), which live in the bulk and on the SM brane, respectively. 

In the case of SU(2)l x U(l)y symmetry, the charged-current (CC) and neutral-current (NC) interactions, 
after compactifying the fifth dimension, are given by []225[] : 



C 



CC 



+ 



^NC ^ ^ 



+ 



Wf + cos^ PY^W^'^^f + 2\/2 sin^ pWi ^ wi""^ + 2 sin^ /3 

n=l n=l 

I Y n'MUw[-^f - g{W^ + ^2 W^f )'^) + (1 - 2) , 

n=l 

00 00 / 00 

+ cos2 /3 + 2v^sin2 /3Z ^ Z^") + 2 sin^ (3 ( ^ Z^") 

n=l 




(1) 



gv 
8? 



n=l 
2 



n=l 



\n=l 



n=l 



S0C0 



Z^^^Y ^^"^^ J -^/f - e ( + V2 ^ 



rem 



(2) 



ra=l 



n=l 



where the fermion currents are: 



71,2 



ri,2 
2 



■01 



J,1 



rem 



and 



= V cos /?, {(P2) = v sin (3; g and g' are the gauge coupUngs of the SU(2)i and U(l)y , respectively; 
gy = T3L/2 — SgQ and ga = T^l/'^- Here, we used the following short-hand notations: sq = s'mOw 
and C0 = cos 6w, where 9w is the weak-mixing angle. The tree-level (non-physical) W and Z masses are 
Mw = gv/2 and Mz = Mw/cg. Since the compactification scale Mc is expected to be in the TeV range, we 
therefore ignore in the above equations the mass of the zeroth-state gauge boson in the expression for the mass 

of the n-th KK excitation: M„ = y'M?+^^^^^^ ^ nMc, n = 1,2, .... 

Using the above Lagrangians we can describe the two major effects of the KK states: mixing with the 
SM gauge bosons and virtual exchanges in high-energy interactions. 
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1.1 Mixing with the SM Gauge Bosons 

The first few terms in the Eqs. (|l|) and (Q) imply the existence of mixings among the SM boson (V) and its 
KK excitations {V'-'^\ V^^\ ...) where V = W,Z. There is no mixing for the fields because of the U(1)em 
symmetry. These mixings modify the electroweak observables (similar to the mixing between the Z and Z'). 
The SM weak eigenstate of the Z-boson, Z^^\ mixes with its excited KK states Z^") (n = 1, 2, ...) via a series 
of mixing angles, which depend on the masses of Z^'^\n = 0, 1, ... and on the angle /?. The Z boson studied at 
LEP 1 is then the lowest mass eigenstate after mixing. The couplings of the Z^^^ to fermions are also modified 
through the mixing angles. The observables at LEP 1 can place strong constraints on the mixing, and thus on 
the compactification scale Mc- Similarly, the properties of the W boson are also modified. 

The effects of KK excitations in the low-energy limit can be included by eliminating their fields using 
equations of motion. From the Lagrangians given by Eqs. (|l|) and the W, Z masses and the low-energy CC 
and NC interactions are given by [l225p : 



Ml 



M, 



M|(l-si4x) 



C 



CC 
int 



NC 
int 



X 



2M| 



XJU^^ + (1^2) 



sece 



In the following, we summarize the results presented in Refs. [223-229]. Nath and Yamaguchi [223] 



used data on Gp, Mw, and Mz and set the lower limit on Mc ^ 1-6 TeV. Carone [227] studied a number 
of precision observables, such as Gp, p, Qw^ leptonic and hadronic widths of the Z. The most stringent 
constraint on Mc comes from the hadronic width of the Z: Mc > 3.85 TeV. Sti^umia [226] obtained a 
limit Mc > 3.4 — 4.3 TeV from a set of electroweak precision observables. Casalbuoni et al. [225] used 



the complete set of precision measurements, as well as Qw and i?i^'s from i^-N scattering experiments, and 



obtained a limit Mc > 3.6 TeV. Rizzo and Wells [224] used the same set of data as the previous authors and 



obtained a limit Mc > 3.8 TeV. Cornet et al. [229] used the unitarity of the CKM matrix elements and were 
able to obtain a limit Mc > 3.3 TeV. Delgado et al. [228] studied a scenario in which quarks of different 



families are separated in the extra spatial dimension and set the limit Mc > 5 TeV in this scenario. 



1.2 Virtual Exchanges 

If the available energy is higher than the compactification scale the on-shell production of the Kaluza-Klein 
excitations of the gauge bosons can be observed [|l6, 230 1. However, for the present collider energies only 
indirect effects can be seen, as the compactification scale is believed to be at least a few TeV. These indirect 
effects are due to virtual exchange of the KK-states. 

When considering these virtual exchanges, we ignore a slight modification of the coupling constants to 
fermions due to the mixings among the KK states and so we use Eqs. ([I]) and (^) without the mixings.0 This 
implies that any Feynman diagram which has an exchange of a W, Z, 7, or g will be replicated for every 

**For Ale >> Mz the mixings are very small. Furthermore, they completely vanish for f3 — 0. 
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corresponding KK state with the masses nMc, where n = 1,2, .... Note that the couphng constant of the 
KK states to fermions is a factor of \/2 larger than that for the corresponding SM gauge boson, due to the 
normalization of the KK excitations. 

The effects of exchanges of KK states can be easily included by extending reduced amplitudes. In the 
limit Mc ^ \/s, \/\t\, \/\u\, the reduced amplitudes take the form: 



M^i( ^ - , dig} ( 9^9% \ vr^ 



c 



based on which, the high energy processes can be described. 



2. High Energy Processes and Data Sets 

Before describing the data sets used in our analysis, let us first specify certain important aspects of the analysis 
technique. Since the next-to-leading order (NLO) calculations do not exist for the new interactions yet, we use 
leading order (LO) calculations for contributions both from the SM and from new interactions, for consistency. 
However, in many cases, e.g. in the analysis of precision electroweak parameters, it is important to use the 
best available calculations of their SM values, as in many cases data is sensitive to the next-to-leading and 
sometimes even to higher-order corrections. Therefore, we normalize our leading order calculations to either 
the best calculations available, or to the low-Q^ region of the data set, where the contribution from the KK 
states is expected to be vanishing. This is equivalent to introducing a -dependent K-factor and using the 
same ii'-factor for both the SM contribution and the effects of the KK resonances, which is well justified by 
the similarity between these extra resonances and the corresponding ground-state gauge boson. The details of 
this procedure for each data set are given in the corresponding section. Wherever parton distribution functions 



(PDFs) are needed, we use the CTEQ5L (leading order fit) set | |231| ] 
2.1 HERA Neutral and Charged Current Data 



ZEUS [ |232| , [233[ ] and HI [ [234{ , |235[ ] have published results on neutral-current (NC) and charged-current (CC) 



deep-inelastic scattering (DIS) in e^-p collisions at « 300 GeV. The data sets collected by HI and ZEUS 



correspond to integrated luminosities of 35.6 and 47.7 pb^^, respectively. HI [234, 235] has also published a 
NC and CC analysis for the most recent data collected in e^p collisions at k, 320 GeV with an integrated 
luminosity of 16.4 pb~^. We used single-differential cross sections da/dQ'^ presented by ZEUS [232, 233 ] 
and double-differential cross sections d'^a /dxdQ^ published by HI [ ]234| , |235 |. 



We normalize the tree-level SM cross section to that measured in the low-Q^ data by a scale factor C (C 
is very close to 1 numerically). The cross section a used in the fitting procedure is given by 

0" = C ((TSM + CTintcrf + O-Rk) , (3) 

where o"intcrf is the interference term between the SM and the KK states and ukk is the cross section due to the 
KK-state interactions only. 

2.2 Drell-Yan Production at the Tevatron 



Both CDF [ ]236| ] and Dj9 [ )237| ] measured the differential cross section da/dMu for Drell-Yan production, 
where Mu is the invariant mass of the lepton pair. (CDF analyzed data in both the electron and muon channels; 
Dp analyzed only the electron channel.) 

We scale this tree-level SM cross section by normalizing it to the Z-peak cross section measured with 
the data. The cross section used in the fitting procedure is then obtained similarly to that in Eq. (|3|). 



87 



2.3 LEP 2 Data 

We analyze LEP 2 observables sensitive to the effects of the KK states of the photon and Z, including hadronic 
and leptonic cross sections and forward-backward asymmetries. The LEP Electroweak Working Group com- 
bined the gq, /i^/i", and r+r" data from all four LEP collaborations [238] for the machine energies between 
130 and 202 GeV. We use the following quantities in our analysis: (i) total hadronic cross sections; (ii) to- 
tal /i+/x^,T+r~ cross sections; (iii) forward-backward asymmetries in the /x and r channels; and (iv) ratio 
of 6-quark and c-quark production to the total hadronic cross section, and R^. We take into account the 



correlations of the data points in each data set as given by [238] 



For other channels we use various data sets from individual experiments. They are [23S-254]: (i) 
Bhabha scattering cross section cr(e^e^ e~^e~); (ii) angular distribution or forward-backward asymmetry 
in hadroproduction e~^e~ —>■ qq; (iii) angular distribution or forward-backward asymmetry in the e~^e~, n^fJ.^, 
and r^r^ production. 

To minimize the uncertainties from higher-order corrections, we normalize the tree-level SM calculations 
to the NLO cross section, quoted in the corresponding experimental papers. We then scale our tree-level results, 
including contributions from the KK states of the Z and 7, with this normalization factor, similar to Eq. (^). 
When fitting angular distribution, we fit to the shape only, and treat the normalization as a free parameter. 



2.4 Kaluza-Klein states of the Gluon in Dijet Production at the Tevatron 

Since the gauge bosons propagate in extra dimensions, the Kaluza-Klein momentum conservation applies at 
their self-coupling vertices. Because of this conservation, the triple interaction vertex with two gluons on the 
SM 3-brane and one KK state of the gluon in the bulk vanishes. (However, the quartic vertex with two gluons 
on the SM 3-brane and two gluon KK states in the bulk does exist.) The cross sections for dijet production. 



including the contributions from KK states of the gluon, are given in Ref. [ ]214| ]. 



Both CDF [255, 256] and Dp [257, 258| published data on dijet production, including invariant mass 



Mjj and angular distributions. In the fit, we take into account the full correlation of data points in the data sets, 
as given by each experiment. 



2.5 Kaluza-Klein States of the Gluon in tt Production at the Tevatron 



In Ref. [ 259 1 , it was shown that tt production in Run 2 of the Tevatron can be used to probe the compactification 
scales up to ~ 3 TeV. In this paper, we consider the sensitivity from the existing Run 1 data by using the tree- 
level tt production cross section, including the contribution of the KK states of the gluon in the qq tt 
channel. (The gg — > tt channel does not have the triple vertex interaction with two gluons from the SM 3-brane 
and one KK state of the gluon in the bulk, as explained in the previous subsection.) 

The latest theoretical calculations of the ti cross section, including higher-order contributions, at ^/s = 
1.8 TeV correspond to 4.7 - 5.5 pb [|260|, 261 1. The present data on the tt cross sections are [262, 263] 



crtt (CDF) 



: 6.5 til pb; 



a^f (D0) = 5.9±L7 pb, 



and the top-quark mass measurements are 

mt (CDF) = 176.1 ± 6.6 GeV; mt (D^) = 172.1 ± 7.1 GeV. 

In our analysis, we normalize the tree-level SM cross section to the mean of the latest theoretical predictions 
(5.1 pb), and use this normalization coefficient to predict the cross section in the presence of the KK states of 
the gluon (similar to Eq. (||)). 
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3. Constraints from High Energy Experiments 

Based on the above individual and combined data sets, we perform a fit to the sum of the SM prediction and 
the contribution of the KK states of gauge bosons, normalizing our tree-level cross section to the best available 
higher-order calculations, as explained above. The effects of the KK states always enter the equations in the 
form rj = 7r^/(3M^) [ |214| ]. Therefore, we parameterize these effects with a single fit parameter r]. In most 
cases, the differential cross sections in the presence of the KK states of gauge bosons are bilinear in rj. 

The best-fit values of r/ for each individual data set and their combinations are shown in Table |l|. In all 
cases, the preferred values from the fit are consistent with zero, and therefore we proceed with setting Umits on 
r]. The one-sided 95% C.L. upper limit on t] is defined as: 

dr] P{ri) 

95 = — 111 (4) 

where P{ri) is the fit likehhood function given by P{ri) = exp(— (x^(r/) — Xmin)/^)- The corresponding upper 
95% C.L. limits on i] and lower 95% C.L. limits on Mc are also shown in Table |l|. 



Table 1: Best-fit values of 77 = tt /(3M^)andthe95%C.L. upper limits on -q for individual data set and combinations. Corresponding 
95% C.L. lower limits on Mc are also shown. 





7] (TeV-2) 


r/95 (TeV-2) 


M^s (TeV) 


LEP2: 








hadronic cross section, ang. dist., ii^ c 


-0.33 


0.12 


5.3 


/i, r cross section & ang. dist. 


0.09 


0.42 


2.8 


ee cross section & ang. dist. 


-0.62 tlf. 


0.16 


4.5 


LEP combined 


n 90 +0.092 
— U.ZS _o.092 


0.076 


6.6 


HERA: 








NC 


-2.74 t\f. 


1.59 


1.4 


cc 


-0.057 


2.45 


1.2 


HERA combined 


, r,o +0.98 
— i.Zd _o.99 


1.25 


1.6 


TEVATRON: 








Drell-Yan 


-o.87ll:Si 


1.96 


1.3 


Tevatron dijet 


n +0-37 
0.46 Io.58 


1.0 


1.8 


Tevatron top production 


— U.OCi _o.49 


9.2 


0.60 


Tevatron combined 


-0.38 


0.65 


2.3 


All combined 


n on +0.090 

— u.zy -0.090 


0.071 


6.8 



4. Sensitivity in Run 2 of the Tevatron and at the LHC 

At the Tevatron, the best channel to probe the KK states of the photon or the Z boson is Drell-Yan production. 
In Ref. |264|, we showed that using the double differential distribution d^a /MudcosO can increase the 
sensitivity to the KK states of the graviton compared to the use of single-differential distributions. Similarly, 
we expect this to be the case for the KK states of the photon and the Z boson. 



We follow the prescription of Ref. [264] and use the Bayesian approach, which correctly takes into 
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account both the statistical and systematic uncertainties, in the estimation of the sensitivity tor] = vr^/ (3Af|H,).[^ 
Due to the high statistics in Run 2 and particularly at the LHC, the overall systematics becomes dominated by 
the systematics on the s-dependence of the K-factor from the NLO corrections. (Systematic uncertainties on 
the integrated luminosity and efficiencies are not as important as before, because they get canceled out when 
normalizing the tree level SM cross section to the Z-peak region in the data.) The uncertainty on the i^-factor 
from the NLO calculations for Drell-Yan production [ |265[ ] is currently known to a 3% level, so we use this as 
the correlated systematics in our calculations on Mc- For the LHC we quote the limits for the same nominal 
3% uncertainty and also show how the sensitivity improves if the uncertainty on the ii'-f actor shape is reduced 
to a 1% level. It shows the importance of higher-order calculations of the Drell-Yan cross section, which we 
hope will become available by the time the LHC turns on. 

In the simulation, we use a dilepton efficiency of 90%, a rapidity coverage of \7]\ < 2.0, and typical 
energy resolutions of the Tevatron or LHC experiments. The simulation is done for a single collider experiment 
in the combination of the dielectron and dimuon channels. 

As expected, the fit to double-differential cross sections yields a ~ 10% better sensitivity to Mc than 
just using one-dimensional differential cross sections. We illustrate this by calculating the sensitivity to Mc in 
Run 1, which is slightly higher than the result obtained from the fit to the invariant mass spectrum from CDF 
and Dp. The sensitivity, at the 95% C.L., to Mc in Run 1 (120 pb^i). Run 2a (2 fb^i). Run 2b (15 fb^i), 
and at the LHC (100 fb-i) is given in Table ^ While the Run 2 sensitivity is somewhat inferior to the current 
indirect limits from precision electroweak data, LHC would offer a significantly higher sensitivity to Mc, well 
above 10 TeV. 

Table 2: Sensitivity to the parameter rj = tt^/SMc in Run 1, Run 2 of the Tevatron and at the LHC, using the dilepton channel. The 
corresponding 95% C.L. lower limits on Mc are also shown. 





7?95 (TeV-2) 


95% C.L. lower hmit on Mc (TeV) 


Run 1 (120 pb" ^) 


1.62 


1.4 


Run 2a (2 fb^^) 


0.40 


2.9 


Run 2b (15 fb"^) 


0.19 


4.2 


LHC (14 TeV, 100 fb~\ 3% systematics) 


1.81 X 10"^ 


13.5 


LHC (14 TeV, 100 fb-\ 1% systematics) 


1.37 X IQ-^ 


15.5 
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Part XIV 

Kaluza-Klein Excitations of Gauge Bosons in the 
ATLAS Detector 

G. Azuelos and G. Polesello 

Abstract 

Kaluza-Klein excitations of the gauge bosons are a notable feature of theories with 
"small" (~ 1 TeV) extra dimensions. The leptonic decays of the excitations of 7 and 
Z bosons provide a striking signature which can be detected at the LHC. We investi- 
gate the reach for these signatures through a parametrized simulation of the ATLAS 
detector. With an integrated luminosity of 100 fb~^ a peak in the lepton-lepton invari- 
ant mass will be detected if the compactification scale (Afc) is below 5.8 TeV. If no 
peak is observed, with an integrated luminosity of 300 fb^^ a limit of Mc < 13.5 TeV 
can be obtained from a detailed study of the shape of the lepton-lepton invariant mass 
distribution. If a peak is observed, the study of the angular distribution of the two lep- 
tons will allow to distinguish the KK excitations from alternative models yielding the 
same signature. 



1. Introduction 



In models with "large" Extra Dimensions, characterized by compactification radii » 1 /TeV , gravity propa- 
gates in the bulk, and the SM fields are confined to a 3-brane. The presence of the Extra Dimensions could be 
probed by searching for the Kaluza-Klein excitations of the gravitons at the future high energy accelerators, 
and these scenarios have been the subject of many phenomenological studies. An interesting variation of the 



ADD model [218-221] assumes that only the fermions are confined in the 3-brane, whereas the gauge fields 
propagate in a number of additional "small" extra dimensions orthogonal to the brane with compactification 
radius ~ 1 TeV. 

For definiteness we concentrate here on a model with only one "small" extra dimension, and where all 
of the SM fermions are on the same orbifold point {D = 0). The phenomenology of this model, which we 
will label as MI is discussed in some detail in [23C]. For compactification on jZ'^ dimension, the model 
is completely specified by a single parameter M^, the compactification scale, and the masses Af„ of the KK 
modes of the gauge bosons are given by the relation = {nM^)^ + Mq, where Mq is the mass of the 
zero-mode excitation corresponding to the Standard Model gauge boson. The couplings are the same as the 
corresponding SM couplings, scaled by a factor \pi. As an example of variation on our reference model we 



also briefly consider an alternative model, | |21| ] (M2), where quark and leptons are at opposite fixed points. The 
difference between the reference model and M2 is that for model M2 the signs of the quark couplings of the 
bosons are reversed for excitations with n odd, yielding a different interference pattern with the SM Drell-Yan 
production. 

The constraints on the compactification scale from precision electroweak measurements have been eval- 



uated in a number of papers, [ ]223[ - [225[ , [227[ - |230| , |266| ]. These studies estimate an approximate lower limit 
of 4 TeV on the compactification scale for the reference model considered in this analysis. A recent paper 



[ |214[ , |267| ] calculates the limits which can be extracted from precision measurements at present high-energy 
accelerators. A 95% limit of 6.8 TeV is obtained, dominated by the LEP 2 measurements. The limit, however, 
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comes from the fact that for two of the three fits to LEP data, an unphysical negative value for Mc is obtained, 
with a significance of two to three standard deviations. In view of this fact, waiting for a clarification of the 
claimed discrepancy with the Standard Model, we conservatively study the performance of the ATLAS detector 
starting from = 4 TeV. 



2. Signal simulation and data analysis 

We simulate at particle level the production of the gauge boson excitations, including the full interference and 
angular information. We include the full Breit-Wigner shape for the first two excitations of 7 and Z [268], 
and a resummed expression for the higher lying states, for which the approximation M{i) ^ s is used. The 
expression with only the first two resonances, does not alter the results for the reach in the peak region, but it 
significantly underestimates the deviation from the Standard Model in the low mass off-peak region. Since the 
dominant contributions to the low s off -resonance region comes from the interference term between SM 7/Z 
and the KK excitations, the deviation from the SM is approximately proportional to: 



1 °° 1 

— T — 

^ 71=1 



(1) 



Therefore the deviation from the SM increases by ^ - 1 ~ 30% when the full tower of resonances is 
considered instead of just the first two. If we consider model M2, the sum over the tower of resonances gives 
a term proportional to 

00 00 
_|_ v^/_-|\nJ_ ^ _i:_L\^_ 

^ n=l ^ n=l 

Therefore, the summed contribution of the interference terms in model M2 will be of opposite sign and half of 
the one for the reference model. The matrix elements are interfaced to PYTHIA 6.125 [ ]113| ] event generator 
as an external process, and full events have been generated, including the full PYTHIA machinery for QCD 
showering from the initial state quarks, and for the hadronization. 



The events thus generated have been passed through the fast simulation of the ATLAS detector [269]. 
As discussed in the introduction, the lowest Mc considered in this study is 4 TeV, consistent with precision 
electroweak measurements. Therefore we need to detect and measure leptons with momenta in the few TeV 
range. In this range the energy resolution for electrons is dominated by a constant term due to the imperfect 
knowledge of the detector performance. From studies performed on test beam data and on fully simulated 
events, for energies up to a few hundred GeV, this term has been evaluated as a few per mill. Detailed studies 
need to be performed to evaluate how well these results extrapolate to the momentum range of interest for 
this analysis. With this caveat, we use here the standard parametrization included in the ATLFAST program 
which yields a resolution of ~0.7% for the energy measurement of 2 TeV electrons. For muons the ti^ansverse 
momentum measurement of high Pt muons is achieved through the sagitta measurement in the precision drift 
chambers, and for a 2 TeV muon the resolution is of order 20%. Considering the natural widths of the gauge 
excitations, the width of the lepton-lepton invariant mass distribution will be dominated by the natural width for 
electrons, and by the experimental momentum resolution for muons. This is illustrated in Figure [l| where the 
invariant mass spectra for a 4 TeV KK resonance is shown both for electrons (full line) and muons (dashed line). 
Although the muon peak is much broader, both lepton species can be used in order to observe the existence of 
an excess in the peak region with respect to the Standard Model. 

Comparing the two-lepton invariant mass spectrum for Standard Drell-Yan production (full line), and for 
the reference model (dashed line) as shown in Fig.^, two basic features can be observed: 

• A peak centered around M^, corresponding to the superposition of the 7(1) and Z{1) Breit-Wigner 
shapes 
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Fig. 1: Distribution of the lepton lepton invariant mass for electrons (full line) and muons (dashed line). The distribution assumes 4 TeV 
for the mass the lowest lying KK excitation. 



• A suppression of the cross-section with respect to the SM for masses below the resonance. This sup- 
pression is due to the negative interference terms between the SM gauge bosons and the whole tower 
of KK excitations, and is sizable even for compactification masses well above the ones accessible to a 
direct detection of the mass peak. This shape is the consequence of the model choices requiring both the 
leptons and the quarks to be at the same orbifold point (D=0). The different model choices corresponding 
to M2 would yield an enhancement of the off peak cross-section, as shown in the dotted line in Fig. ^. 
We select events with two isolated opposite sign leptons, satisfying the following requirements: 

• rrii+i- > 1000 GeV {£ = e, ft) 

• P^>20 GeV, \r]e\ < 2.5 

The isolation criterion consists in requiring a transverse energy deposition in the calorimeter smaller than 
10 GeV in a (rj, 4>) cone of radius 0.2 around the lepton direction. In the absence of new physics, approximately 
500 events survive these cuts for an integrated luminosity of 100 fb~^, corresponding to one year of high 
luminosity LHC running for each of the lepton flavors. 

The reach for the observation of a peak in the m^+i- distribution can be naively estimated from Table [11, 
which, for both electrons and muons gives the number of signal and background events for an integrated lumi- 
nosity of 100 fb^^ for different values of Mc- As an arbitrary requirement for discovery we ask for the detection 
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ATLAS e'^e" 



1 




2000 4000 6000 

m.| (GeV) 



Fig. 2: Invariant mass distribution of e e pairs for the Standard Model (full line) and for models Ml (dashed line) and M2 (dotted 
line). The mass of the lowest lying KK excitation is 4000 GeV. The histograms are normalized to 100 fb~^. 



above a given M« of 10 events summed over the two lepton flavors, and a statistical significance S/ \fB > 5. 
The reach thus calculated is ~5.8 TeV for an integrated luminosity of 100 fb^^. In order to achieve this reach, 
a good control of high ma background events which might be produced by the mismeasurement of leptons is 
crucial. A handle on these events is however provided by the consideration of the momentum balance of the 
event in the transverse plane, which will allow to reject events with one badly mismeasured lepton. 

Unfortunately, even for the lowest allowed value of M^, 4 TeV, no events would be observed for the 
second resonance at 8 TeV, which would have been the most striking signature for this kind of model. 

In order to fully evaluate the sensitivity of the invariant mass spectrum off-resonance to interference 
effects, a likelihood fit to the expected spectrum can be performed, and will be discussed in the next section. 
As a first approach, one can simply evaluate the variation in number of events within a given m^^ range with 
respect to the SM, as a function of M^. We show the invariant e+e~ mass spectrum between 1000 and 2000 
GeV in Figure ^for Standard Model and for three choices of Mc. 

A naive parameterization of the statistical significance of the cross-section suppression is 
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Me(GeV) 


Cut (GeV) 






Background (e) 


Background (/x) 


4000 


3000 


172 


156 


1.45 


1.8 


5000 


4000 


23 


20 


0.15 


0.22 


5500 


4000 


9 


8 


0.15 


0.22 


6000 


4500 


3.3 


2.8 


0.05 


0.1 


7000 


5000 


0.45 


0.38 


0.015 


0.05 


8000 


6000 


0.042 


0.052 


0.0015 


0.012 



Table 1: Expected number of events in the peak for an integrated luminosity of 100 fb~^, for different values of the mass of the lowest 
lying KK excitation, and Standard Model Drell-Yan background. The peak region is defined by requiring a minimum l~ invariant 
mass as shown in the second column. The results for electrons and muons are given separately. 
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Fig. 3: Invariant mass distribution ofe^e~ pairs in the region below 2 TeV. The Standard Model contribution is shown as a thick line. 
We show the reference model with three different values for the compactification scale Mc: 5000, 7000 and 9000 GeV as dashed, dotted 
and dash-dotted lines respectively. The histograms are normalized to 100 fb~^ 



A relevant variable which should also be considered is the ratio Nev{Mc)/Nev{SM), because the systematic 
uncertainty in our knowledge of the shape of m« sets a limit on the detectable value of this ratio. The choice 
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Mc(GeV) 




SM 


483 


4000 


210 


5000 


295 


5500 


324 


6000 


349 


7000 


381 


8000 


405 


8500 


413 


9000 


419 


10000 


432 


11000 


443 


12000 


450 



Table 2: Expected number of events in the interference region for an integrated luminosity of 100 fb ^, for different values of the 
compactification scale Mc and Standard Model e^e~ background (1 lepton flavor). The considered mass interval is 1000 < e^e~ < 
2500GeV 



of the mass interval is subject to the consideration of the systematical uncertainty, as the statistical signifi- 
cance somewhat increases by lowering the lower limit of the considered mass window, at the price of a worse 
Nev{Mc)/Nev{SM). We choose for this analysis a mass interval 1000 < m« < 2500 GeV. 

From the numbers in Table if we consider both lepton flavors, the ATLAS 5 a reach is ~8 TeV for an 
integrated luminosity of 100 fb^^ and ~10.5 TeV for 300 fb^^. The deviation from the Standard Model will 
be 16% for 8 TeV, and ~ 10% for 10.5 TeV, defining in each case the level of systematic control on the relevant 
region of the lepton-lepton invariant mass spectrum we need to achieve to exploit the statistical power of the 
data. 



3. Optimal reach and mass measurement 

In the previous sections we have evaluated in a naive way by simple event counting the Mc range within 
which LHC will be able to observe a peak generated by the KK gauge excitations, and/or a deviation from 
the Standard Model in the mu distribution off-peak. An optimal estimate of the reach can be obtained by 
performing a likelihood fit to the invariant mass shape expected for different values of Mc- Instead of using just 
the invariant mass, we use the full information contained in the events. Ignoring the transverse momentum of 
the £~^£~ system, the event kinematics for a given event i is defined by the variables x\,X2,cos6''. The values 
of xi, X2 have been evaluated from the four-momenta of the detected electrons, according to the formulas: 



xi - X2, rUf+f- = X1X2S 



For the evaluation of cos 9 we use the CoUins-Soper convention [ |270| ], consisting in the equal sharing of the 
i'^£~ system transverse momentum between the two quarks. A basic problem for the likelihood calculation is 
the fact that, as the LHC is a pp collider, it is not possible to know from which direction the quark in the qq 
hard scattering comes from, so only the absolute value of cos 6 can be measured, but not its sign. Part of this 
information can however be recovered, by using the knowledge of xi and X2 and the fact that in the proton 
the X distribution for valence quarks is harder than for anti-quarks. A detailed discussion of the experimental 
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Mc(GeV) 




RMS 


4000. 


4001.96 


10.91 


5000. 


5003.16 


35.91 


5500. 


5502.19 


77.24 


6000. 


6045.22 


216.61 


7000. 


7129.48 


544.35 



Table 3: Average estimated value (Muk) and RMS of Mcfor ~ 2500 experiments and an integrated luminosity of of 100 fb 



reconstruction of the three variables is given in [ 271 1. 

For the processes under study, the initial state is qq for both signal and background, so the optimal 
result can be obtained by just using the two physical variables sensitive to the dynamics of the hard-scattering 
processes, namely mu = ^/x\X2S and cos 6. However, since for electrons the effect of the experimental 
smearing is small, an effective approach to the problem is to use the theoretical cross-section expression to 
build an unbinned likelihood. In this approach, the use of only two variables would require an integration 
over the third one for each step in the likelihood calculation for each Monte Carlo experiment, making the 
process unacceptably slow. For muons the experimental smearing must be taken into account, and the fit can 
be performed by building an event density grid in the mu — cos 9 plane. 

In the following we will only perform the likelihood fit for electrons, calculating the unbinned likeli- 
hood functions on event samples corresponding to an integrated luminosity of 100 fb^^. In order to evaluate 
the uncertainty on the measurement, for each input Mc value we generated an ensemble of Monte Carlo 
experiments (100 fb^^ each) and for each of them we estimated by maximizing the likelihood function. 

The likelihood fit is performed on the variable 1 /M^, since for mu <^ it is the natural variable for describ- 
ing the deviation of the cross-section from the Standaid Model, as shown in Eq. |I} With the use of this test 
variable, the Standard Model is the limit corresponding to 1/M^ = 0, and it is possible to build a continuous 
likelihood function extending the evaluation to unphysical negative values of 1/M^. 

We show in Figure ^ the distributions of the estimated values of 1/M^ for four input values of Mc- As 
expected, the distributions are gaussian as long as events in the peak exist, and tails start to appear for Mc=6 TeV 
for which, on average, only three events appear in the peak region for the considered statistics. For Mc=l TeV, 
less than 1 event is observed in the peak and the distribution becomes very broad, with an RMS corresponding 
to ~600 GeV, and large tails. The average and RMS of the estimated value of Mc are given in Table |[ The 
statistical error is below the percent level as long as events are observed in the peak region. A small systematic 
shift in the average of the estimated Mc is observed, due to the fact that the likelihood is built using analytical 
test functions neglecting detector smearing and transverse motion of the e+e^ system. 



The experimental sensitivity is defined in [ ]272| ] as the average upper limit that would be attained by 
an ensemble of experiments with the expected background and no true signal. To evaluate the sensitivity, we 
therefore produced an ensemble of Monte Carlo experiments for which only SM Drell-Yan was generated. For 
the evaluation of the 95% CL limit for each MonteCarlo experiment we use the following prescription. For 
each Monte Carlo experiment we build the likelihood function £ as a function of l/M^ as described above. 
We then define as 95% CL limit the value of Mc such that the integral of C between zero and I/M^ is 95% 
of the integral between zero and infinity. The experimental sensitivities for one lepton flavor thus obtained 
are respectively 9.5, 11 and 12 TeV for integrated luminosities of 100, 200 and 300 fb^^. These values are 
pessimistic, since they do not take into account the systematic deviation from zero of the estimated l/M^ value 
due to the approximate test function used to perform the study. Correcting for the deviation from zero yields 
an improvement of approximately 200 GeV on these numbers. If we assume similar sensitivity for electrons 
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Fig. 4: Distribution of the value ofl/M^ estimated through the maximization of the likelihood function for a set of ^1500 Monte Carlo 
experiments for four Mc values . The input Mc are respectively 4, 5, 6 and 7 TeV. The integrated luminosity is 100 fb~^. 



and muons, the sensitivity is ^13.5 TeV for 300 fb' and both lepton species. These figures only express 
the statistical sensitivity of the ATLAS experiment, the possible sources of systematic uncertainty must be 
considered in order to evaluate the final ATLAS sensitivity. 
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4. Systematic uncertainties 

As shown in the the previous section, the effect of KK excitation can be detected for well above the mass 
range for which the direct observation of a peak is possible from a detailed study of the event shape in the 
interference region. The experimental sensitivity in this region crucially depends on our understanding of the 
kinematic distributions of the lepton-lepton system both under the experimental and the theoretical point of 
view. 

As shown in Figure ^, as Mc increases, the difference in shape with respect to the standard model 
becomes less and less significant, and systematic effects in the lepton-lepton invariant mass measurement might 
affect the shape of the distribution, and destroy the experimental sensitivity. We consider for this analysis 
electrons of very high momentum, around 1 TeV. At this energy scale the linearity of the lepton momentum 
measurement, as well as the momentum dependence of the acceptance are difficult to assess using the data. In 
fact very few of the leptons from the decay of high momentum W and Z, which could in principle be used to 
perform the measurement will have high enough momentum. From studies performed for lepton calibration in 
ATLAS, we know that the lepton energy scale will be known to 0.1% at the Z mass. We therefore parametrize 
the deviation from linearity as a logarithmic term which is zero for lepton momentum of 100 GeV, and ±1 or 
±5% for momenta of 2 TeV. We perform the likelihood analysis on all our simulated data samples, modifying 
event by event the reconstructed lepton energy according to the logarithmic formula. For the evaluation of 
Mc between 4 and 6 TeV, the relative deviation from the nominal Mc approximately scales with the deviation 
from linearity for 2 TeV leptons, as shown in Fig.^ for 3 values of Mc- 4, 5 and 5.5 TeV. The variation of the 
sensitivity with the assumed value of the deviation from linearity is shown in Figure ^. As discussed above, 
the systematic uncertainty is reflected in a systematic shift of the average Mc estimate, and an overestimate of 
the lepton calibration is going to yield an optimistic evaluation of the Mc value excluded by the experiment. 
Taking the sensitivity values obtained with a negative deviation from linearity, the sensitivity for 100 fb^^ and 
one lepton species is reduced from 9.5 TeV to 9.3 TeV and 8.75 TeV for 1 and 5% deviation respectively. As 
an approximate rule, the experimental limit should be reduced by ~2% for each percent of uncertainty on the 
energy calibration of 2 TeV leptons. 

An additional uncertainty factor is the theoretical systematics on the likelihood calculation. The likeli- 
hood function is built by weighting real events according to a theoretical cross-section formula. Any discrep- 
ancy between the theoretical formula employed and reahty will induce an uncertainty on the measurement of 
Mc- In the likelihood analysis we are not sensitive to an absolute K-factor, since we do not use the absolute 
normalization, but only to distortions of the kinematic distributions of the lepton-lepton system. Three main 
sources of uncertainty can be identified: 

• QCD higher order corrections; 

• electroweak higher order corrections; 

• the parton distribution function (PDF) for the proton. 

The main effect of QCD higher order corrections is the modification of the Pt distribution of the lepton-lepton 
pair, due to radiation from initial state quarks. This effect is taken into account in a very pessimistic way 
in the study on fully generated events, where the likelihood is built from the leading order 2-to-2 Drell-Yan 
expression, and the events are generated with the full PYTHIA machinery for initial state radiation. Therefore, 
the experimental error quoted in the previous section includes a very pessimistic estimate of this effect. In fact 
in a real experiment a more realistic theoretical modelling will probably be used to build the likelihood. 

Electroweak higher order corrections were recently calculated at NLO |273], and shown to be sizable, 
leading to a reduction of the cross-section which varies with the lepton-lepton invariant mass, and can be as 
large as 35% for pp —>■ and m^+^-. The size of these corrections critically depends on the lepton 

identification and isolation criteria, as a substantial part of the higher order effects yield energetic photons 
produced alongside the leptons. The evaluation of the uncertainties on these corrections is thus a complex 
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Fig. 5: Distribution of the estimated masses (lOOfi-^) in GeV as a function of the allowed deviation from linearity for electrons of 
2 TeV momentum for three values of Mc: 4, 5 and 5.5 TeV (left, center and right plot respectively). 



interplay of experimental and tlieoretical considerations which requires a dedicated study which is outside the 
scope of this analysis. 

The shape of the kinematic distributions of the lepton-lepton system, in particular mu has a strong depen- 
dence on the quark and antiquark PDF's in particular for high values of x. All the events were generated with 
the CTEQ4L PDF's. In order to evaluate the effect of the uncertainty on the structure functions parametrization, 
the likelihood fit was performed on the data set thus generated, using a number of structure function sets. To 
this purpose we have selected the sets providing a leading order parametrization, and which are based on the 
latest available data sets. The distributions of estimated masses are shown in Fig.0 for the eight choices of 
structure function sets used for Mc = 4, 5 and 5.5 TeV. The systematic displacement from the true value is 
between 3 and 4 GeV for 4 TeV, increasing to 10-20 GeV for 5 TeV and 20-40 GeV for 5.5 TeV, and it is well 
below the RMS of the distributions given in Table ^. Another notable effect is that the quality of the likelihood 
fit is worse, giving rise to less Gaussian distributions, and sizable tails start to appear for Mc = 5.5 TeV. The 
experimental reach for 100 fb^^ is shown in Fig. |8|, as a function of the structure function set. In the worst case 
the reach is reduced by ~ 200 GeV with respect to CTEQ4L. 
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Fig. 6: Distribution of the expected sensitivity (100 fb ^ ) in GeV as a function of the allowed deviation from linearity for electrons of 
2 TeV momentum 



5. Spin and Asymmetry Measurement 

If a Kaluza-Klein gauge excitation is discovered, one of the ways of distinguishing the signal from a Z', 
predicted by GUT theories, or from a narrow graviton resonance G* is by the angular distribution of the decay 
products, which should be consistent with the spin 1 nature of the excitation, and by the forward-backward 
asymmetry. By adjusting parameters of the models, the cross sections can be made comparable, but, as shown 
above and in | |274 |, the shape of the mass distribution can provide an additional distinguishing criterion. The 
present study compares these distributions, but does not attempt to distinguish the resonances by the shape of 
their mass distributions, by their relative cross sections, nor by the branching ratios. 



5.1 Cases studied 

We studied the following cases 

a) Z^^^/j^^^: this is the case of gauge excitations, model Ml [230], at mass 4 TeV. The process was 
implemented in PYTHIA 6.201. 

b) Z(i)/7(i)-M2: this case of gauge excitation is with the alternative model M2 [EH], also at 4 TeV. The 
process was implemented in PYTHIA 
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(lower plot). 



c) Z': this is a standard model Z'. The same code as for case a) was used, but the first 7 excitation and 
higher excitations of Z and 7 were removed. 

d) G* : This is the case of a narrow graviton resonance, as was studied by [ 275 ] . The process is implemented 
in PYTHIA. In order to reproduce a resonance of width comparable to the Z^^^ /j^^^ of a) above, the 
dimensionless coupling which enters in the partial widths of the G* (PARP(50) in PYTHIA) was set to 
0.8. The reconstructed width is thus o" ~ 82 GeV. The angular distributions depend on the incoming 
partons. The two processes qq G* ^ £^i~ and gg G* ^ t^l^ were generated and added in 
proportion of their cross section. To their sum was added the Standard Model Drell-Yan background 
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qq Z/7 ^ U. 

The mass distributions normalized to a luminosity of 100 fb^^ are displayed in Figs. ^ for the different 
cases. The cross sections for the different processes are summarized in table H. 



process 


a X BR{Z* e+e") (fb) 




4.05 


Z(i)/7(i)-M2 


11.75 


Z' 


4.65 


qq^G* 


0.20 


gg^G* 


0.13 


qq — > e^e^ 


4.83 



Table 4: Nominal cross sections of the different processes, after a preselection \/i > 1 TeV. 
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Mass distributions, 100 fb-1 
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Fig. 9: Mass distributions of the resonances considered: left: electron channel; right: muon channel; (i) black full line, Z'^^^ /')''^\ (ii) 
red dashed line, ^'^^77*^^' model M2, (Hi) dotted green line, Z' , (iv) dash-dotted blue line, G* + SM Drell-Yan 



5.2 Angular Distributions 

As mentioned above, because the colliding particles at LHC are both protons, the forward-backward asymmetry 
is measured with some ambiguity. Assuming that the resonances are produced by qq fusion, the third component 
of the reconstructed momentum of the dilepton system is taken to be the quark direction, since the quark in the 
proton is expected to have higher energy than an antiquark from the sea. 

Events around the peak of the resonance were selected: 3750 GeV < rriee < 4250 GeV or 3250 GeV 
< < 4750 GeV. For these events, the cosine of the angle of the lepton (e^ or /i^), with respect to the 



beam direction, in the frame of the decaying resonance, is shown in Figs. 10 and |ri|. The positive direction is 
defined by the sign of the reconstructed momentum of the dilepton system. Since we will be interested only in 
the shape, and not in the cross sections, the angular distribution histograms have been normalized, to a total of 
138 events, corresponding to the number of events predicted with an integrated luminosity of 100 fb~^ for the 
reference case 

To compare the shape of these distributions, a set of 1000 angular distributions from the different types 
of resonances was generated by sampling from the expected distributions of Figs. [TO and 1 1 . A Kolmogorov 



test was then applied]^ between the expected Z^^"^ /j^^'' distribution and distributions sampled from the other 
resonances. The output of the test is expected to be a uniform distribution between and 1 if they come from 



the same parent distribution. The histogram of the outputs is displayed in Fig. [12|. No significant difference is 
found between models Ml and M2 of Z^^^ /^^^\ as expected. However, the Kolmogorov test, applied to the 
distributions obtained for the e+e" channel, will give an average probability of consistency between Z(^)/7(^) 
and Z' or between Z^^^ /l*-^^ and G* of 0.105 and 0.015 respectively and will reject, at 95% confidence level, 
the hypothesis that the distributions derive from the same parent distribution 52% and 94% of the times. For 
higher resonance masses the statistical significance quickly decreases: at 5 TeV, with only 18 events in the peak 
region, no discrimination becomes possible. However, for this mass but with an integrated luminosity of 300 

**In principle, the Kolmogorov test should be applied on unbinned data, but the application on binned data should still provide a 
valid test, in the present case since the bins are narrower than the expected features 
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Angular distributions, electrons, 138 events 
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Fig. 10: Electron channel: angular distributions for the different types of resonances considered. The distributions are normalized to 
a total of 138 events in the peak. 
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Angular distributions, muons, 139 events 
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Fig. 11: Muon channel: angular distributions for the different types of resonances considered. The distributions are normalized to a 
total of 138 events in the peak. 
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fb , the Kolmogorov test would reject the hypothesis, at 95% CL, about 20% of the times. Similar results are 
obtained for the /i^/x^ channel. 

A test was also performed between these distributions, leading to the same conclusions. Here, also, a 
histogram of the calculated was produced from a sample of 1000 pseudo experiments with 138 events each. 
The average x^/d.f. are 0.998, 1.50 and 2.10 for the cases of model M2, Z' and G* respectively. The goodness 
of fit test between the Z^^^ and the Z' or G* angular distributions would yield a confidence level below 
5% respectively 38% and 84% of the times. 



5.3 Forward-backward asymmetry 

as a 



From the angular distributions, the forward-backward asymmetry is obtained and shown in Figs. 13 and 14 ; 
function of the reconstructed dilepton mass. It allows a clear distinction between a resonance due to Z^^j^^^l 
and either a or a G* resonance. Indeed, the asymmetry is expected to be close to at the mass peak of the 
Z', if the couphngs are those of the SM, because sin^ 9w ~ 1/4: 

A% = ^AgA, (2) 

with (3) 

Ivgai _ 2(1 - 41(5^1 sin^6'i4') 
vj + aj ~ l + (l-4|Q£|sin2 Owf 



= -T^ = ^-. 1 ^ , ■ .; \o -0 (4) 



(5) 



For masses below, but close to the resonance, the FB asymmetry can also serve as a distinguishing 
criterion between the Z' and the Z^^^ /'y^^\ For large masses, however, the discrimination power becomes 
quickly limited by statistics. 
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Fig. 13: Electron channel: measured forward-backward asymmetry at LHC, for different types of resonances, centered at m = 4 TeV. 
The error bars are representative of a sample having 138 events in the peak region, or 100 fb^^ for Z'^^'^ /')^^\ 



6. Conclusions 

We have performed a detailed study of the leptonic signatures for the production of the KK excitations of the 7 
and Z in models with TeV-scale extra-dimensions. 

The production and decay of the excitations were fully simulated, including initial state QCD radiation, 
and the resulting particles were passed through a parametrized simulation of the ATLAS detector. 

We found that with an integrated luminosity of 100 fh^^ ATLAS will be able to detect a peak in the 
lepton-lepton invariant mass if the compactification scale (Mc) is below 5.8 TeV. Even in the absence of a peak, 
a detailed study of the shape of the lepton-lepton invariant mass will allow to observe a deviation from the 
standard model due to the interference of the KK excitations with the SM bosons. From a study based on a 
maximum likelihood estimation of the compactification mass, ATLAS will be able to exclude at 95% CL a 
signal from the models considered in this work for Mc < 13.5 TeV with an integrated luminosity of 300 fh~^. 
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Fig. 14: Muon channel: measured forward-backward asymmetry at LHC, for different types of resonances, centered at 4 TeV. The 
error bars are representative of a sample having 138 events in the peak region, or 100 fb^^ for Z^^"' /j^^K 



We have performed an evaluation of the influence of experimental and theoretical uncertainties on this result. 
A 1% deviation from linearity in lepton momentum measurement yields a 2% reduction in sensitivity. The 
maximum effect observed from the consideration of various sets of PDF's is a reduction of order 200 GeV on 
the achievable limit. 

Once a peak is observed, an important question is the assessment of the model which has produced the 
signal. We show that for resonances of mass up to ^ 5 TeV, and with an integrated luminosity of 300 fb~^, the 
KK excitations can be distinguished from mass peaks produced by SM-like Z' or graviton resonances from the 
study of the polar angle distribution of the leptons in the peak region. The forward-backward lepton asymmetry 
as a function of invariant mass can provide a more general distinguishing criterion among the different models. 
For invariant masses around the peak, it will allow to distinguish the KK excitations from alternative models 
yielding the same signature. 
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Part XV 

Search for the Randall Sundrum Radion Using the 
ATLAS Detector 

G. Azuelos, D. Cavalli, H. Przysiezniak, L. Vacavant 

Abstract 

The possibility of observing the radion (i;^) using the ATLAS detector at the LHC is 
investigated. This scalar, postulated by Goldberger and Wise to stabilize brane fluctu- 
ations in the Randall-Sundrum model of extra dimensions, has Higgs-like couplings. 
Studies on searches for the Standard Model Higgs with the ATLAS detector are re- 
interpreted to obtain limits on radion decay to 77 and ZZ^*). The observability of 
radion decays into Higgs pairs, which subsequently decay into 77 + bb or rr + bb 
is then estimated. 



1. Introduction 

Theories with extra dimensions have recently received considerable attention. One of the most interesting 
incarnations was formulated by Randall and Sundrum (RS) [Q], who postulate a universe with two 4-d surfaces 
(branes) bounding a slice of 5-d spacetime. The SM fields are assumed to be located on one of the branes (the 
TeV brane), while gravity lives everywhere: on the TeV brane, on the Planck brane and in the bulk. The metric 
is exponentially warped in the fifth dimension, allowing for a natural resolution of the hierarchy problem. 

The theory admits two types of massless excitations: the usual 4-d graviton and a graviscalar. In order 



to stabilize the size of the extra dimension without fine tuning of parameters, Goldberger and Wise (GW) [ ]19| ] 
have proposed a mechanism which requires a massive bulk scalar (p, the radion, expected to be lighter than 
the J=2 Kaluza Klein excitations. The presence of the radion is one of the important phenomenological conse- 
quences of these theories of warped extra dimensions [ 276| -279]. The study of this scalar therefore constitutes 
a crucial probe of the model. 

1.1 Radion branching ratios and width 



The radion couplings to fermions and bosons are similar to those of the Standard Model (SM) Higgs [ |276| ]. 
They are expressed as a function of three parameters: the mass of the radion m^, the vacuum expectation value 



of the radion or scale, A^, and ^, the radion-SM Higgs mixing parameter [ |276[ , |280| ]. 

In the following study, we assume that the SM Higgs has been discovered and that its mass has been 
measured. The branching ratios of the radion are calculated using those of the SM Higgs as calculated in 



HDECAY [ |150| ], and using the ratio of the radion to Higgs branching ratios given by [276] 



Figure |I| shows the principal branching ratios as a function of scalar mass for decays of the SM Higgs 
(top plots) and of the radion when m-h = 125GeV/c^ and =1 TeV, for ^ = when there is no (j)-h mixing 
(middle plots), and for ^ = 1/6 when (p and h are heavily mixed (bottom plots). We note the following: 

• BR{(p gg) is greatly enhanced with respect to the Higgs and is close to unity for > 500 GeV / 
and^ = 1/6 

• the radion decays into two SM Higgs for > 2mh 

• BR{(p — > 77) is enhanced for ^ = 1/6 and ~ 600 GeV/c^. 
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Fig. 1: Log(BR) versus the mass of the scalar for the SM Higgs (top), and for the radion when ^ = (middle) and ^ = 1/6 (bottom) 
when A,^ = 1 TeV. The Higgs mass in the lower curves is set to mh = 125GeV/c^. A smaller (larger) mass range is shown on the 
left(right)-hand side. 



For ^ = 1/6 and for a radion with mass close to that of the Higgs, a strong interference produces a strong 
suppression of decays to vector bosons. 
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Fig. 2: Log(r) for the SM Higgs and for the radion, for ^ =0 and 1/6 and for = 1 TeV. 



The radion has a very narrow natural width. Figure § shows the total width as a function of mass, for the 
SM Higgs and for the radion with ^ = and 1/6, for A,^ =1 TeV. The width is inversely proportional to the 
square of A,^. 

The aim of the present study is to investigate the possibility of observing a RS radion with the ATLAS 
detector through the following decays: cf) 77, cf) ZZ^*^ Al, cf) ^ hh. ^ bb77 and — > hh — > bbr+r". 
Only the direct production of the radion gg — > 1?!) is considered. 



2. (/) ^ 77 and ZZ^*) ^ 4£ 
For the 77 (m^ < 160 GeV/c^) and ZZ^*^ (m^ > 100 GeV/c^) decay channels, the radion signal significance 



is determined from the SM Higgs results obtained in the ATLAS TDR [97|, for 100 fb ^ (one year at high 
luminosity 10'^^ cm~^s^^). The ratio of the radion S/^/B over that of the SM Higgs is given by [276]: 



S/VBi^P) r((/<^gg)BR(0^77,ZZ) max(rh 



tot' ' 



S/VB{h) r(h ^ gg)Bi?(h ^ 77, ZZ) Y max(rtt, (J, 



where the mass resolutions are given by am = 0.10-y/m + 0.005m and a"^ = y^{T/2.36)'^ + (0.02m)2. 
Using the ATLAS TDR SM Higgs signal significance results, the radion signal significance is determined and 
shown versus the mass of the radion, in Figure ^, for the 77 channel (top) and for the ZZ^*) channel (bottom), 
for =1,10 TeV, ^=0,1/6, and for an integrated luminosity of 100 fb~^. 



3. (/) — > hh ^ 77bb 

The radion, unlike the SM Higgs but similarly to the ones in the Minimal Supersymmetric Standard Model 
(MSSM), decays into Higgs pairs with relatively high BR (see Figure |l|). As shown in Figure ^, the total width 
of the radion is a factor of 10 (100) smaller for ^ =0 (1/6) than that of the Higgs, such that it is completely 
negligible with respect to the reconstructed mass resolutions. 
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Fig. 3: Signal significance versus tlie mass of tlie radion, for tlie 77 cliannel (top) and for tfie ZZ*-*' ctiannel (bottom). In botfi plots, the 
values for =1,10 TeV and ^ =0,1/6 are shown, for an integrated luminosity of 100 fb~^. 



The specific decay ciiannel ^ hli ^ 77bb offers an interesting signature, with two high-py isolated 
photons and two b-jets. The background rate is expected to be very low for the relevant mass region > 115 
GeV/c^ and > 2mh. In addition, triggering on such events is easy and the diphoton mass provides very 
good kinematical constraints for the reconstruction of m^. 

The decay hh — > 77bb was studied in the context of the MSSM Higgs [ ]281p , although at that time the 
mass ranges investigated were lower. The approach and the selection we use in this study are very similar. 
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3.1 Signal 

Signal events were generated with PYTHIA 6.158 |113|. The modified version of HDECAY described in 
Section |0] is used at the initialization phase to input the correct parameters and at the end of the run to re-scale 
the cross-section. Note that the default PYTHIA parameters, as opposed to the HDECAY ones, are used for the 
light Higgs couplings and parameters. 

The heavy Higgs production process via gluon-gluon fusion (in the framework of the Minimal 1- 
Higgs doublet Standard Model, process ISUB=152)is used to produce the radion. The mass of the was set 
to reflect that of the radion, and the light Higgs mass was set to ruh = 125 GeV/c^. In addition, since the width 
of the radion is much narrower than what is usually encountered in a Higgs scenario, a specific correction was 
implemented [ |282| ] and the branching ratio corrected appropriately. Two samples of 100k events each were 
generated, for = 300 GeV/c^ and for = 600 GeV/c^. 



3.2 Background 

The backgrounds for this channel are 77bb (irreducible), 77CC, 77bj, 77cj and 77jj (reducible with b-tagging). 
The events were generated with PYTHIA 6.158. The main production process is the box diagram gg —>■ 77 
(process ISUB=11 4), where the jets arise from initial state radiation, eventually combined with gluon splitting 
for heavy flavor jets. The rates are therefore very low. However large uncertainties apply to these backgrounds 
since the jet part comes only from radiation and not from the hard-scattering. Generating a background sample 
of a sensible size turns out to be very CPU time consuming, and some cuts had to be applied at the event 
generation: the sample was generated in different bins of (50, 100, 200, 400, 800, 1600 and 3200 GeV/c). 
For each bin, ten million events were generated. 

Single photon production in the hard process 7 j where either the photon or jet is misidentified represents 
another reducible background. This background was studied in the context of the SM H ^ 77 channel, and was 
found to increase the total background by a factor of two. In the context of the radion where the backgrounds 
are negligible, this would not affect the final results. 



3.3 Fast-simulation 



The detector effects on the signal and background events are simulated with ATLFAST 2.53 | |269| ]. While most 
procedures and parameters are the standard ATLFAST ones for low luminosity operation (10^^ cm~^s~^), a 
few improvements are applied for this study: 

• jets are recalibrated using a detailed parameterization 

• the photon reconstruction efficiency is assumed to be 80% 

• a p J- -dependent b-tagging parameterization is used with an average efficiency of £{, = 60% and a rejection 
of approximately 93% for light-quark jets and 7% for c-jets. 



3.4 Selection 

To extract the signal, two isolated photons with pT > 20 GeV/c and \rj\ < 2.5, and two jets with px > 15 
GeV/c, 1 77 1 < 2.5 are required. At least one of the jets has to be tagged as ab. 

The diphoton and the dijet invariant masses are then formed. Figure ^ shows the reconstructed invariant 
masses for rrifj, = 300 GeV/c^, ^ = and = 1 TeV. Subsequently, two mass window cuts are applied by 
requiring that: 

• m-y^ = mji ± 2 GeV/c^ 

• mbj = mh ± 20 GeV/c^. 
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Fig. 4: Diphoton (left) and dijet (right) invariant mass distributions, for = 300 GeV/c^, ^ = 0, = 1 TeV and 30 fb"^ (three 
years at low luminosity 10^^ cm~'^s~^). The right-hand plot shows the impact of requiring two b-tagged jets instead of one. 



The photons and jets fulfilling these requirements are combined to form the m^^bj invariant mass as shown 
in Figure |5[ The mass resolution is improved down to a value of 5 GeV/c^ by constraining the reconstructed 
masses mbj and m^^ to the light Higgs mass mh, as shown on the right-hand plot of Figure ^ The signal 




Fig. 5: Reconstructed 77bj invariant mass distribution, for m0 — 300 GeMIc?, ,f = 0, = 1 TeV and for three years at low 
luminosity. The left-hand plot shows all the combinations and the ones fulfilling the mass window cuts (cf. text). The right-hand 
distribution is obtained by constraining the reconstructed masses mbj and m^.^ to the light Higgs mass mh, after the mass window 
cuts. 
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Cuts 


= 300 GeV/c^ 


7710 = 600 GeV/c^ 


photons kinematics selection 


46% 


51% 


jets kinematics selection 


36% 


28% 


b-tagging 


76% 


78% 


m^^ window cut 


83% 


85% 


mbj window cut 


49% 


53% 


total 


5% 


5% 



Table 1: Acceptance for the signal, for ^ = 0, = 1 TeV and for the two radion masses studied. For each cut the acceptance is 
defined with respect to the previous one. 



X 10 
'Lq 0.8 
8 0.7 
I 0.6 
0.5 



0.4 
0.3 
0.2 
0.1 




80 



100 



120 



140 



160 



m^(GeV/cO 



X 10 




Mean 
RMS 
UDFLW 
OVFLW 



85.01 
51.36 
0.000 
0.2260E-02 



bj 
bb 




50 



100 



150 200 



250 



m^j (GeV/c^) 



Fig. 6: Diphoton (left) and dijet (right) invariant mass distributions for the background sample, for three years at low luminosity. The 
right-hand plot shows the impact of requiring two b-tagged jets instead of one. 

acceptances after the various cuts described above are given in Table 

The same analysis procedure is applied to the background sample. The resulting 771^^ and 77ibj distribu- 
tions are shown in Figure ^ 

Since there are some uncertainties concerning the level of the background, a more conservative approach 
is also tried: the mass window cuts are loosened to keep events fulfilling: 

• m^^ = ni\^ lb 30 GeV/c^ 

• 77ibj = 777h ± 40 GeV/c^ 

The 777^^bj invariant mass distributions for this conservative approach are shown in Figure ^ 

3.5 Results 

The final number of events selected is obtained by counting the candidates after all cuts within a mass window 
of < 77i^^bj > il-ScTm^^bj signal and background. The results are shown in Table ^ Since this channel 
is practically background free, the usual significance S / \fB is not relevant. A signal discovery, defined to be a 
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Fig. 7: Reconstructed 77bj invariant mass distribution for the background around the expected peak for the signal (top: = 300 
GeV/c^; bottom: = 600 GeV/c^), for three years at low luminosity and for the conservative approach (cf. text). The left-hand 
plot shows all the combinations and the ones fulfilling the mass window cuts (cf. text). The right-hand distribution is obtained by 
constraining the reconstructed masses rribj and m-y^ to the light Higgs mass mh, after the mass window cuts. 



minimum of ten events, is straightforward for low values of A,^ early in the physics program of the LHC. This is 
shown in Table ^ where the minimum integrated luminosities needed for discovery are shown. Approximately 
1 fb~^ is needed for A<^ ~ 1 TeV while requiring S > 10 and S/\fB > 5. 

In the special case where ^ = 0, the cross-section is proportional to A^^. Therefore the maximum 
reach in A,^ is derived from this study. This is obtained using the prescription of [ ]272| ] : for a known mean 
background of zero, the signal mean is larger than 10 with 95% CL if the number of observed events is larger 
than 18. The corresponding reach in A^ is 4.6 TeV for = 300 GeV/c^ and 5.7 TeV for = 600 GeV/c^. 
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= 300 GeV/c^ 


= 600 GeV/c 


Z 


background 





background 





background (conserv.) 


1.42 10"'' 


background (conserv.) 





e = 0,A^ = 1 TeV 


380.1 


e = 0,A^ = 1 TeV 


575.8 


^ = 0, = 10 TeV 


3.8 


^ = 0, A0 = 10 TeV 


5.9 


e = l/6,A<^ = l TeV 


680.4 


e = 1/6,A<^ = 1 TeV 


439.9 


^ = 1/6, A,^ = 10 TeV 


5.5 


^ = 1/6, A0 = 10 TeV 


5.9 



Table 2: Number of events selected for background and for signal, for =300 and 600 GeV/c^, for three years at low luminosity 
and for mh — 125 GeV/c^. 





= 300 GeV/c^ 


= 600 GeV/c^ 


e = 0, A^ = 1 TeV 


0.8 


0.5 


^ = 0, A,^ = 10 TeV 


80 


50 


e = 1/6,A^ = ITeV 


0.4 


0.7 


i = 1/6, A^ = 10 TeV 


55 


55 



Table 3: Minimum integrated luminosity (fb ^) needed for discovery. 



4. </) ^ hh ^ bbr+T~ 

The channel (/> ^ hh — > bbr+r" provides another potentially interesting signal for radion discovery, although 
the background is higher and the reconstructed mass resolutions are poorer than in the ^ hh — > 77bb chan- 
nel. 

In order to provide a trigger, a leptonic decay of the r is required. Here, only the case when one r 
decays leptonically and the other hadronically is considered. As above, the events were generated by appro- 
priately adapting the process of MSSM decay of the heavy Higgs H*^ into two light Higgs bosons (h) in Pythia 
6.158 [n3]. The effect of the ATLAS detector on the resolution and efficiency of reconstruction of these 
events was simulated with the ATLAS fast simulation package (ATLFAST 2.53). The efficiency for hadronic r 
reconstruction is assumed to be 40%. For b-jet tagging, an efficiency of 60% is assumed, with a rejection factor 
of 10 for c jets and 100 for light jets 0. 



4.1 Signals and backgrounds 

As in the previous section, the radion mass values investigated are 300 and 600 GeV/c^, while the Hi ggs mass 
is set to 125 GeV/c^. 

The fast simulated samples are: 

• hh ^ bb r+T~ with one r decaying leptonically and the other hadronically (10 000 events) 

• tt ^ bW+ bW^ with one W decaying leptonically and the other hadronically (10^ events) 

• Z ^ r+T^ with one r decaying leptonically and the other hadronically (10^ events). Initial and final 
state radiation provide additional jets which can fake the signal. 

• inclusive W bosons decaying leptonically (2x10^ events). 
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Fig. 8: Reconstructed rr invariant mass for the signal (dots, 
arbitrary normalization) and for the background: tt (full line, 30 
fb~^) and Z ^ rr (dashed line, 30 fb"^). 



Fig. 9: Reconstructed bb invariant mass for the signal (dots, 
arbitrary normalization) and for the background: tt (full line, 30 
fb"^) and Z ^ rr (dashed line, 30 fb"^). 



4.2 The selection 

The study is performed assuming conditions of low luminosity (10'^'^ cm~^s~^) since, at high luminosity 
(10^^ cm^^s^^), the reconstructed tt mass resolution is seriously compromised by pile-up effects. [97]. 
The events are selected if they satisfy the following criteria: 

• A lepton is required with > 25 GeV, |r/^| < 2.5 (this lepton serves as a trigger). 

• The transverse mass — p™*^** is required to be < 40 GeV/c. This cut rejects background events 
containing W bosons. 

• The TT invariant mass is determined by combining the lepton with a tagged T-jet having p!^^ > 15 GeV/c, 
|r/j'^*| < 2.5 (see Figure P). If more than one jet is tagged as a tau-jet, the combination with the mass 
nearest to m-h is chosen. 

• A pair of b-tagged jets with pT > 15 GeV/c and \7]\ < 2.5 is required and their jet-jet mass reconstructed 
(see Figure^. If more than two jets are tagged as b-jets, the pair having the invariant mass closest to the 
Higgs mass is chosen. 

• Cuts on the reconstructed tt mass and bb mass are applied: 

110 < rrirr < 140 GeV/c^ and 90 < m^^ < 140 GeV/c^ in the case of the 300 GeV/c^ radion, and 
110 < rrirr < 150 GeV/c^ and 85 < m^^ < 130 GeV/c^ in the case of the 600 GeV/c^ radion. 



4.3 Results 

Although the signal efficiency is low, the background rejection is high. The expected cross sections for signal 
and background before the event selection are given in Table ^ for A,^ = 10 TeV and ^ = 0. The branching 
ratios account for leptonic decays into a muon or an electron. 

Figures [I^ and ^ show the reconstructed masses for signal when m<^=300 and 600 GeV / respectively, 
for 30 fb~^, = 1 GeV and ^ = 0. The shape for a 300 GeV/c^ radion is similar to the one for background 
(mostly tt), therefore systematic errors will most probably be dominated by the understanding of the level of 
this background. 
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Signal 


m<^ = 300 GeV/c^ 


m^= 600 GeV/c^ 


o-(gg 4>) 


290 fb 


60 fb 


BR((/) hh) 


0.30 


0.25 


BR(hh TT bb) 


2 X 0.06 X 0.8 


2 X 0.06 X 0.8 


BR(rr — > £ + hadrons) 


2 X 2 X 0.17 X 0.65 


2 X 2 X 0.17 X 0.65 


crxBR = 


3.98 fb 


0.652 fb 


tt ^ WbWb 1 + hadrons 


~180 pb 




~40000 pb 


7j ^ TT ^ I + hadrons 


~ 730 pb 



Table 4: Expected cross sections for =10 TeV and ^ = for signal and background before the event selection cuts. {I = e, fj.) 




Fig. 10: Reconstructed mass of the radion (m^ = 300 GeV/c^) 
(dark: background and light: signal), for 30 fb^^, = 1 TeV 
and f = 0. 



Fig. 11: Reconstructed mass of the radion (m^ — 600 GeV/c^) 
(dark: background and light: signal), for 30 fb^^, A^ = 1 TeV 
and ^ = 0. 



The expected number of events for an integrated luminosity of 30 fb (three years at low luminosity) 
are given in Table 5 for the two radion masses and for the backgrounds, when ^ = and = 1 TeV. 

Requiring a minimum of 10 events and a S/Vb > 5, the maximum reach in is 1.05 and 1.4 TeV 





m0=3OO GeV/c^ 


m<^=600 GeV/c^ 


Signal 


43 


22 


tt 


58 


~ 6 


Z ^ TT 


~ 


~ 


w 


~ 


~ 




5.6 


9.0 



Table 5: Expected number of events for signal and background, for an integrated luminosity of 30 fb ^ for =300 and 600 GeV/c^, 
^ = and A^ = 1 TeV, after all cuts. 
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for m^=300 GeV/c^ and m0=6OO GeV/c^ respectively, but the uncertainties in background subtraction may 
affect considerably the observability of this channel in the first case. 

5. Conclusion 

We have studied the possibility of observing the radion using the ATLAS detector at the LHC. The radion has 
couplings similar to those of the SM Higgs, and mixes with it, but it has also a large effective coupling to 
gluons. Re-interpreting results of previous studies on the search for a SM Higgs in ATLAS, a significance for 
observing a radion decaying into 77 or ZZ^*'^ has been determined as a function of its mass (see Figure 
For an integrated luminosity of 100 fb~^, the values S/^/B (0.1) are obtained for the 77 channel, with 
a mixing parameter ^=0 and a scale K^=\ (lO)TeV, in the range 80 GeV/c^ < < 160 GeV/c^. For the 
ZZ(*) channel, S/V^ ~100 (1) for 200 GeV/c^ < < 600 GeV/c^ for the same conditions. Because the 
couplings are similar to those of the SM Higgs, a good measurement of the branching ratios will be necessary 
to discriminate between the two scalars. 

The radion can also decay into a pair of Higgs scalars, if the masses permit. Two cases were examined: 
(j) ^ hh ^ 'yjbb and (p ^ hh ^ Trbb, for radion masses of 300 and 600 GeV/c^, for nih = 125 GeV/c^ 
and for an integrated luminosity of 30 fb~^. Limits on the maximal reach in were obtained for these two 
channels. For the 7766 channel, the background is negligible and the reach in is 4.6 (5.7) TeV for = 300 
(600) GeV/c^, when ,^=0. For the Trbb channel, the similarity between the signal and background shapes make 
the observation of a radion of mass m0=3OO GeV/c^ difficult, and the reach for A^ is about 1.4 TeV for 
m<^=600 GeV/c2, when ^=0. 
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Part XVI 

Radion Mixing Effects on the Properties of the 
Standard Model Higgs Boson 

J.L.Hewett and T.G. Rizzo 

Abstract 

We examine how mixing between the Standard Model(SM) Higgs boson, h, and the 
radion of the Randall-Sundrum model modifies the expected properties of the Higgs 
boson. In particular we demonstrate that the total and partial decay widths of the 
Higgs, as well as the h gg branching fraction, can be substantially altered from 
their SM expectations, while the remaining branching fractions are modified less than 
<^ 5% for most of the parameter region. 



The Randall-Sundrum(RS) model [Q] offers a potential solution to the hierarchy problem that can be 
tested at present and future accelerators (for an overview of RS phenomenology, see [[13, 16, 283 1). In this 
model the SM fields lie on one of two branes that are embedded in 5-dimensional AdS space described by the 
metric ds"^ = e^^^l^'r/^j^dx'^dx'^ — dy"^, where k is the 5-d curvature parameter of order the Planck scale. To 
solve the hierarchy problem the separation between the two branes, Tc, must have a value of krc ~ 11 — 12. That 
this quantity can be stabilized and made natural has been demonstrated by a number of authors [|lS, 284-288] 
and leads directly to the existence of a radion (r), which corresponds to a quantum excitation of the brane 
separation. It can be shown that the radion couples to the trace of the stress-energy tensor with a strength A 
of order the TeV scale, i.e., Ceff = -r T^/A. (Note that A = \/3A^ in the notation of Ref. |[T3|,p^|283|].) 



This leads to gauge and matter couplings that are qualitatively similar to those of the SM Higgs boson. The 
radion mass (m, ) is expected to be significantly below the scale A implying that the radion may be the lightest 
new field predicted by the RS model. One may expect on general grounds that this mass should lie in the range 
of a few X 10 GeV < rrir < A. The phenomenology of the RS radion has been examined by a number of 



authors [276-278, 289-295] and in particular has been reviewed by Kribs [279]. 

On general grounds of covariance, the radion may mix with the SM Higgs field on the TeV brane through 
an interaction term of the form 



(1) 



where H is the Higgs doublet field, R^^^ [g^] is the Ricci scalar constructed out of the induced metric g^ on the 
SM brane, and ^ is a mixing parameter assumed to be of order unity and with unknown sign. The above action 
induces kinetic mixing between the 'weak eigenstate' tq and /iq fields which can be removed through a set of 
field redefinitions and rotations. Clearly, since the radion and Higgs boson couplings to other SM fields differ 
this mixing will induce modifications in the usual SM expectations for the Higgs decay widths. To make unique 
predictions in this scenario we need to specify four parameters: the masses of the physical Higgs and radion 
fields, mh^r, the mixing parameter ^ and the ratio v/ K, where v is the vacuum expectation value of the SM 
Higgs ~ 246 GeV. Clearly the ratio v/ K cannot be too large as A^ is already bounded from below by collider 
and electroweak precision data [p^, |l^ , |283| ] ; for definiteness we will take v/A < 0.2 and — 1 < ^ < 1 in 
what follows although larger absolute values of ^ have been entertained in the literature. The values of the two 
physical masses themselves are not arbitrary. When we require that the weak basis mass-squared parameters of 
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the radion and Higgs fields be real, as is required by hermiticity, we obtain an additional constraint on the ratio 
of the physical radion and Higgs masses which only depends on the product |^| Exphcitly one finds that either 

^ > 1 + 2 sin^ p + 2\ sin p\ \J\ + sin^ p or ^ < 1 + 2 sin^ p - 2| sin p| y^l + sin^ p where p = tan~^(6^^). 
This disfavors the radion having a mass too close to that of the Higgs when there is significant mixing; the 
resulting excluded region is shown in Fig. [ij These constraints are somewhat restrictive; if we take m/j = 115 
GeV and i\ = 0.1(0.2) we find that either m,. > 189(234) GeV or m,. < 70(56) GeV. This lower mass range 
may be disfavored by direct searches. 

Following the notation of Giudice et al. [ ]296| ], the coupling of the physical Higgs to the SM fermions 
and massive gauge bosons V = W, Z is now given by 

— 1 - V 
L = — (mfff — myVuV '^)[cos pcos 9 + — (sin0 — sin/? cos 6)]h , (2) 
V A 

where the angle p is given above and 9 can be calculated in terms of the parameters ^ and v /A and the physical 
Higgs and radion masses. Denoting the combinations a = cospcos^ and (3 = sin 9 — sinpcos^, the cor- 
responding Higgs coupling to gluons can be written as Cg^G^uG^'^h with Cg = -^[{a + jiP)Fg — 2h^(5j^ 
where 63 = 7 is the 5C/(3) /3-function and Fg is a well-known kinematic function of the ratio of masses of 
the top quark to the physical Higgs. Similarly the physical Higgs couplings to two photons is now given by 
c^^F^^Ft'^h where = ^[(62 + 6y)/5^ - (a + ll3)F^], where 62 = 19/6 and 5y = -41/6 are the 
SU{2) X U{\) /3-functions and F^ is another well-known kinematic function of the ratios of the W and top 
masses to the physical Higgs mass. (Note that in the simultaneous limits a — > 1, /3 ^ we recover the 
usual SM results.) From these expressions we can now compute the change of the various decay widths and 
branching fractions of the SM Higgs due to mixing with the radion. 

Fig. ^ shows the ratio of the various Higgs widths in comparison to their SM expectations as functions of 
the parameter ^ assuming that = 125 GeV with different values of ?ti,. and ^. We see several features right 
away: {i) the shifts in the widths to ff/VV and 77 final states are very similar; this is due to the relatively 
large magnitude of Fj while the combination 62 + W is rather small, (ii) On the other hand the shift for the 
gg final state is quite different since Fg is smaller than Fy and 63 is quite large, (iii) For relatively light radions 
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Fig. 2: Ratio of Higgs widthis to their SM values, _Rr, as a function of ^ assuming a pliysical Higgs mass of 125 GeV: red for fermion 
pairs or massive gauge boson pairs, green for gluons and blue for photons. In the left panel we assume rrir = 300 GeV and w/A = 0.2. 
In the right panel the solid(dashed) curves are for = 500(300) GeV and v/ A. = 0.2(0.1). 

with a low value of A the width into the gg final state can come close to vanishing due to a strong destructive 
interference between the two contributions to the amplitude for values of ^ near - 1 . {iv) Increasing the value of 
rrir has less of an effect on the width shifts than does a decrease in the ratio ^ . 

The deviation from the SM expectations for the various branching fractions, as well as the total width, 
of the Higgs are displayed in Fig. |3| as a function of the mixing parameter ^. We see that the gluon branching 
fraction and the total width may be drastically different than that of the SM. The former will affect the Higgs 
production cross section at the LHC. However, the 77, //, and VV, where V = W, Z branching fractions 
receive small corrections to their SM values, of order <, 5% for most of the parameter region. Observation of 
these shifts will require the accurate determination of the Higgs branching fractions available at an e+e^ Linear 
ColUder. 

In summary, we see that Higgs-radion mixing, which is present in some extra dimensional scenarios, 
can have a substantial effect on the properties of the Higgs boson. These modifications affect the widths and 
branching fractions of Higgs decay into various final states, which in turn can alter the Higgs production cross 
section at the LHC and may require the precision of a Linear Collider to detect. 
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Fig. 3: The deviation from the SM expectations for the Higgs branching fraction into 77, gg, ff, and VV final states as labeled, 
as well as for the total width. The black, red, and blue curves correspond to the parameter choices rrir — 300, 500, 300 GeV with 
v/S, = 0.2, 0.2, 0.1, respectively. 
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Part XVII 

Probing Universal Extra Dimensions at Present and 
Future Colliders 

Thomas G. Rizzo 

Abstract 

In the Universal Extra Dimensions model of Appelquist, Cheng and Dobrescu, all of 
the Standard Model fields are placed in the bulk and thus have Kaluza- Klein (KK) ex- 
citations. These KK states can only be pair produced at colliders due to the tree-level 
conservation of KK number, with the lightest of them being stable and possibly hav- 
ing a mass as low as ~ 350 — 400 GeV. We investigate the production cross sections 
and signatures for these particles at both hadron and lepton colliders. We demonstrate 
that these signatures critically depend upon whether the lightest KK states remain sta- 
ble or are allowed to decay by any of a number of new physics mechanisms. These 
mechanisms which induce KK decays are studied in detail. 



1. Introduction 

The possibility that the gauge bosons of the Standard Model (SM) may be sensitive to the existence of extra 



dimensions near the TeV scale has been known for some time [ |20| , [297| - |302| ]. However, one finds that the 
phenomenology of these models is particularly sensitive to the manner in which the SM fermions (and Higgs 
bosons) are treated. 

Perhaps the most democratic possibility requires all of the SM fields to propagate in the ~ TeV^^ bulk 



[|22|], i.e., Universal Extra Dimensions (UED). In this case, the conservation of momentum in the extra dimen- 
sions is restored and one obtains interactions in the 4-d Lagrangian which take the form ~ gCijkf^^^ 1 ^if'^^^ C^j^y 
which for flat space metrics vanishes unless i + j + A; = 0, as a result of the afore mentioned momentum con- 
servation. Although this momentum conservation is actually broken by orbifolding, one finds, at tree level, 
that KK number remains a conserved quantity. (As we will discuss below this conservation law is itself further 
broken at one loop order.) This implies that pairs of zero-mode fermions, which we identify with those of 
the SM, cannot directly interact singly with any of the excited modes in the gauge boson KK towers. Such 
a situation clearly limits any constraints arising from precision measurements since zero mode fermion fields 
can only interact with pairs of tower gauge boson fields. In addition, at colliders it now follows that KK states 
must be pair produced, thus significantly reducing the possible direct search reaches for these states. In fact. 



employing constraints from current experimental data, Appelquist, Cheng and Dobrescu (ACD) [22] find that 
the KK states in this scenario can be as light as ~ 350 — 400 GeV. If these states are, in fact, nearby, they will 
be copiously produced at the LHC, and possibly also at the Tevatron, in a variety of different channels. It is 
the purpose of this paper to estimate the production rates for pairs of these particles in various channels and to 
discuss their possible production signatures. This is made somewhat difficult by the apparent conservation of 
KK number which appears to forbid the decay of heavier excitations into lighter ones. 

Now although the KK number is conserved at the tree level it becomes apparent that it is no longer so 



at loop order [303]. Consider a self-energy diagram with a field that has KK number of 2n(2n + 1) entering 
and a zero(n = 1) mode leaving the graph; KK number conservation clearly does not forbid such an amplitude 
and constrains the 2 particles in the intermediate state to both have KK number n{n and n + 1). The existence 
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of such amplitudes implies that all even and odd KK states mix separately so that the even KK excitations can 
clearly decay to zero modes while odd KK states can now decay down to the KK number=l state. Thus it is KK 
parity, (—1)", which remains conserved while KK number itself is broken at one loop. Since the lightest KK 
excited states with n = 1 have odd KK parity they remain stable unless new physics is introduced. As we ai^e 
only concerned with the production of pairs of the lightest KK particles in our discussion below, we are faced 
with the possibility of producing heavy stable states at colliders. 

2. Collider Production 

Due to the conservation of KK parity, the first KK excitations of the SM fields must be pair-produced at 
colliders. At 77 and lepton colliders the production cross sections for all the kinematically accessible KK 
states will very roughly be of order 100 fb (1 TeF/\/i)^ which yields respectable event rates for luminosities 
in the 100 — 500 fb^^ range. A sample of relevant cross sections at both 77 and e"'"e~ colliders are shown in 
Fig. H} In the case of 77 collisions we have chosen the process 77 — > W^^^^W~^^^ as it the process which 
has the largest cross section for the production of the first KK state. Similarly, gauge boson pair production in 
e'^e~ collisions naturally leads to a large cross section. Clearly, such states once produced would not be easily 
missed for masses up to close the kinematic limit of the machine independently of how they decayed or if they 
were stable. To directly probe heavier masses we must turn to hadron colliders. 
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Fig. 1: Cross section for 77 — > K^'^'+K^'^' (top panel) for different electron and laser polarizations for -^5^^ = 1 TeV. Cross section 
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Fig. 2: Cross section for the pair production of the lightest colored KK states at the ^/s = 2 TeV Tevatron (left) and the LHC (right). 
In the left panel, from top to bottom on the left-hand side, the curves correspond to the processes ii, v, Hi, i and iv, respectively. In the 
right panel, from top to bottom on the left-hand side, the curves correspond to the processes ii, i, iv. Hi and v, respectively. Antiquark 
contributions are included in reactions ii and iv. 



Since both QCD and electroweak exchanges can lead to KK pair production at hadron colliders there are 
three classes of basic processes to consider. Clearly the states with color quantum numbers will have the largest 
cross sections at hadron machines and there are a number of processes which can contribute to their production 
at order []304|] several of which we list below: 



(0 99 - 




(ii) qq' - 




in) gg + qq - 




(iv) qq - 




(v) qq - 





(1) 

where the primes are present to denote flavor differences. Fig. ^ shows the cross sections for these five 
processes at both the ^/s = 2 TeV Tevatron and the LHC summed over flavors. It is clear that during the 
Tevatron Run II we should expect a reasonable yield of these KK particles for masses below ~ 600 GeV if 
integrated luminosities in the range of 10-20 fb~^ are obtained. Other processes that we have not considered 
may be able to slightly increase this reach. For larger masses we must turn to the LHC where we see that 
significant event rates should be obtainable for KK masses up to ~ 3 TeV or so for an integrated luminosity of 
100 / 6^^. As one might expect we see that the most important QCD processes for the production of KK states 
are different at the two colliders. 

The real signature of the UED scenario is that all of the SM fields have KK excitations. Thus we will 
also want to observe the production of the SM color singlet states. Of course color singlet states can also be 
produced, with the largest cross sections being for associated production with a colored state at order aa^; these 
rates are of course smaller than for pairs of colored particles as can be seen in Fig. ^ Here we see reasonable 
rates are obtained for KK masses in excess of ~ 1.8 TeV or so. Lastly, it is possible to pair produce color 
singlets via electroweak interactions which thus lead to cross sections of order a^. Due to the large center of 
mass energy of the LHC these cross sections can also lead to respectable production rates for KK masses as 
great as ~ 1.5 TeV as can be seen from Fig. 0. It is clear from this analysis that the LHC will have a significant 
search reach for both colored and non-colored KK states provided that the production signatures are reasonably 
distinct. 
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Fig. 3: Cross sections for the associated production of the lightest color singlet KK states at the LHC: in the left panel, from top to 
bottom, for W^^^^^, g'^^Z'^' and s^^'t'^' final states; in the right panel, from top to bottom, for g<^>VF<^>*, q^^'Z*^' and g*^^*^' 
final states. Anti-quark contributions are included. 
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Fig. 4: Cross sections for the production of the lightest color singlet KK states at the LHC: in the left panel, from top to bottom, for 
2Zi\7(i'z(^) and 27'^' final states; in the right panel, from top to bottom, for Wf^'+W'^'", W'^'^Z'^) and M/(^'=^7'^' final states. 



3. Signatures 

When examining collider signatures for KK pair production in the UED there are two important questions to 
ask: (i) Are the lightest KK modes stable and (ii) if they are unstable what are their decay modes? From the 
discussion above it is clear that without introducing any new physics the n = 1 KK states are stable so we must 
consider this possibility when looking at production signatures. 



In their paper ACD [ 22 1 argue that cosmological constraints possibly suggest that KK states in the TeV 
mass range must be unstable on cosmological time scales. (Of course this does not mean that they would appear 
unstable on the time scale of a collider experiment in which case our discussion is the same as that above.) This 
would require the introduction of new physics beyond that contained in the original UED model. There aie 
several possible scenarios for such new physics. Here we will discuss three possibilities in what follows, the 
first two of which were briefly mentioned by ACD [22|. 

Scenario I: The TeV~^-scale UED model is embedded inside a thick brane in a higher ((^+4)-dimensional 
space, with a compactification scale Rg » i?o in which gravity is allowed to propagate in a manner similar 
to the model of Arkani-Hamed, Dimopoulos and Dvali [|3[ [ij, |l^]. Since the graviton wave functions are 
normalized on a torus of volume {2ttRgY while the KK states are normalized over 2ttRc the overlap of a 
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Fig. 5: Width for the decay of the first excited KK state (even-top left panel, odd-top right panel) into the corresponding zero mode and 
a graviton tower as a function of the mass of the KK state. The solid(dashed) lines are for Md = 5(10) TeV and from top to bottom 
in each case the curves correspond to 5 = 2, 4, 6, respectively. The lower panel shows the missing energy distribution for these decays 
for these same cases assuming a KK mass of 1 TeV. 



KK zero mode with any even or odd KK tower state n and a graviton will be non-zero. In a sense, the brane 



develops a transition form-factor analogous to that described in p05| ]. This induces transitions of the form 
KK{n = 1) ^ KK{n = 0) + G where G represents the graviton field which appears as missing energy in the 
collider detector. This means that production of a pair of KK excitations of, e.g. , quarks or gluons would appear 
as two jets plus missing energy in the detector; the corresponding production of a KK excited pair of gauge 
bosons would appear as the pair production of the corresponding zero modes together with missing energy. We 
can express this form-factor simply as 



TtRc 



dye*'^»^(cos ny/Rc, sin ny/Rc 



(2) 



for even and odd KK states, respectively, where nig is the graviton mass. Here we have assumed that the thick 
brane resides at = for all i / 1. Given these form-factors we can calculate the actual decay rate for 
KK{n = 1) ^ KK{n = 0) + G, where we now must sum up the graviton towers by following the analyses in 
Ref. [ ]296| , [306| ]; this result should be relatively independent of the spin of the original KK state. Performing the 
necessary integrations numerically we obtain the results shown in Fig. ^. This figure shows that this mechanism 
provides for a very rapid decay over almost all of the parameter space. For light KK states with both 6 and Md 
large the decay rate is suppressed and may lead to finite length charged tracks in the detector. (In particular 
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Table 1: Individual branching fractions in per cent for the first excited fermion KK modes when KK level mixing occurs as in Scenario 

n. 



the production of a charged KK state with a long lifetime would yield a kink-like track structure due to the 
decay to the graviton tower.) Although not a true two-body decay, Fig. || also shows that the typical missing 
energy in the gravitational decay of a KK state will be close to half its mass, which is quite significant for such 
heavy objects. It is clear that events with such a large fraction of missing energy should be observable above 
background given sufficient event rates. These events will not be confused with SUSY signals since they occur 
in every possible channel. 

Scenario II: KK decays can be induced in the UED model by adding a 'benign' brane at some y = yo 
which induces new interactions. By 'benign' we mean that these new interactions only do what we need them 
to do and do not alter the basic properties of the UED model. The simplest form of such interactions are just 
the four dimensional variants of the terms in the the 5-d UED action. For example, one might add a term such 
as 

d^x I dy 6{y - yo) ^V^t^^aV' , (3) 

where A is some Yukawa-like coupling and Mg is some large scale. Note that the brane is placed at some 
arbitrary position y = yo and not at the fixed points where only even KK modes would be effected. These new 
interactions result in a mixing of all KK states both even and odd and, in particular, with the zero mode. Thus 
we end up inducing decays of the form KK^^^ KK^^^ KK^'^). For KK fermions the decay into a fermion plus 
gauge boson zero mode is found to be given by 

r(/(i) ^ /(O)y(o)) = ||,2^j^ .N^.PS, (4) 

where s,^ is the induced mixing angle, N^. is a color factor, gy, the relevant gauge coupling and PS is the 
phase space for the decay. It is assumed that the mixing angle is sufficiently small that single production of 
KK states at colliders remains highly suppressed but is large enough for the KK state to decay in the detector. 
For A ~ 0.1 and Mg ~ a few Mc this level of suppression is quite natural. (Numerically, it is clear that 
the KK state will decay inside the detector unless the mixing angle is very highly suppressed.) The resulting 
branching fractions can be found in Table 1 where we see numbers that are not too different than those for 
excited fermions in composite models with similar decay signatures. However, unlike excited SM fields, single 
production modes are highly suppressed. For KK excitations of the gauge bosons, their branching fractions 
into zero mode fermions will be identical to those of the corresponding SM fields apart from corrections due to 
phase space, i.e., the first excited Z KK state can decay to while the SM Z cannot. 

Scenario III: We can add a common bulk mass term to the fermion action, i.e., a term of the form 
m{DD + SS), where D and S represent the 5-dimensional isodoublet and isosinglet SM fields. We chose 
a common mass both for simplicity and to avoid potentially dangerous flavor changing neutral currents. The 
largest influence of this new term is to modify the zero mode fermion wavefunction which is now no longer 
flat and takes the form ~ e""^'^' and thus remains Z2-even. Clearly there is now a significant overlap in 
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Fig. 6: The form factors C7(upper curve) and Q'(lower curve) as discussed in the text for the case n = 1. 

the 5-d wavefunctions between pairs of fermion zero modes and any Z2-even gauge KK mode which can 
be represented as another form-factor, Q{x), where x = mRc. This form factor then describes the decay 
Qin) _^ /(0)/(0) where G represents a generic KK gauge field. Similarly we can obtain a form-factor that 
describes the corresponding decay j(o) which we denote by G'{x) where x is as above. It is clear 

that the decays of KK states in this scenario will be essentially identical to Scenario II above although they are 
generated by a completely different kind of physics. Fig. ^ shows the shape of these two form-factors as a 
function of the parameter x. The natural question to ask at this point is 'what is the value of m relative to Mc?'. 
It seems natural to imagine that the bulk mass would be of order the compactification scale, the only natural 
scale in the action, which would imply that x ~ 1 so that large form-factors would be obtained. While this 
scenario works extremely well for the decay of Z2-even states it does not work at all for the case of Z2-0AA 
states. 

4. Summary and Conclusions 

In this paper we have begun a detailed examination of the predictions of the Universal Extra Dimensions 
model for future colliders. Since indirect searches for such states give rather poor reaches, direct searches are 
of greater importance in this model. To obtain interesting search reaches requires a hadron collider such as 
the Tevatron or LHC. Based on counting events we expect the reach at the Tevatron Run II (LHC) for KK 
states to be ~ 600(3000) GeV. Within the UED model itself these lightest KK states are stable even when loop 
corrections are included unless new interactions are introduced from elsewhere. If these states are indeed stable, 
the production of a large number of heavy stable charged particles would not be missed at either collider. It is 
more likely, however, that new physics does indeed enter, rendering the KK modes unstable. In this paper we 
have examined three new physics scenarios that induce finite KK lifetimes and compared their decay signatures. 
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Part XVIII 

Black Hole Production at Future Colliders 

S. Dimopoulos and G. Landsberg 

Abstract 

If the scale of quantum gravity is near a TeV, the CERN Large Hadron Collider will 
be producing one black hole (BH) about every second. The decays of the BHs into 
the final states with prompt, hard photons, electrons, or muons provide a clean signa- 
ture with low background. The correlation between the BH mass and its temperature, 
deduced from the energy spectrum of the decay products, can test Hawkings evapora- 
tion law and determine the number of large new dimensions and the scale of quantum 
gravity. We also consider BH production at the proposed future high-energy collid- 
ers, such as CLIC and VLHC, and describe the Monte Carlo event generator that can 
be used to study BH production and decay. 



1. Introduction 



An exciting consequence of TeV-scale quantum gravity |^, 11, 12] is the possibility of production of black 
holes (BHs) [23, 24, 307-309] at the LHC and beyond. This paper summarizes and extends our pioneer work 



on this subject [ 24 ] to the post-LHC future and discusses additional aspects of black-hole phenomenology left 
out from | |24l ] due to lack of space. Since this work has been completed, numerous follow-up publications on 
this exciting subject have appeared in the archives, focusing on both the collider [26, 3^ 4^, 310| ] and cosmic 
ray [35, 38,|9] production. We hope that this new branch of phenomenology of extra dimensions will flourish 
in the months to come, as black hole production might be the very first evidence for the existence of large extra 
dimensions. 

Black holes are well understood general-relativistic objects when their mass Mbh far exceeds the fun- 
damental (higher dimensional) Planck mass Mp '^TeV. As A/bh approaches Mp, the BHs become "stringy" 
and their properties complex. In what follows, we will ignore this obstacle and estimate the properties of light 
BHs by simple semiclassical arguments, strictly valid for Mbh ^ Mp. We expect that this will be an adequate 
approximation, since the important experimental signatures rely on two simple qualitative properties: (i) the 
absence of small couplings and (ii) the "democratic" nature of BH decays, both of which may survive as aver- 
age properties of the light descendants of BHs. Nevertheless, because of the unknown stringy corrections, our 
results are approximate estimates. For this reason, we will not attempt selective partial improvements - such 
as time dependence, angular momentum, charge, hair, and other higher-order general relativistic refinements - 
which, for light BHs, may be masked by larger unknown stringy effects. We will focus on the production and 
sudden decay of Schwarzschild black holes. 



2. Production 



The Schwarzschild radius Rs of an (4 + n) -dimensional black hole is given by [ |31l| ], assuming that extra 
dimensions are large (S> Rs)- 

Consider two partons with the center-of-mass (c.o.m.) energy \/j = Mbh moving in opposite di- 
rections. Semiclassical reasoning suggests that if the impact parameter is less than the (higher dimensional) 
Schwarzschild radius, a BH with the mass Mbh forms. Therefore the total cross section can be estimated from 
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geometrical arguments |312|, and is of order 



Ml 



Mbh /8r(^^)' 



Mp 



n + 2 



2 

n+1 



(seeFig.[T]a,d) 

This expression contains no small coupling constants; if the parton c.o.m. energy \/I reaches the fun- 
damental Planck scale Mp TeV then the cross section is of order TeV~^ « 400 pb. At the LHC or VLHC, 
with the total c.o.m. energy ^/s = 14 TeV or 100-200 TeV, respectively, BHs will be produced copiously. To 
calculate total production cross section, we need to take into account that only a fraction of the total c.o.m. en- 
ergy in a pp collision is achieved in a parton-parton scattering. We compute the full particle level cross section 



using the parton luminosity approach (see, e.g., Ref. [314]) 



da{pp ^ BH + X) (IL 



-a{ah BH) 



dMBH dMBH 

where the parton luminosity dL/dM-Qn is defined as the sum over all the initial parton types: 

dL 2Mbh 



dMBH 



SXn 



and fi{xi) are the parton distribution functions (PDFs). We used the MRSD— ' [315] PDF set with the Q scale 



taken to be equal to Mbh [316], which is within the allowed range for this PDF set, up to the VLHC kinematic 
limit. Cross section dependence on the choice of PDF is 10%. 

The differential cross sections da /c?Mbh for the BH produced at the LHC and a 200 TeV VLHC ma- 
chines are shown in Figs. [T]b and [|f, respectively, for several choices of Mp. The total production cross section 
at the LHC for BH masses above Mp ranges from 0.5 nb for Mp = 2 TeV, n = 7 to 120 fb for Mp = 6 TeV 
and n = 3. If the fundamental Planck scale is 1 TeV, the LHC, with the peak luminosity of 30 fb^^/year 
will produce over lO'^ black holes per year. This is comparable to the total number of Z's produced at LEP, and 
suggests that we may do high precision studies of TeV BH physics, as long as the backgrounds are kept small. 
At the VLHC, BHs will be produced copiously for their masses and the value of the fundamental Planck scale 
as high as 25 TeV. The total production cross section is of the order of a millibarn for Mp = 1 TeV and of 
order a picobam for M p = 25 TeV. 

Similarly, the black holes can be produced at future high-energy lepton colliders, such as CLIC or a muon 
collider. To a first approximation, such a machine produces black holes of a fixed mass, equal to the energy of 
the machine. The total cross section of such a BH produced at a 3 TeV and a 5 TeV machine, as a function of 
Mp and n, is shown in Fig |^ and Fig. |2|f, respectively. For more elaborated studies of the BH production at 
electron colliders, one should take into account machine beamstrahlung. The beamstrahlung-corrected energy 
spectrum of the machine plays the same role as the parton luminosity at a hadron collider, except that for the 
e"'"e~ machine it is peaked at the nominal machine energy, rather than at small values of characteristic of 
a hadron collider. Using typical beamsstrahlung spectra expected for a 3 TeV or a 5 TeV CLIC machine, we 
show the differential cross section da / dM^n of the black hole production at a 3 and a 5 TeV CLIC in Fig. ^ 
and Fig. B, respectively. 




Fig. 1: Black-hole properties at the LHC a)-d),h) and VLHC d)-h). a,e) Parton-level production cross section; b,f) differential cross 
section da / dAfeH ; c,g) Hawking temperature; d) average decay multiplicity for a Schwarzschild black hole; and h) black-hole lifetime. 
The number of extra spatial dimensions n = 4 is used for a)-c), e)-h). The dependence of the cross section and the Hawking temperature 
on n is weak and would be hardly noticeable on the logarithmic scale. The lifetime drops by about two orders of magnitude for n 
increase from 2 to 7. 

3. Decay 

The decay of the BH is governed by its Hawking temperature Th, which is proportional to the inverse radius, 
and given by [ ]31l[ ] : 

Tn = MJ^^l±iy^'^-±l = ^ (1) 

(see Figs. |I]c,g andp5,g). As the collision energy increases, the resulting BH gets heavier and its decay products 
get colder. 

Note that the wavelength A = corresponding to the Hawking temperature is larger than the size of 
the black hole. Therefore, the BH acts as a point-radiator and emits mostly s-waves. This indicates that it 
decays equally to a particle on the brane and in the bulk, since it is only sensitive to the radial coordinate and 
does not make use of the extra angular modes available in the bulk. Since there are many more particles on 
our brane than in the bulk, this has the crucial consequence that the BH decays visibly to standard model (SM) 



particles [gOg^ 

The average multiplicity of particles produced in the process of BH evaporation is given by: {N) = 
/Mbr\ where E is the energy spectrum of the decay products. In order to find (A^), we note that the BH 



\ ^ 

evaporation is a blackbody radiation process, with the energy flux per unit of time given by Planck's formula 



^ ~ i^f^' where x = E/Th, and c is a constant, which depends on the quantum statistics of the decay 



products (c = — 1 for bosons, +1 for fermions, and for Boltzmann statistics). 




Fig. 2: Black hole properties at high-energy lepton colliders. Plots a)-d) and f)-i) correspond to the properties (production cross-section, 
temperature, lifetime, and average decay multiplicity) of a fixed-mass 3 TeV and 5 TeV black hole produced at a 3 TeV or a 5 TeV 
machine, respectively. Plots e),j) show the differential cross section of BH production for n = 4, as a function of the BH mass at a 3 
TeV or a 5 TeV CLIC e^e~ -collider, respectively. 



The spectrum of the BH decay products in the massless particle approximation is given by: ^ ~ 
~ For averaging the multiplicity, we use the average of the distribution in the inverse particle 

energy: 

/ 1 \ 1 r°°(ixi— 



E / Th Jq°° dx 



where a is a dimensionless constant that depends on the type of produced particles and numerically equals 
0.68 for bosons, 0.46 for fermions, and ^ for Boltzmann statistics. Since a mixture of fermions and bosons 
is produced in the BH decay, we can approximate the average by using Boltzmann statistics, which gives the 
following formula for the average multiplicity: (N) ^2T,^ ■ Using Eq. (|l|) for Hawking temperature, we 
obtain: 
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Eq. (Q) holds for Mbh ^ Th, i.e. {N) » 1; otherwise, the Planck spectrum is truncated at E ^ Mbh/2 
by the decay kinematics [ |318 |. The average number of particles produced in the process of BH evaporation is 



shown in Figs. |T]d and ^,i. 

The lifetime of the BH can be estimated by using the Stefan's law of thermal radiation. Since BH 
evaporation occurs primarily in three spatial dimensions, the canonical 3-dimensional Stefan's law applies, and 
therefore the power dissipated by the Hawking's radiation per unit area of the BH event horizon is p = crTf^, 
where a is the Stefan-Boltzmann constant and Th is the Hawking temperature of the BH. Since the effective 
evaporation area of the BH is the area of a 3-dimensional sphere with the radius equal to the BH Schwarzschild 
radius Rs, the total power dissipated by the BH is given by: 

P = AttRIp = ATTRlaT^ = aT^-^ ^ 



47r 

The BH lifetime r is given by: 

and using Eq. (|l|), as well as the expression for a in natural units (h = c = k = I), a = 7r^/60 [ ]319| ], we find: 

_ 3840 /Mbh \ ^ / ST (^) ^ ^ 



Mp{n + l)'^\Mp J y n + 2 J 

The lifetime of a black hole as a function of its mass and the fundamental Planck scale is shown in Figs. |l|h 
and |2|3,h. A typical lifetime of a BH is ~ 10^^^ s, which corresponds to a rather narrow width of the BH state 
~ 10 GeV, i.e. typical for, e.g., a W or Z' resonance of a similar mass. 

We emphasize that, throughout this paper, we ignore time evolution: as the BH decays, it gets lighter 
and hotter and its decay accelerates. We adopt the "sudden approximation" in which the BH decays, at its 
original temperature, into its decay products. This approximation should be reliable as the BH spends most 
of its time near its original mass and temperature, because that is when it evolves the slowest; furthermore, 
that is also when it emits the most particles. Later, when we test the Hawking's mass-temperature relation by 
reconstructing Wien's displacement law, we will minimize the sensitivity to the late and hot stages of the BHs 
life by looking at only the soft part of the decay spectrum. Proper treatment of time evolution, for Mbh ~ Mp, 
is difficult, since it immediately takes us to the stringy regime. 

4. Branching Fractions 

The decay of a BH is thermal: it obeys all local conservation laws, but otherwise does not discriminate between 
particle species (of the same mass and spin). Theories with quantum gravity near a TeV must have additional 
symmetries, beyond the standard SU{3) x SU{2) x U{1), to guarantee proton longevity, approximate lepton 
number(s) and flavor conservation [020]. There are many possibilities: discrete or continuous symmetries, four 



dimensional or higher dimensional "bulk" symmetries [ |321| ]. Each of these possible symmetries constrains 
the decays of the black holes. Since the typical decay involves a lai^ge number of particles, we will ignore the 
constraints imposed by the few conservation laws and assume that the BH decays with roughly equal probability 
to all of the w 60 particles of the SM. Since there are six charged leptons and one photon, we expect ~ 10% of 
the particles to be hard, primary leptons and ~ 2% of the particles to be hard photons, each carrying hundreds 
of GeV of energy. This is a very clean signal, with negligible background, as the production of SM leptons 
or photons in high-multiplicity events at the LHC occurs at a much smaller rate than the BH production (see 
Fig. These events are also easy to trigger on, since they contain at least one prompt lepton or photon with 
the energy above 100 GeV, as well as energetic jets. 
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5. Test of Hawking Radiation 

Furthermore, since there are three neutrinos, we expect only ~ 5% average missing transverse energy (^) 
per event, which allows us to precisely estimate the BH mass from the visible decay products. We can also 
reconstruct the BH temperature by fitting the energy spectrum of the decay products to Planck's formula. 
Simultaneous knowledge of the BH mass and its temperature allows for a test of the Hawking radiation and can 
provide evidence that the observed events come from the production of a BH, and not from some other new 
physics. 

There are a few important experimental techniques that we will use to carry out the numerical test. First of 
all, to improve precision of the BH mass reconstruction we will use only the events with ^ consistent with zero. 
Given the small probability for a BH to emit a neutrino or a graviton, total statistics won't suffer appreciably 
from this requirement. Since BH decays have large jet activity, the Mbh resolution will be dominated by the 
jet energy resolution and the initial state radiation effects, and is expected to be ~ 100 GeV for a massive BH. 
Second, we will use only photons and electrons in the final state to reconstruct the Hawking temperature. The 
reason is twofold: final states with energetic electrons and photons have very low background at high and 
the energy resolution for electrons and photons remains excellent even at the highest energies achieved in the 
process of BH evaporation. We do not use muons, as their momenta are determined by the track curvature 
in the magnetic field, and thus the resolution deteriorates fast with the muon momentum growth. We also 
ignore the r-lepton decay modes, as the final states with r's have much higher background than inclusive 
electron or photon final states, and also because their energies can not be reconstructed as well as those for the 
electromagnetic objects. The fraction of electrons and photons among the final state particles is only ~ 5%, 
but the vast amount of BHs produced at the LHC allows us to sacrifice the rest of the statistics to allow for 
a high-precision measurement. (Also, the large number of decay particles enhances the probability to have a 
photon or an electron in the event.) Finally, if the energy of a decay particle approaches the kinematic limit 
for pair production, Mbh/2, the shape of the energy spectrum depends on the details of the BH decay model. 
In order to eliminate this unwanted model dependence, we use only the low part of the energy spectrum with 
E < Mbh/2. 

The experimental procedure is straightforward: we select the BH sample by requiring events with high 
mass (> 1 TeV) and multiplicity of the final state {N > 4), which contain electrons or photons with energy 
> 100 GeV. We smear the energies of the decay products with the resolutions typical of the LHC detectors. 
We bin the events in the invariant mass with the bin size (500 GeV) much wider than the mass resolution. The 
mass spectrum of the BHs produced at the LHC with 100 fb^^ of integrated luminosity is shown in Fig. || for 
several values of Mp and n. Backgrounds from the SM Z{ee)+ jets and 7+ jets production, as estimated with 
PYTHIA [|lT|], are small (see Fig. |). 



To determine the Hawking temperature as a function of the BH mass, we perform a maximum likelihood 
fit of the energy spectrum of electrons and photons in the BH events to Planck's formula (with the coefficient 
c determined by the particle spin), below the kinematic cutoff (A/bh/2). This fit is performed using the entire 
set of the BH events (i.e., not on the event-by-event basis), separately in each of the Mbh bins. We then use 
the measured Mbh vs. Th dependence and Eq. (|l]) to determine the fundamental Planck scale Mp and the 
dimensionality of space n. Note that to determine n we can also take the logarithm of both sides of Eq. (|l|): 

log(rH) = log(MBH) + const, (4) 

n + 1 

where the constant does not depend on the BH mass, but only on Mp and on detailed properties of the bulk 
space, such as the shape of extra dimensions. Therefore, the slope of a straight-line fit to the log(T// ) vs. 
log(MBH) data offers a direct way of determining the dimensionality of space. This is a multidimensional 
analog of Wien's displacement law. Note that Eq. (0) is fundamentally different from other ways of determining 
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Fig. 3: Number of BHs produced at the LHC in the electron or photon decay channels, with 100 fb~^of integrated luminosity, as a 
function of the BH mass. The shaded regions correspond to the variation in the number of events for n between 2 and 7. The dashed 
line shows total SM background (from inclusive Z(ee) and direct photon production). The dotted line corresponds to the Z{ee) + X 
background alone. 
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Fig. 4: Determination of the dimensionality of space via Wien's displacement law at the LHC with 100 fb ^ of data. 



the dimensionality of space-time, e.g. by studying a monojet signature or virtual graviton exchange processes, 
also predicted by theories with large extra dimensions. (The properties of the latter two processes always 
depend on the volume of the extra-dimensional space, i.e. they cannot yield informaton on the number of extra 
dimensions without specific assumptions on their relative size.) 

A test of Wien's law at the LHC would provide a confirmation that the observed e + X and 7 + X event 
excess is due to BH production. It would also be the first experimental test of Hawking's radiation hypothesis. 
Figure ^ shows typical fits to the simulated BH data at the LHC, corresponding to 100 fb^^ of integrated 
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Table 1: Determination of Mp and n from Hawking radiation. The two numbers in each column correspond to fractional uncertainty 
in Mp and absolute uncertainty in n, respectively. 



Mp 


1 TeV 


2 TeV 


3 TeV 


4 TeV 


5 TeV 


n = 2 


1%/0.01 


l%/0.02 


3.3%/0.10 


16%/0.35 


40%/0.46 


n = 3 


1%/0.01 


1.4%/0.06 


7.5%/0.22 


30%/1.0 


48%/1.2 


n = 4 


1%/0.01 


2.3%/0.13 


9.5%/0.34 


35%/1.5 


54%/2.0 


n = 5 


l%/0.02 


3.2%/0.23 


17%/1.1 






n = 6 


l%/0.03 


4.2%/0.34 


23%/2.5 


Fit fails 


n = 7 


l%/0.07 


4.5%/0.40 


24%/3.8 







luminosity, for the highest fundamental Planck scales that still allow for determination of the dimensionality 
of space with reasonable precision. The reach of the LHC for the fundamental Planck scale and the number of 
extra dimensions via Hawking radiation extends to Mp ~ 5 TeV and is summarized in Table |l| | |322 |. 

Similar tests can be performed at the VLHC and CLIC machines. While the VLHC case is identical to 
that at the LHC, with appropriately scaled energies, CLIC is complementary to the LHC in many ways, as the 
maximum number of BH produced at CLIC is found at the highest accessible masses. This has the advantage 
that the stringy effects, as well as the kinematic distortion of the Planck black-body spectrum, decrease with 
the increase of the BH mass. Thus the Mbh vs. Th fit at CLIC is less affected by these unknown effects. 
Preliminary studies show that statistical sensitivity to the number of extra dimensions and the value of the 
fundamental Planck scale at CLIC is similar to that at the LHC. 

Note, that the BH discovery potential at the LHC and VLHC is maximized in the e/ ^ + X channels, 
where background is much smaller than that in the 7 + X channel (see Fig. |3|). The reach of a simple counting 
experiment extends up to Mp w 9 TeV at the LHC and Mp w 50 TeV at the VLHC {n = 2-7), where one 
would expect to see a handful of BH events with negligible background. 



6. Black Hole Monte Carlo Event Generator 

A Monte Carlo package, TRUENOIR, has been developed for simulating production and decay of the black 



holes at high-energy colliders. This package is a plug-in module for the PYTHIA [ |113| ] Monte Carlo generator. 
It uses a heuristic algorithm and conservation of baryon and lepton numbers, as well as the QCD color, to 
simulate the decay of a black hole in a rapid-decay approximation. While the limitations of such a simplistic 
approach are clear, further improvements to this generator are being worked on. In the meantime, it provides a 
useful qualitative tool to study the detector effects and other aspects of the BH event reconstruction. Figure |5| 
shows a display of a typical BH event at a 5 TeV CLIC collider, produced using the TRUENOIR code. A 
characteristic feature of this event is extremely large final state multiplicity, very atypical of the events produced 
in e+e~ collisions. 



7. Summary 

Black hole production at the LHC and beyond may be one of the early signatures of TeV-scale quantum gravity. 
It has three advantages: 

• (i) Large Cross Section: No small dimensionless coupling constants, analogous to a, suppress the pro- 
duction of BHs. This leads to enormous rates. 

• (ii) Hard, Prompt, Charged Leptons and Photons: Thermal decays are flavor-blind. This signature has 
practically vanishing SM background. 
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Fig. 5: A typical black-hole event at a 5 TeV CLIC accelerator. The two views correspond to the end and side-views of a CLIC detector. 
Detector simulation by Albert De Roeck. 

• (iii) Little Missing Energy: This facilitates the determination of the mass and the temperature of the black 
hole, and may lead to a test of Hawking radiation. 

It is desirable to improve our primitive estimates, especially for the light black holes (Mbh ~ Mp); 
this will involve string theory. Nevertheless, the most telling signatures of BH production - large and growing 
cross sections; hard leptons, photons, and jets - emerge from qualitative features that are expected to be reliably 
estimated from the semiclassical arguments of this paper. 

Perhaps black holes will be the first signal of TeV-scale quantum gravity. This depends on, among other 
factors, the relative magnitude of M p and the (smaller) string scale Ms. For Ms <C Mp, the vibrational modes 
of the string may be the first indication of the new physics. 

Studies of the BH properties at future facilities, including very high-energy lepton and hadron colliders, 
would make it possible to map out properties of large extra dimensions, to measure the effects of quantum 
gravity, and to provide insight into other quantum phenomena, such as Hawking radiation, the information 
paradox, etc. 
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